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Abstract. The limiting absorption principle asserts that if is a suitable 
Schrodinger operator, and / lives in a suitable weighted space (namely 
lj0,i/2+a foj, gQjj^g ^ > 0), then the functions R{\ + ie)f := {H - \ - ie)-^ f 
converge in a another weighted L space unique solution u 

of the Helmholtz equation {H — X)u = / which obeys the Sommerfeld outgo- 
ing radiation condition. In this paper, we investigate more quantitative (or 
effective) versions of this principle, for the Schrodinger operator on asymptot- 
ically conic manifolds with short-range potentials, and in particular consider 
estimates of the form 

|-R(A + ie)f\\jjo,-i/2-a < C(A, /i')||/||^o,i/2+<T ■ 
We are particularly interested in the exact nature of the dependence of the 
constants C{X,H) on both A and H. It turns out that the answer to this 
question is quite subtle, with distinctions being made between low energies 
A ^ 1, medium energies A ~ 1, and large energies A 2> 1, and there is also 
a non-trivial distinction between "qualitative" estimates on a single operator 
H (possibly obeying some spectral condition such as non-resonance, or a geo- 
metric condition such as non-trapping), and "quantitative" estimates (which 
hold uniformly for all operators H in a certain class). Using elementary meth- 
ods (integration by parts and ODE techniques), we give some sharp answers 
to these questions. As applications of these estimates, we present a global-in- 
time local smoothing estimate and pointwise decay estimates for the associated 
time-dependent Schrodinger equation, as well as an integrated local energy de- 
cay estimate and pointwise decay estimates for solutions of the corresponding 
wave equation, under some additional assumptions on the operator H. 



1. Introduction 



1.1. The limiting absorption principle for the free Schrodinger operator. 

The limiting absorption principle for the free Schrodinger operator Hq := —A on 
the Euchdean space R" describes the behavior of the resolvents Ro{X ± ie) :~ 
(Hq — (A ± ie))^^ in the limit e ^ in the weighted Sobolev spaces H^'^CR"-), 
defined for s = 0, 1, 2 and m G R by 

s 

UWh^-'^W :=^|1(x->'"|VVIIIl^(R"), 

where (x) := (1 -I- 
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Proposition 1.2 (Limiting absorption principle for Hq). Let X > 0, n > 3, e,a > 
0, and ± be a sign. Then for any f e iJ°'^/^+'^(R") the problem 

{Ho-{\±ie))u = f 

has the unique solution u±e = -Ro(A ± ie)f G if^(R") and we have the estimate^ 

(1) ||-Ro(A ± ie)/||H0,-l/2-<,(R,n) < CA~^/^||/||jjO,l/2+a(R,n) 

for all e > 0, and we have the variant estimate 

(2) ||i?o(A ± «£)/||^2.-x/2-.(R„) < C{1 + X)X-'/^\\f\\H0..f2+,^^^y 

Furthermore, as e ^ 0, 

w±e = i?o(A ± ie)f ^ ii± = i?o(A ± iO)f in H^'-^^^'" {BJ"). 

The function u± is the unique solution in ff^'~^/^~'^(R") to the Helmholtz equation 
{Hq — X)u = f which obeys the Sommerfeld radiation condition 

where 5^ := jfy • V^^ is the radial derivative. Moreover, 

{dr T iVX±iI)u±e {dr T iVX)u± m ijO,-l/2+<T' ^ny 
Finally, if o < \, we may replace (1) with the estimate 

(3) ||i?o(A ± i£)/||^o,-3/.+.(R„) < C(l + X)-''^f\\H0.U2^.^^.y 

A similar version, with the factor o/(l + A)~^/^ instead o/A~^/^, holds in place of 
(2). 

The result is well-known for the free Schrodinger Hamiltonian Hq = — A. It was 
established for a more general class of differential operators with constant coeffi- 
cients by Agmon [1]. We will reprove this result in the course of this paper (see 
in particular Section 6.1 for an elementary argument). The sharpness of the above 
estimates can be easily verified from the behavior of the kernel 



giVAibie \x-y\ 

Ro{X ± ie){x, y) ^ C |^ _ ^|„_2 

of the resolvent Ro{X ± ie) for |a; — ?/| > CA~^/^ and the bound Ro{X ± ie){x,y) < 
C\x - for \x-y\< CX-^I"^ . 

The limiting absorption principle can be viewed as a qiiantitative formulation of 
the fact that the operator Hq has no embedded eigenvalues or resonances in its ab- 
solutely continuous spectrum [0, +00). As is well known, it is also closely connected 
to the limiting amplitude principle for the wave equation, and the local smoothing 
estimate for the Schrodinger equation. We shall return to these connections later 
in this introduction. 



^Here and in the sequel, we allow all absolute constants C to depend on the dimension n and 
on exponents such as a. 
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1.3. Quantitative limiting absorption. In this paper wc study the problem of 
extending the hmiting absorption principle (and its applications) in a quantitative 
manner to other Schrodinger operators H. More specifically, we shall limit our 
attention to operators of the form H = —Am + V, where M is an asymptotically 
conic manifold and ^ is a short-range potential; we now make these concepts more 
precise. Henceforth n > 3 and a,ao > are fixed, and we allow all constants to 
depend on n, a, and ctq. 

Definition 1.4 (Asymptotically conic manifold). An asymptotically conic manifold 
is a smooth connected n-dimensional Riemannian manifold (M, g) which, outside 
of a compact set Kq C M, can be parameterized as M\Kq = [i?o,oo) x dM = 
{(r, w) : r > Ro,u) £ dM} for some Rq > 1, where {dM,h) is a smooth compact 
n— 1-dimensional Riemannian manifold, and the metric g takes the scattering form^ 



for all {r,y) S M\Ko. We also make the qualitative assumption that the function 
eah(l/s,y) extends smoothly to a function on [0, 1/Ro] x dM. We define a weight 
(x) on M by setting {x) := r on M\Kq, and extending the weight smoothly to Kq 
so that it stays between Rq and Ro/2. Wc define the space L?'{M) by using the 
measure dg{x) :— ^/g dx induced by the Riemannian metric, and let Am := V"Va 
be the usual Laplace-Beltrami operator, where Vq, denotes covariant differentiation 
with respect to the Levi-Civita connection, raised and lowered in the usual manner. 
We define the weighted Sobolev spaces iJ*''"(M) for s = 0, 1, 2 and m e R by the 
formula 



Important convention. Henceforth, all constants C are allowed to depend on 
the manifold M and the quantities RQ,C,aQ appearing above, as well as on the 

dimension n. 

Remark 1.5. The class of asymptotically conic manifolds includes the class of 
asymptotically flat manifolds (as studied for instance in [30], [40], [12], [69], [98]) as 
a sub-class, in which dM is the unit sphere S*""^ with the standard metric. How- 
ever in general we do not assume any topological flatness of M, nor do wc assume 
that dM is topologically a sphere. In particular it is certainly possible for M to 
contain trapped geodesies. The condition (6) seems to indicate that we only require 
control on our metric. However, this is because of our decision to use normal 
co-ordinates (4). A metric in a more general long-range form, e.g. g = go + r~'^''°e, 
where g^ = dr'^ ~\- r'^hab and e is a smooth function obeying symbol- type estimates 
of order 0, can be placed in the normal form (4) (sec [52, Section 10.5]) but as 

^We use the indices a, b to parameterize the n — 1-dimensional space dM, and the indices i,j 
to parameterize the n-dimensional space M. 




s 



(7) 




fe=0 
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is well known, the use of normal forms costs one degree of regularity, so that one 
would need type bounds on the metric g in the original co-ordinates in order to 
get the C^-type control in (4). 

We define a short-range potential : M — >■ R on M to be any real- valued function V 
such that {xY^'^'^°V{x) e L°°{M) for some ctq > 0, and then define the Schrddinger 
operator H := — Am + V. 

Remark f .6. The spectral theory of the operator H in the short-range case is well 
understood; indeed, it is known, see e.g. [49], that H is essentially self-adjoint, so 
that the resolvents i?(A ± ie) := {H — (A ± ie))~^ are well-defined for e 7^ and are 
bounded operators on L'^{M) (with bounds depending on £, of course). We note 
that the essential self-adjointness of H, without symbol behavior assumptions on 
the potential V and the metric coefficients of g, follows from a lower bound 

(8) / Huu>-C [ \u\^ 

JM JM 

and uniqueness of sohitions of the wave equation 

utt + Hu^O, u{0) = uo, dtu{0)=ui, 

see [6], [26], [85] and references therein. The lower bound (8) easily follows from 
the L°° bound on the potential V and positivity of the Laplace-Beltrami operator 
—Am- The uniqueness statement is a consequence of the energy identity 

/ {\dMtt + \Vu{t)\l + V\u{t)f)dg= [ {\dMO)f + \^u{0)\l + V\u{0)f)dg. 

JM JM 

Furthermore, the spectrum <t{H) consists of the right half-line [0, -|-oo) (where 

it is piircly absolutely continuous, except possibly at where one may have an 
eigenvalue), together with a finite number of eigenvalues Ai < A2 < . . . < Afe < 
on the negative half-line (—00, 0); the finite number of bound states is implied by the 
Cwikcl-Lieb-Rozcnblum inequality for manifolds, see e.g. [67], [68], since V G _L"/^. 
In particular one can define spectral multipliers m{H) of H for any bounded (or 
polynomially growing) m in the usual manner. 

Our first main result establishes a limiting absorption principle with explicit control 

on the constants in this general setting: 

Theorem 1.7 (Quantitative limiting absorption principle). Let A > 0, n > 3, 
cr, (To > 0, and ± be a sign. Let M be an asymptotically conic manifold, let V : 
M — >■ R 6e a short-range potential obeying the pointwise bounds 

(9) \V{x)\ < A (^{x)-^-^^° + X^/^{x)-^-^^°^ 

for all X e M, and let H := —Am + V with resolvents R{X±is) := {H — {X±ie))~^ . 
Let f €H°'^/^+''{W). 

• (Limiting absorption near infinity) There exists a compact set K = K{M, A) 
such that for s = 0, 1, 2 

(10) ||i?(A ± i£)/||H3,-i/.-.(M\K) < A^/2C(M, A)(A-^(^'^) + 1)||/||ho,i/2+.(m)- 
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• (Global limiting absorption) We have the bound 

(11) 

||i?(A±i£)/||H2.-i/2-.(M) < C{M,A) (a-^(^'^) + e^(^'^)^) ||/||^,o,v2+.(M). 

The constant C{M,A) can be explicitly computed in terms of M, A, n, u, gq. 

Remark 1.8. A key feature of this estimate is that it is quantitaMve (or effective), 
in the sense that the bounds in (10), (11) depend only on the underlying manifold 
M and on the bounds enjoyed by the potential V. In particular, no spectral as- 
sumptions on H (e.g. involving cigcnfunctions or resonances at zero) arc assumed. 
Later on we shall also discuss qualitative (or ineffective) limiting absorption prin- 
ciples, in which the bounds are obtained indirectly (via compactness arguments or 
Fredholm theory) and are allowed to depend on the potential V (and in particular 
on the spectral behavior of V). In order to obtain effective bounds in Theorem 1.7, 
we shall avoid the use of compactness methods, and instead rely on "elementary" 
methods such as integration by parts and ODE analysis. 

Remark 1.9. The condition (9) on the potential demands l/(a;)^+ type decay at 

low energies, but only l/{x)^^ decay at high energies. It appears to be essentially 
optimal in the class of pointwise bounds. It is almost scale invariant under the 
transformation x jjLX, A The limiting absorption principle in the Eu- 

clidean case H = — Ar^ + V was established for A > c > by Agmon [1] essentially 
under the assumption that |l^(a;)| < A{x)~^~'^'^° . The global- in-time local smooth- 
ing estimate and Strichartz estimates for the Schrodinger group e'*^~^R^+^) were 
proved in [5] and [45] respectively under the assumption that |V^(.t)| < A(x)"'^^^'^" . 
These results are qualitative as their proof relies on Fredholm theory. In [19] (see 
also the followup papers [90], [48]), these results were extended to the class of inverse 
square potentials V{x) = A/jxj^; due to the explicit nature of the cigcnfunctions in 
this case, Fredholm theory can be avoided, and the results here are quantitative. 

It would be natural to replace pointwise conditions of (9) corresponding integral 
assumptions on V , i.e., on the scale- invariant L"/^-norm of V. In the flat case there 
is a number of recent qualitative results, see [44] , [58] . 

The second term in (9) suggests that one could weaken the decay hypotheses on 
V in exchange for more regularity. It is well known however that the limiting 
absorption principle can fail for potentials decaying like l/|.x| or slower. As the 
standard example with the Wigner-Von-Neumann potential shows, such potentials 
can create embedded eigenvalues in the continuous spectrum, which would clearly 
destroy the limiting absorption principle in this case. 

Remark 1.10. The bound (10) shows that (at least for high energy A) geometry of 
the "black box" compact region K essentially does not affect the behavior of the 
resolvent restricted to the complement of K. Such bounds for sufficiently large A 
with the constant C(M, A, A) < C(M, A)A~^/^ have been proven to hold in [22], see 
also [15]. It is also interesting to note that in this context, local in time Strichartz 
estimates, restricted to the complement of a compact set, for solutions of a time- 
dependent Schrodinger equation have been established in [10]. 
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Remark 1.11. The high-energy case A ^ 1 of the above theorem is essentially 
contained in earlier work of Burq [13] and Cardoso- Vodcv [22], [24]. Thus the main 
novelty here is the ability to treat medium energies A ~ 1 and low energies A <C 1, 
as well as the universality of the technique developed to treat all ranges of energies. 



The above theorem only reprises one component of Proposition 1.2, namely the 

estimates (1), (2). Using those estimates, however, one can fairly easily obtain 

convergence properties of the resolvent as £ 0. 

Proposition 1.12 (Limiting values of resolvent). Let the notation and assumptions 
be as in Theorem 1.7. LetO<a < min(l,(7o) and f € H°''^/'^+''{M). Let\±ie = 
with z = a ± b such that b > 0, Then the functions u±s := i?(A ± ie)f obey the 
Sommerfeld radiation condition 

\\{dr T *^)w±ellffO,-i/2+<T(M\/fo) + II ^"""ie II (M\Ko) 

^^^^ < C(M,Ko,A)(A-^(^'^) +e<^(^^^^)^)||/|!^o.v.+.(M) 

where \Vuiu\'^ = r^(|Vu|g — |5rUp). Note that the constants are uniform in the 
choice of e. Furthermore, for a fixed sign ± and fixed A > 0, the functions 
u±e converge in _ff°'~^/^~'^(M) to a limit u±, which is the unique solution in 
if2 -i/2-T(-^j) the Helmholtz equation {H - \)u = f such that (5^ T iX^/'^)u± 
lies m H°-'^/'^+'^'{M\Ko) for at least one <a' <a. 



We prove Theorem 1.7 and Proposition 1.12 in Section 3. 



1.13. Refinements at high energy. At high energies A ^ 1, Theorem 1.7 suf- 
fers an exponential loss in the constants. We now recall the standard quasimode 
example that shows that this loss is necessary for certain manifolds M: 

Proposition 1.14 (Quasimode construction). Let Co > and a > be arbitrary. 
There exists a smooth Riemannian manifold M , which is equal to Euclidean space 
R" outside of a compact set, and a sequence A; — )■ -l-oo andei — >■ 0+ of real numbers, 
as well as functions fi € if°'^/^+'^(M) such that if H := — Am (i.e. V = 0) then 

\\R{Xl ± «£0/dlHO.-i/2-<T(M) > e'^°^ll/;||H0.i/2+<T(M) 

for all I. 



We give this standard counterexample (based on the trapped geodesies in a sphere) 
in Section 13. It shows that we cannot hope to eliminate the exponential fac- 
tor e'^i'^'^)^ from Theorem 1.7, at least for manifolds which contain trapped 
geodesies. 

However, as is well known, the limiting absorption properties improve substantially 

at high energies if one assumes a non-trapping condition on the manifold M. We 
give a (standard) result in this direction, with somewhat stronger assumptions on 
the manifold and potential than is strictly necessary: 
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Theorem 1.15 (Quantitative high-energy non-trapping limiting absorption prin- 
ciple). Let the notation and assumptions be as in Theorem 1.7. Suppose also that 
V obeys the bouncP 

\V{x)\ < 

for some A> and all x G M, and that the metric perturbation error function Cab 
obeys the symbol estimates 

(13) \{rdrrV^yeab{r,y)\<Co.0 

for all multi-indices a, /3 > 0. Suppose also that the manifold M is non-trapping 
in the sense that every unit speed geodesic t i— > ^{t) in M escapes to the infinite 
boundary dM as t ^ ±oo. Then, for A sufficiently large depending on M and A, 

we have 

(14) ||i?(A±i£)/||H..-i/2-.(M) < C(M,A)A('^-i)/2||/||^o,v2+.(M) 
for all0<s<2. 

We prove this Theorem in Section 7 using the positive commutator method and 
the scattering pseudodifferential calculus. 

Remark 1.16. Such results for non-trapping metrics have a long history starting 
with the work [100], [66] for compact perturbations of the free Hamiltonian Hq = 
—Ann. In the general case a related result is contained in [102], see also [13] for 
a qualitative version, and [82]. For the closely related task of establishing local 
smoothing estimates, results analogous to those above were established in [30]. 

Remark 1.17. In between the two extremes of non-trapping manifolds (for which one 
has the optimal A~^/^ decay in the limiting absorption principle at high energies) 
and the highly trapping counterexamples (for which one has exponential growth in 
y/X) there is a complicated range of intermediate behaviour, for instance if there is 
just one hyperbolic trapped orbit. We will not attempt to survey the substantial 
literature on this topic, but see [55], [56], [76], [16], [38], [27], [28], [29] for some 
work in this area. 

Remark 1.18. If we replace the assumption on the potential by the weaker bound 
|V(a;)| < AX^/'^{x)~^~'^'^° then one would have to modify the non-trapping condition 
as the potential term V has the same order as the kinetic term —Am and one would 
have to consider the Hamilton flow associated to the symbol of —Am + V rather 
than the geodesic flow. We omit the (standard) formalization of this fact. 

Note that by combining this Theorem with Theorem 1.7 one can weaken the hy- 
pothesis that A is sufficiently large, instead requiring that A > Aq for some fixed 
Ao > (but then the constant C{M,A) in (14) may depend on Aq). 

Remark 1.19. The proof of Theorem 1.15 is microlocal. The result therefore can 
be strengthened to include the situation where there is trapping. In this case we 
consider the conical set 

A^nontrap = {{x,0 ■ li'^) 7(0) = X, ■j{0) = ^ escapcs to dM, t ±00} 

^Alternatively we may assume that |V(a;)| -|- (a;)|Vy(a;)|g < A(x)~^''o , which allows V to be a 
long range potential. 
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invariant under the geodesic flow in the phase space and a corresponding \l/DO 
-Pnontrap with support in A^nontrap- Then we can replace (14) with 

||i?(A ± i£)Pnontrap/||H'^.-i/2-a(M) ^ ^'(M, A)A*~"^/^ ||/|| jjo,i/2+<7(m) 

1.20. Refinements at low energy. At the low energy limit A 0, the constant 
in (11) (or (10)) blows up like a large negative power of A. For the purposes of 
effective bounds, which are uniform for all potentials V obeying the bound (9), 
such blowup is necessary, as the following result shows: 

Proposition 1.21 (Bcssel equation matching construction). Let Co > and 
cr, (To > be arbitrary. Then there exists a sequence Vm ■ R" — > R o/ spherically 
symmetric potentials obeying the uniform bounds 

sup ||(x)i+'^°V^y„|Uoo(Rn) < oo for all a>0, 

m 

sequences Xm 0+ and Sm — >■ 0+ of real numbers with < Sm < Xm, as well 
as functions fm € iJ°'^/^+'^(R") such that if ■= —A + Vm then Hm has no 
resonance or eigenfunction at zero, and 

||-R(Ato ± «em)/ni||H0,-i/2-<T(M) > A^*^" ||/to||ho,i/2+<,(m) 

for all m. 

We give this standard example in Section 14. This proposition shows that the large 

negative power in (11) is necessary. One can also show that a negative power is nec- 
essary in (10) by combining this example with the Carleman estimate (Proposition 
2.6); we omit the details. 

Despite this example, one can eliminate (or at least reduce) the low energy blowup 
in certain cases. Firstly, in the flat case M = R" the the work of Jensen-Kato [60] 
and Jensen [59] provides small energy resolvent expansions 

(15) (if-A)-i= ^ A^/2B,+a(A^/2) 

J=-2 

for Hamiltonians H = — Ar^ + y in odd^ dimensions n, under appropriate con- 
ditions on V, dependent in particular on the index £, even in the presence of zero 
eigenvalues and resonances. The expansion (15) shows that for a fixed potential V 
the limiting absorption constant for the resolvent {H — A)^-'^ blows up at worst as 
C{M,V)X~^ as A — >■ 0. Note however that (15) is qualitative in a sense that the 
bounds on the norms Bj, as operators between weighted Sobolev spaces, depend 
on the potential V and not just on the bound A for that potential; thus there is 
no inconsistency between (15) and Proposition 1.21. Note however subtle distinc- 
tion between effective bounds (which hold uniformly for all V obeying (9)) and 
ineffective bounds (which are not uniform in V). 

It is quite likely that these expansions could be extended to asymptotically conic 
manifolds. We will not do so in full generality here, but just give the standard 



In even dimensions resolvent expansions involves additional terms containing log A. 
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qualitative limiting absorption estimate in the case that there are no eigenfunctions 
or resonances at zero: 

Definition 1.22 (Resonances at zero). Let A = for some k > 0. We say that 

H has an outgoing (rcsp. incoming) resonance at X if the equation Hu — Xu 
has a solution u G iJ'''~'^/^"'"'^(R") for some a > with the property that u ^ 
U^>o ilO'-V2+<^(R") and drU - inu e /fO.-V2+a' (j-esp. drU + mu e 
for some a' > 0. Note that when A = 0, there is no distinction between incoming 
and outgoing resonances, and we shall talk simply about resonances at zero. 

Proposition 1.23 (Qualitative low-energy limiting absorption principle). Let the 
notation and assumptions be as in Theorem 1.7. Suppose also that V is strongly 
short-range in the sense that (.x)^+^'^"V <E L°°{M). Assume also that the operator 
H := — Am + V contains no eigenfunction or resonance at zero. Then for A, £ ^ 
sufficiently small depending on M and V, we have 

(16) \\R{X ± i£)/||^.,-v2-.(M) < C{M, V^)|A|-i/2||/||^„.,/,+,(^^. 
Furthermore, if a > 1/2, then we can strengthen this further, to 

(17) ||i?(A ± ie)/||H2.-i/2-.(M) < C{M, 1^)||/||ho.i/2+.(m)- 

The proof of this result is essentially an application of the Fredholm alternative^ 
and the theory for the free case V = 0;we give this proof at the end of this section. 
The estimates here should be compared with (11); the constants do not blow up 
as fast as A — > 0, but on the other hand the bounds depend directly on V and not 
just on A. 

Remark 1.24. As is well known, in asymptotically Euclidean case in high dimensions 
n > 5 no resonances at zero can occur; see [59]. Moreover, zero eigenfunctions and 
resonances are non-generic relative to a family of perturbations = —Am + kV; 
see [60]. 

Remark 1.25. The above proposition is only stated for sufficiently small A. But by 
combining this result with Theorem 1.7 we see that it in fact holds for all bounded 
A. If we also assume non-trapping, then by Theorem 1.15 we can extend (16) and 
(17) to arbitrary energies A > 0. 

Remark 1.26. The subject of qualitative limiting absorption principles has a long 
and rich history. Most of the results have been fociiscd on a high or a fixed non-zero 
energy regime in various geometric settings, such as the asymptotically Euclidean, 
conic, hyperbolic and obstacle cases; see for instance [66], [72], [57], [54], [53] for 
a representative set of results. The low energy results have been developed for 
example in [46], [47], [9]. For the operator — AR,n -|- V resolvent behavior have 
analyzed in particular in [1], [60]. [59]. [44]. Recently there have been a lot of 
interest in establishing sub-c!xponcntial. in fact logarithmic, boimds on the resolvent 
at high energies for geometries with sufficiently "small" or filamentary trapped sets, 
see [18], [76], [103], [39]. 



Strictly speaking, wc do not apply the Fredholm alternative directly, as we will need to 
uniformly invert a family 1 + -R()(A ± ie)V of compact perturbations of the identity, rather than 
a single such perturbation, but instead will be using the proof of the Fredholm alternative in our 
arguments. 
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In some cases we can make the absence of eigenfunctions or resonances at zero 
quite quantitative. For instance, when the potential is mostly positive we have the 
following result. 

Proposition 1.27 (Quantitative low-energy limiting absorption principle for mostly 
positive potentials). Let the notation and assumptions be as in Theorem 1.7. Sup- 
pose also that V obeys the bounds 

\V{x)\ < A{x)-^-^^° 

for some A > 0, and suppose also that the negative part V- := max(— V, 0) of V 
obeys the smallness condition 

[ \V.{x)r/^ < I3{M) 
Jm 

for some sufficiently small /3 > depending only on M . Then H has no eigenfunc- 
tions or resonances at any energy A G R (and in particular at X = 0). Furthermore, 
if X,s ^0 are sufficiently small depending on M and A, then 

(18) \\R{X±i6)f\\H2,-u^-.^M) < C(M,A)|A|-i/2||/||^o,i/.+.(M). 

Finally, for all sufficiently small \X\,s and a < min((To,l) we have the following 

uniform bound 

(19) ||i?(A±ze)/||^..-3/.+.(M) < C(I(M),^)||/||^o.i/.+.(M), 

where the constant C{I{M),A) depends on the manifold {M,g) only through the 
isoperimetric constant I{M) . 

We prove this result in Section 12. 

Remark 1.28. Observe that the Euclidean limiting absorption principle, Proposition 

1.2, now follows from Theorem 1.7, Proposition 1.12, Theorem 1.15 (since the 
Euclidean space R" is non-trapping), and Proposition 1.27 (with V = 0). 

Note that Proposition 1.27 already applies in the free case V = 0. Prom this and 
Predholm theory we can now prove Proposition 1.23: 

Proof of Proposition 1.23. We begin with the proof of (17). We write H = Hq-\-V 
where Hq = —Am is the free Laplacian, and denote the resolvents for Hq by 
i?o(A ± ie). We have the resolvent identity 

R{X ± ie) = (1 + i?o(A ± ie)V)-'^RQ{X ± ie). 

From (19) (and duality and interpolation), the resolvent i?o already obeys (17) for 
any a > 1/2, and so it suffices to show that the operators 1 + i?o(A ± i£)V are 
uniformly invertible on iJ^'~^/^~'^(M) for a sufficiently close to (but larger than) 
1/2. Prom the Rellich embedding theorem and the strongly short-range nature of 
V, we know that V maps -^/^-'^(M) compactly to H°^'^/'^+" [M), and so by 
(17) for the free operator, the operators i?o(A ± ie)V are uniformly compact on 
^2,-i/2-a(^) for A sufficiently small. 
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Suppose that uniform invcrtibility fails, then we can find a sequence An ± iSn going 
to zero, functions u„ of unit magnitude in H'^'~^/'^~'^{M), and /„ going to zero 
strongly in H'^-^/'^-'' {M) such that 

(1 + Ro(\n ± ien)V)Un = fn- 

If wc then write {t„ := i?o(A ± is)Vun = fn — and /„ := Ro(\± ie)Vfn, then 
Un has norm 1 + o(l) in H^'~^^^~'^ (M) and lies in a fixed compact subset of that 
space, /„ goes to zero strongly in H'^~^I'^~'^{M), and 

(1 + i?o(A„ ± ien)V)Un = fn- 

Applying Hq — (A„ ± ie„) to both sides, we conclude that 

Hun = (A„ ± ien)un + {Hq - (A„ ± ie„))/„. 

By compactness, we may pass to a subsequence such that {<„ converges strongly in 
jj2 -i/2-<T(^) to a limit u of norm 1. Taking distributional limits, we conclude 
that Hu = in the sense of distributions. But this contradicts the hypothesis that 

H has no resonance or eigenvalue at zero. 

The proof of (16) is similar, but one takes cr close to zero instead of to 1/2, and 
uses iJ0'3/2-'^(M) instead of H°''^/'^+'^{M) (and relies on (18) and the adjoint of 
(19)) to mediate between Rf){X ± ie) and V; we omit the details. □ 



1.29. Applications. We now give some applications of the above limiting absorp- 
tion estimates. We begin with some (well-known) spectral consequences. 

Proposition 1.30 (Absence of embedded eigenvalues or resonances). Let the no- 
tation and assumptions be as in Theorem 1.7, and let X > 0. Then there are no 
non-zero eigenf unctions or resonances (incoming or outgoing) at A. 

Proposition 1.31 (Absence of embedded singular continuous spectrum). The spec- 
trum of H on (0, +oo) is purely absolutely continuous. 



We prove these results in Section 15.1. 

In Section 15.5 we give a version of the celebrated RAGE theorem (cf. [83], [2], 
[42]). 

Proposition 1.32 (RAGE theorem). Let the notation and assumptions be as in 
Theorem 1.7. Let f G L^{M) be orthogonal to all eigenf unctions of H. Then for 
any compact set K, we have 

(20) lim ||e^*^/||i2(^) = 0. 



Similarly, we have 



lim \\u{t)\\HHK) + \\dMt)\\mK)=0 



for a solution of the wave equation dfu + Hu = with initial data (w(0), dtu{0)) € 
H^{M) X L'^{M), orthogonal to the eigenfunctions of H. 
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Next, we use the limiting absorption principle and the RAGE theorem, together 
with a closely related result that gives Holder continuity bounds on the resolvent, 
to derive the limiting amplitude principle for the wave equation. 

Proposition 1.33 (Limiting amplitude principle). Let the notation and assvmp- 
tions be as in Theorem 1.7, let f € -L^(M) be compactly supported, and let uq G 
H^{M), u\ e L'^{M), and fi > 0. Assume that f,uo,ui are all orthogonal to all 
the eigenf unctions of H. Let u : R x M — )• C be the solution to the inhomogeneous 
wave equation 

(21) dfU- Amu = e'"*/, u\t=o = uo, dtu\t=o = ui. 

Then for any compact set K c M we have 

\\u{t) - e''^*v\\Hi{K) as i +oo 
where v is the outgoing solution of the Helmholtz problem 

{H-i^^)v = f, 

i.e. V IS such that {dr - i\iJL\)v G iJ0'-i/2+'^'(R»). 



We prove this result in Section 15.6. 

Now wc give a global-in-time local smoothing estimate for the Schrodingcr equation. 

Proposition 1.34 (Global-in-time local smoothing estimate for H). Let the no- 
tation and assumptions be as in Theorem 1.7. Suppose that M is a non-trapping 
manifold, that V obeys the bounds 



\V{x)\ < A{x) 



-2-2(70 



and that H := —Am + V has no eigenfunction or resonance at zero. Let u : 
R X R" C be a solution to the forced Schrddinger equation 

iut — Hu = F. 

and let Pr denote the projection on the continuous spectrum of H. Then we have 

the estimate 

(22) 

JjH'^'PHu{t)r„o.-y.-,^M) dt < C{M,V,A) {\\ummn(M)+ jjnmm,y.+^ 
If furthermore a > 1/2, then we have the variant estimate 

(\\H^'^PHu{t)\\m,-u^-^(M) + \\PHu{t)\\U-u^-^m) ^ C{M,V,A)x 

If we assume that V obeys the hypotheses of Proposition 1.27 then we may eliminate 
the projection Pr and replace the constants C{M,V,A) here by C{M,A). 
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Finally, if we assume symbol estimates on both the metric coefficients Cab and the 
potential V then we can replace the last estimate with the family of bounds 

I ||Pffu(f)||^3+i/2,-i/2-„(M) dt<C{s,M,V,A)x 

X (ll«(0)|lH3(M)+^l|P(i)llH-V2,i/2+.(M) dt\j 

for all s > 0. 



Wc prove this proposition in Section 15.3. 

Remark 1.35. The nontrapping assumption of Theorem 1.34 can be removed if one 
replaces the projection Ph by the projection on the continuous spectrum with 
energies A < A < oo. 

Remark 1.36. The close connection between limiting absorption principles and lo- 
cal smoothing (they arc essentially Fourier transforms of each other, with t being 
the dual variable to A) was first observed by Kato [63]. See [87] for some further 
discussion. In Euclidean space, the global-in-time local smoothing estimates were 
first established in [88], [101], [31]. In order to obtain estimates which are global 
in time, it is necessary (by the uncertainty principle) to establish limiting absorp- 
tion principles at very low energies; high-energy analysis alone is only suSicient to 
establish local-in-time local smoothing estimates. 

Remark 1.37. It is very likely that the above global-in-time local smoothing esti- 
mate will imply global-in-time Strichartz estimates, by adapting the arguments in 
[89], [94], [10], [11], [69]. This would allow one to create shghtly more quantitative 
formulations of some of the Strichartz estimates in [11] and [69], however these 
improvements seem to be rather minor and so we will not detail them here. 



The limiting absorption principle for the Hamiltonian H also leads to the integrated 
local energy decay for the wave equation. 

Proposition 1.38 (Integrated local energy decay). Suppose that M is a non- 
trapping manifold, that V obeys the bounds 

\Vix)\ < A{x)-^-^^°, 

and that H := —IS.M+V has no eigenfunction or resonance at zero. Let u : R"+^ — >• 
C be a solution of the wave equation 

(23) d^u + Hu = F. 

and let Ph denote the projection on the continuous spectrum of H. Then we have 
the estimate 



(24) 



/ (||5,PH«(i)||^o,-i/2-.(M) + ||VPH«(i)||^o,-i/2-.(M) dt) < C{M,V,A)x 

X (l|Vw(0)||2,(^,) + ||7Xt(0)||i2(M) + JjF{t)\\lo.^f2^.^M) dt^ ■ 
For a < min(l,ao) the retarded solution 

.,..«):=/ ^M^n*')^*- 

Jt'<t V rl 
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of the intiomogeneous problem (23) obeys additional bounds 
(25) 

- dr)PHUret{t)\\%o.-i/2+^(^M) + 1^ ^Va,PHUret(0llHO.-i/2+<,(M) dt^ 

<CiM,V,A) [ ||F(t)||^o,i/.+.(M) dt. 

J R, 

A similar estimate holds for the advanced solution 

r sin((t - t')VH) 

'Wadv(t) := - / 7^ F{t ) dt 

Jt'>t y ti 

with {dt + dr) in place of {dt — dr). 

Furthermore, even if M does not satisfy a non-trapping condition there exists a 
compact set K G M such that 

(^\dtPHu{t)\\]j,,.,,,.^^^^j,) + ||VPHu(i)||^o,-i/2-.(M\;f) dt) < C{M,V,A) 

X (l|Vu(0)||l,i(M) + II «t (0)11 1,1 (M) + _^II^^WIlHi,im.(M) dt^ . 

If we assume that V obeys the hypotheses of Proposition 1.27 then we may eliminate 
the projection Ph and replace the constants C{M,V,A) by C{M,A). 

We prove this in Section 15.3. 

Remark 1.39. The statement of integrated local energy decay for _ffo = ^ Ar^ 
goes back to Morawetz [74] . The proof of such estimates for solutions of the wave 
equation on black hole spacetimes, with geometries which are quite different from 
the ones considered in this paper, have been instrumental in understanding their 
global behavior. See [7], [33], [69] for Schwarzschild, [34], [96], [3], [36] for slowly 
rotating Kerr and [37] for the general sub-extremal Kerr case. 

Our final results are pointwise decay estimates for the solutions of the Schrodinger 
and wave equations, obtained by commuting these equations with a Morawetz-type 
operator, applyign energy estimates, and using an iteration argument to amplify 
the resulting decay. 

Proposition 1.40 (Decay for Schrodinger). Let M be a non-trapping manifold 
with metric g given in (4) with ctq > 1/2. We assume that h[r\ := h-\- e satisfies 
the following assumptions. 

\{rdrf{VZ)h[r\\<Cka, A;<3, |q|<2. 

Let ^ : R X R" ^ C fee a solution to the Schrodinger equation 

idti) + AmV' = 0. 
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Then for any t > we have the dispersive estimates 

\m)\\L^(M)<Ce{l+t)-i+'\mO)\\H2,HM), V£>0, Tl = 3, 

||V'(i)llLp(M) < C(l + t)-2+f||^(0)||^2,i(M), V2<p<oo, n = 4, 

wmw^^^^^ < c{i + tr^mwH-^HM), n>5. 

Remark 1.41. It is well known that a solution of the free Schrodinger equation 
idfip + Art.^ = satisfies the dispersive estimate 

||V'(i)IU-(R„) <Ci-"/2||Vo|Ui(Rn). 

As a consequence, we do not believe the rates of decay given by Proposition 1.40 to 
be sharp, especially in higher dimensions. Nevertheless, they appear to be novel in 
such general setting, and we give them here to illustrate an application of how the 
limiting absorption principle (or more precisely, the global-in-time local smoothing 

estimate) can bo used to obtain dispersive inequalities. 

Remark 1.42. Dispersive estimates for solutions of the Schrodinger equation with 
H = — Arti + V, projected on the continuous spectrum of H and assuming absence 
of zero eigenvalues and resonances have been proved in [78], [61], [45], [4]; see also 

the survey [84] and the references therein. 

Proposition 1.43 (Decay for wave). Let M satisfy the assumptions of Proposition 
1.40 and let u : R^^^ C be a solution of the wave equation 

d^u — Amu = 0. 

Then in dimension n = 3 we have 

\Ht)\\Loo^M) < a{l + t)-^+' (||V«(0)||hi,1(M) + ||Wt(0)||Hi.=(M)) , 

for all £ > 0. 

Remark 1.44. A modification of our arguments also gives the variant estimate 

ii^wii ^ < c{i+tr' (iivt.(o)iiHi.i(M) + \\utm\m.HM)) 

L "-3 (M) 

in higher dimensions n > 3. 

Remark 1.45. Decay estimates for solutions of the wave equation in Minkowski 
space, M = R", have, of course, a very long history, including the Huygens prin- 
ciple in odd dimensions and the uniform t-decay with the rate of A 
quantitative approach to decay in Minkowski space has been developed by Klain- 
erman, in what is known as the vector field method, [64]. Qualitative decay results 
in non-trapping geometries have been obtained in the pioneering works [65], [99], 
[75], [70]. For the problem with — Art, +V see e.g. [43]. Quantitative decay results 
for solutions of the wave equation on black hole spacetimes have been obtained 
in [7], [33], [34], [35], [3], [36], [95], [37]. A general approach to the derivation of 
decay from the integrated local energy decay statements have been developed in 
[35] , where the arguments can be in particular adapted to the asymptotically conic 
case, and in [95] for stationary asymptotically Euclidean spacetimes. See also [104] 
for applications to nonlinear problems. The results here follow an earlier approach 
of [33] and serve as an illustration. 

We prove these propositions in Section 16 and Section 17 respectively. 
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2. Key estimates 



Throughout this paper we use the notation and assumptions of Theorem 1.7. All 
constants are henceforth allowed to depend on M, n, <to,A, and a. 



Fix A, £ > 0, and suppose that we have a solution to the resolvent equation 

(26) {H-{X±is))u = f 

(or equivalently that u = R{X ± ie)f) and the closely related Helmholtz equation 

(27) {-Am - z'^)u = F, z'^ = X±ie 

Of course, the former equation can be viewed as a special case of the latter with F := 
/ — Vu. To avoid technicalities, we shall always make the qualitative assumptions 

uJ,FGL\M). 

In this section we lay out the fundamental estimates for these equations which we 
shall repeatedly use in our analysis. The proof of these estimates will be deferred 
to later sections. We begin with a simple charge estimate. 

Lemma 2.1 (Charge estimate). Let X,s,a > 0, let ± be a sign, let f e if°'^/^+'^(M), 
and let u = R{X ±ie)f G L^{M). Then we have 



(28) e |wP (^5 < ll/llH0,i/2+a(M)||M||jf0,-i/2-<,(M). 

Jm 

Moreover, for any a G R we have 

(29) [ {xf'^\u\^dg<C{a) ( {{xn\^ + {x)-^+^''\Vu\l) dg, 

J(x)>2R J(x)>R 



This easy estimate can be established by integration by parts, and is proven in 
Section 4.2. Roughly speaking, this estimate allows us to handle most of the terms 
in our analysis which contain a factor of e, and which do not have derivatives on u. 

Another integration by parts gives the following standard elliptic estimates, proven 
in Section 5: 

Lemma 2.2 (Elliptic regularity). Let H = —Am + V where V is a short-range 

potential obeying (9). Let m G R, and let u, f G H'^'™{M) which satisfies the 
Helmholtz equation {H — z^)u = f in the distributional sense (at least) with z^ = 
A ± ie and A > 0. 
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• We have u G with the elliptic regularity estimate 

(30) \\u\\H2.m.(M) < C("l)(||/|!ffO,™(M) + (1 + A)||?i||HO,™(M)). 

• For any R > 2Ro, we have the localised elliptic regularity estimates 

(31) / \Vu\ldg<c[ (A + i?-2)|«|2 + 

JR<{x)<2R Jr/2<{x)<3R A -\- U 

and 

(32) / \Wu\ldg<c[ {\ + R-^)\u\^+R^\f\^ dg 

Jr<(x)<2R J R/2<{x)<3R 

and 

(33) 

/ \V%\l{xf"' rfff < C(m,i?)[||/||l,o..(M) + (1 + A) / \u\l{xf"^ dg]. 

J{x)>R J{x)>R/2 

• We have the energy estimate 

(34) A^||m||^o,-1/2-^ < C (||Vm||^0,-1/2-^ + ||u||jJ0,-3/2-^ + 11/11^0,1/2 + ^) . 

• If e < CX for some constant C > then we have the charge-type estimate 

(35) 

£||V"IU2(M)||w||l,2(M) < C||/||jJ0,l/2+<, (||Vu||jJ0,-l/2-a + ||u||h0,-3/2-<, + || / 1| JJO, 1/2+,, ) . 

as well as the localised estimates 

(36) 

[ \Vu\l [ \u\l<C [ {x)'+"^\f\' 

J{x)>2R J{x)>2R J{x)>R 

X / {{x}-'-'^\Vu\l + {x}-'-'^\u\' + . 

J{x)>R 

and 

(37) 

A / {x)-'-'^\u\' <C I {{x)-^-^^\Vu\l + {x)-'-'^\uf + {xy+'^lff) . 

J{x)>2R J{x)>R 

Very roughly speaking, these estimates ahow us to exchange derivatives on u with 
factors of A^/^ or 1 + A^/^ whenever necessary. 

Next, by applying the positive commutator method^ to a first-order differential 
operator (a variant of dr) we obtain a useful estimate which allows us to control 
the portion of u in the far region r 1 by the portion of u in the intermediate 
region r ^ 1. 

Lemma 2.3 (Pohozacv-Morawetz type estimate). Let the notation and assump- 
tions be as in Theorem 1.7. In addition we require that e < CX for some positive 



^One can also view this method as another application of integration by parts. 
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constant C. If ro > Ro is a sufficiently large number (depending only on M and 
A) then we have 

[ (a;)-3-2-|M|2 dg < 

J(x)>2ro 



C{ro)( f (l + A)H'rfff + ||/||^o,v2+.(M) 

\Jro<{x)<2ro j 

and 

[ {{x)-'-''^\Vu\l dg + X{x)-'-'^\u\') dg 

J{x)>2ro 



< C{ro) if (1 + A)|«p dg + ||/||^o.v.+.(m) J • 

\Jro<(x><2ro / 

We prove this result in section 6.3. The positive commutator method, when applied 
to a suitable pseudodifferential operator, will also give Theorem 1.15; see Section 
7. 

We also need the following variant of the above estimate, which we prove in Section 
11.3, using a more refined analysis based on spherical energies rather than on the 

positive commutator method. 

Lemma 2.4. Let u € H^{M) be a solution to the resolvent equation 

(iJ - (A ± is))u = f, X±is = z"^ = {a± ihf 

for some A, £ > with the Hamiltonian H = — Am + V satisfying the assumptions 
of Theorem 1.7. Then for any < a < min(l,o-o) and a sufficiently large ro, 

(38) 

\\u\\H0.-3/2+a(^M2ra) - C(^o)(||/||ho,i/2+<.(m^o) + (1 + '^^^^)II"IU=(M^o\M4^(,)), 

l^l^^^ll"llHO.-l/2-<T(M2ro) + ll^"llH<'.-l/2-<.(M2ro) - '^('"o) II / II 1/2 + <t (M^g ) 

(39) + C(ro)(l + AV2)||„||^2(m,^\m..„) 

with Mra = {x & M : {x) >ro}. Furthermore, 

|k"VwU||HO,-i/2+a(M2^o) 11"^ T«2:M||HO.-i/2+<.(M2ro) < C'(^o)||/||ff0.i/2+.(M,„) 

+ C(ro)(l + Ai/2)||^|| 

Remark 2.5. The result of Lemma 2.4 should be compared with Lemma 2.3. First, 
the analog of (38) in Lemma 2.3 provides control of the A^/^||u||jjo,-3/2-<t(M2^^) 
norm. This improvement will be important in our analysis of the low energy regime. 
Furthermore, in Lemma 2.3 validity of the estimate (39) is restricted to the region 
£ < ex. Removal of this condition is crucial for establishing a uniform (in the 
domain £,A) limiting absorption principle. Finally, (40) is the key ingredient in 
establishing the Sommerfeld radiation condition of Proposition 1.12. 



We will need to complement Lemma 2.3 by the following estimate, which controls 
the portion of u in the near region r <C 1 by the intermediate region r ^ 1. 
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Proposition 2.6 (Unique continuation estimate). Let K be a compact subset of 
M, and let K' be a compact set contained in the interior of K . Suppose that the 
potential V obeys the bound \V{x)\ < A for x £ K. Then if u £ H^{M) solves the 
Helmholtz equation (i? — A =F i£)u = f on K, then we have the estimate 

[ \u\'+\'^u\l dg < C{K, K', A)e^(^'^' ( j \f\^ dg + j \u\^ + | V«|^ dg 
Jk \Jk Jk\k' j 



One immediate consequence of this proposition is that any sohition of the equation 
[H — X)u = which vanishes near the boundary of K, also vanishes on the interior 
of K, which explains the terminology "unique continuation". Proposition 2.6 is 
essentially due to Burq [12], but in the interest of self-containedness we provide 
a proof in Section 8 using the Carleman inequality method (which, again, can be 
viewed as a type of integration by parts). Notice the exponential dependence on 
the energy here, which will eventually lead to the exponential factors in Theorem 
1.7. 



To conclude the proof of limiting absorption for general manifolds, we will need to 
analyse how solution u decays at infinity. To this end we introduce the spherical 
masses 

(40) A^[r]:=r"-^/ \u{r,uj)f dh[r]{oj) 

JdM 

for r > i?Oj where the angular metric h[r] was defined in (5). 

Remark 2.7. In the Euclidean case, at least, we expect u to decay like 0{r~^"~^')^'^), 
and so we expect A4[r] to stay bounded. 

Lemma 2.8 (Dichotomy). Let Ci ^ Ro be a large number, and then let C2 ^ Ci 
be an even larger number. Let Mci/2 := {x & M : r > Ci/2}. Let u e H^{Mci/2) 
be a solution to the resolvent equation (26) for some A, e > and sign ±, so that u 
is also a solution to the Helmholtz equation (27) with F := f — Vu ± isu. Suppose 
also that we have the normalization 

(41) ll'"|l_fJ-0,-l/2-<T(Mcj/2) = 1; ll.f IIhO.V2 + <t(Mcj/2) = ^ 

for some 5 > 0. Assume that a < 1/2, that S < 1 and that e < X. If C\ is 
sufficiently large (but not depending on X), and C2 is sufficiently large depending 
on Ci (but not on X ), then one of the following must be true: 

• (Boundedness) There exists a radius Ci < rg < C{Cx,C2) such that 

(42) M[r] < C(Ci, C2)(A-^(<^^'^=) + 1),5 for all ro/2 < r < 4ro. 

• (Exponential growth) For all Ci < r < lOCi, we have 

(43) -^M[r] < -C2(l + X''/^)M[r]. 



This lemma is proven by an ODE analysis of the equations of motion obeyed by the 

mass A4 [r] and several other related quantities; it is the main technical innovation of 
this paper and is proven in Section 10, after some important preliminaries in Section 
9. Observe that this is a "black box" result, in that no assumptions whatsoever 



20 



IGOR RODNIANSKI AND TERENCE TAO 



are made eoncerning u in the region r < Ci/2. The dichotomy between (42) and 

(43) may seem strange, but one way to motivate it is by considering the free case 
M R", F 0, / = 0, e with an ansatz u{rLu) = r-^"-'^'>/'^v{r)Yi{uj) for some 
spherical harmonic Y; of order I, in which case (26) simphfies to the Bessel equation 

Vrr + (A - ^^ )V = 

where L := I + ^^y^; the quantity Ai[r] is essentially just |w(r)p in this case. The 
solutions to this ODE can be described in terms of Bessel and Hankel functions. 
By analysing the asymptotics of such functions one can observe that all solutions 
are either bounded or grow exponentially as r decreases from infinity'^, which helps 
explain the above dichotomy. 

3. Proof of main theorem 

We now show how the estimates in the previous section can quickly be used to 
deduce the main limiting absorption principles in our paper, namely Theorem 1.7 
and Proposition 1.12. (The remaining limiting absorption principles will require 
slight modifications of the above estimates, which we shall address in later sections.) 

The strategy will be as follows. Proposition 10.1 shows us that the spherical mass 
A4[r] exhibits either boundedness (42) or exponential growth (43). The exponential 
growth scenario (43) will turn out to be incompatible with the unique continuation 
estimate in Proposition 2.6, and so we in fact have boundedness for u in the inter- 
mediate region r ^ 1. We will then use Lemma 2.3, Proposition 2.6, and Lemma 
2.2 to recover control of u in other regions of space (and with higher derivatives), 
thus leading to the full limiting absorption principle. 

3.1. Proof of Theorem 1.7. We begin by proving Theorem 1.7. Let the notation 
and assumptions be as in the theorem. We allow all constants to depend on M, A. 
Write u := R{X ± ie)f, thus 

{H-{X± te))u = f 
and u e iJ2(M) (by (28) and (30)). 

We begin by proving (10). We wish to show that 

(44) lhllH..-i/2-.(MW) < C(A-'^ + X^/')\\f\\HO^^;.+.^M) 

for s = 0, 1, 2 and for sufficiently large K. By interpolation it suffices to establish 
the cases s = 0, 2. Prom (33) we see that the s = 2 case follows from the s = 
case, so it suffices to take s = 0. 

Henceforth s = 0. We may now assume £ < A, since the claim (44) follows from 

(28) otherwise. 



^Of course, exponential growth as r decreases from infinity is the same tiling as exponential 
decay as r — >■ oo. But in our application, it is best to think of r = oo as the "initial" value of r 
for equations such as the Bessel equation. 
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Let ro be a large constant, and take M^^j^ := {x € M : (x) > ro/2}. We then 
normalize 

(45) ll"llHO.-l/2-<T(Mro/2) = II/IIh0.1/2+<'(M) = 

We may assume that 5 < 1, since the claim (44) is immediate otherwise. 

Let ro <C Ci <C Ci be large numbers (depending on M and A, but independent 
of A) to be chosen later. If one then applies Lemma 2.8, we conclude that we are 
either in the boundedness scenario (42) or the exponential growth scenario (43). 

Suppose we arc not in the boundedness scenario (42), so that we are necessarily in 
the exponential growth scenario (43). Then we see in particular from Gronwall's 
inequality that 

/ \u\^ dg > e^^(i+^) / |w|2 dg. 

JCi<{x)<2Ci ^3Ci<{x><10Ci 

Using the energy inequality (32) and the normalization (45) we have 

/ \Vu\^ dg<CiCi)l [ {1 + X)\uf dg + sA . 

J4Ci<{x)<8Ci \J3Ci<{x)<10Ci J 

We may assume that the first-term on the right-hand side dominates the second, 
since we are not in the boundedness scenario (42). We then conclude that 



(46) 



/ d5 > e^^^(i+^) / {\u\^+\Vu\l) dg. 

JCi<{x)<2Ci JiCi<{x)<8Ci 



The constant C2 is sufficiently large compared to Ci. Then using Proposition 2.6 
with K = {x: {x) < 8C1} and K' = {x : (x) < 4Ci} wc obtain 

/ (|«r + IVt.1^) dg < C(COe^(^^)^ / \f\' dg 

J{x}<8Ci J{x}<8Ci 

Applying (46) again we conclude (again taking C2 large compared to Ci) 

\uf dg < 6-^^^(1+^)52. 



/ 

J4C 



l'lCi<{x)<8Ci 

Applying Lemma 2.3 we conclude 

A(x)-i-2^|u|2 dg < C(Ci)(l + A)e-3C^=(i+^)52 

'{x)<8Ci 

and the claim (44) follows from (45). 



/ 

J(x 



Now suppose instead that we are in the boundedness half of the dichotomy, thus 
there exists a radius Ci < ro < C(Ci, C2) such that 



/ 



dg < C(Ci,C2)(A-^(^i'^=) + 1)6. 



lra/2<{x)<iro 

Using Lemma 2.3 as before, wc conclude 

/ A(x)-i-2-|w|2 d5<C(ro,Ci,C2)(l + A)(A-^(^^'^^) + l)5 

J{x)>'lro 
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and (44) follows. 

Finally, we prove (11). In view of (30), it suffices to show that 

||m||jj0,-1/2-<t(m) <C{X *^ + e'^^)||/||^0,l/2+a(M)- 

But this follows from (10) and Proposition 2.6. This completes the proof of Theorem 
1.7. 

3.2. Proof of Proposition 1.12. We now turn to the proof of Proposition 1.12. 

The Sommerfeld radiation condition (12) follows immediately from Theorem 1.7 
and (40). Their combination gives the estimate 

lk"^Va;U||i/o,-i/2+.(M2,^) + ||Wr-i\/Au||^o,-i/2+.(M2,^) < C(A-'^+e'^^)||/||^o.i/2+.(M)- 

and completes the proof of the first part of Proposition 1.12. 

To show that for fixed A > the the functions u±^ = R{X ± ie)f converge in 
^0,-1/2-CT unique solution u± of the Helmholtz equation {H — X)u = f selected 
by the requirement that u e iJ2,-i/2-<T(^) {drTiX^^^)u G H°'-'^/'^+''' {M\Ko) 
for some cr' > we consider the difference 

with = ^ + i^i,2 and Iin(2;i,2) ~ £1,2^"^^^ > (we assume that A ^ £1^2) and 
ei > £2 (convergence for u-^ can be treated in a similar fashion.) The function s{x) 
is assumed to be smooth and obeys the requirement that s{x) = r ior x & M\Kq. 
The function w verifies the equation 

{H -\-iei)w = G 

where G is given by the formula 

G := i(^2-^i)e^("^-"=)^(") (2(VfesV'=-i^2 |Vs|2)+AMS+i(^i+^2)(|Vs|2-l))«+,^. 

For values of £2 < £1 <C A we have that 

\zi -Z2\< 4Im(zi - Z2) ~ \-^'^{ei - £2). 

Therefore the first part of Proposition 1.12 and the choice of function s{x), which 
in particular gives that |AMs(a;)| < C{x)~^, imply that for any < a' < a and 
any a > we have the estimates 

||gi(.i-.2).(x)(y^gyfc - iZ2 |Vs|2)«+,J|^o,i/2W(M) 

< C\z^ - ^21"^+'^-'^' (A-^ + e^^)||/||H0,i/2+.(M), 

||gi(..-.2Mx)(|^^g| ^ Iz^+Z2\ |VS|2 - l)w+,J|^o,i/2+.'(M) 

<C||e'(^-^=)^(-)(a.)-i«+,J|^o,v2+.'(M) 

< C\zi - Z2r'-"{X-'' + e^^)||/||H0,i/2+.(M)- 
Thus for any 0<a'<a<l — awe have 

\\G\\jjo,i/2+,^'(^M) ^ (^1 ~ £2)'^"'^ C(A~'^ + e'^^)||/||jj-o,i/2+<,(M), 
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By Theorem 1.7 and the first part of Proposition 1.12, the solution of the problem 

{H — X — ie2)w = G obeys 

(47) 

\\{dr-\^'^)w\\HO.-u^+.'(M\K,) + \MH-.-u^-^m ^ (ei-£2)"-"'C(A-^'+e^^)||/||^o,i/2+. 
Since by Theorem 1.7, 

< C(ei - £2)'^"'^ ll"+e2llH"'-i/2-».'(M) 

< (£1 - £2)^-'^'C(A-^ + e^^)||/|Uo.i/2+.(M). 

we obtain 

which implies that is a Cauchy sequence in ff^'~^/^"'^(M) converging to 
the limit u+ e i?2'~^/^"'^(M). Moreover, (47) also implies {dr - iX^/^)u+ e 
^o,-i/2+<r'(^ \ i^o) for any a' < a. 

To prove uniqueness of solutions of the Hclmholtz equation 

{H-iX± te))u = f 

with the property that u e H'^''^/^-'' (M) and {dr T iX'^^^)u e iJ0'-V2+'^'(M \ Kq) 
for some ct' > we simply repeat the argument above with the functions 

v{x) = e^'^^^-^>^''^u{x) - u±^{x) 
and conclude that u±^ converges to u. 

This concludes the proof of Theorem 1.7 and Proposition 1.12, except for the proofs 
of the estimates in Section 2. These estimates will be the focus of the next few 
sections. 

4. Conservation laws 

All of the results stated in Section 2 rely, fundamentally, on integration by parts 
arguments. It will be convenient to present these arguments in a unified framework, 
namely that of exploiting the conservation laws for the Helmholtz equation (27) 
using the abc method of Priedrichs. While this framework is in practice too general 
to use directly, it does provide a convenient way to display the interrelationships 
between the more specialized methods we study below. Thus we shall devote this 
section to the general abc method for the Helmholtz equation (27). As a quick 
application, we will be able to establish the charge estimate. Lemma 2.1. 



4.1. Densities and currents. Let u € H'^{M) be a (complex-valued) solution to 
the Helmholtz equation (27) for some A > 0, so that F G L'^{M). We introduce a 
number of (real) tensor fields on M, namely 
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• the charge density q := [itp, 

• the energy density e := |Vu|g := Re(V'^uV/cu), 

• the charge current j'^ := Im(u V'^u), 

• the charge gradient v'^ := Re(u V'^u), 

• the energy-momentum tensor Q™*^ := Re(V'"wV'^«) — ^g"^''{e — \z\^q). 

Note that from the assumption u G H^{M) that all of these tensor fields are ab- 
solutely integrable. A direct computation (recalling that the Levi-Civita covariant 
derivative V respects the metric g and is torsion free) yields the charge gradient 

identity 

(48) Vfeq = 2vfe, 
the charge conservation law 

(49) Vfej'' + Im(z2)q = -Im(l<F) 
the identity 

(50) Vfev'^ = i Awq = e - Re(z2)q - Re(uF) 
and the conservation of energy-momentum 

(51) VftQ"'= + Re ((z^ - \zf)W"'uu) = -Re(V^i^). 

The energy density e does not directly obey a useful equation of the above type (the 
expression Vaq contains expressions involving other double derivatives of u than the 
Laplacian A^u). However, this energy density is clearly related to the other tensor 
fields, for instance we have the Cauchy-Schwarz inequalities \j\g, |v|g < q^/^e^/^ and 
also that Q = 0(e) + 0(Aq). 

The conservation laws (49), (50) and (51) are directly related to the conservation 
law of the energy-momentum tensor 

Q"'^ := Re(V^V^u) - ^g"''(V^ V^^u) 

associated with the wave equation DgU = with the metric ga^dx^dx^ = —dt^ + 
gijdx'^dx^ . Conservation of energy- momentum is expressed in the form 

S/mQ''"' = 0. 

The relationship between conservation laws for the Helmholtz and wave equations 
becomes apparent under the formal rule 

d . 



4.2. Proof of Lemma 2.1. As a quick application of the charge conservation law 
(49) to the resolvent equation (26), we now prove Lemma 2.1. 

Let the notation and assumptions be as in that lemma. Then u G L'^{M). From the 
resolvent equation (26) (and the boundedness of V) we deduce that Amw G L'^ (M) 
also. From elliptic regularity we conclude the qualitative fact that u € H'^{M). 
In particular from Cauchy-Schwarz we see that the charge current j" is absolutely 
integrable. Also if we set F := f — Vu then F e L^{M) and so the charge source 
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term Im(uF) is also absolutely integrable. We may thus integrate (49) (by inserting 
a large cutoff function xr, adapted to the ball of radius R, and integrating by parts) 
to obtain the identity 

/ XR\u\'^dg+ / XKlm(t7F)rf5= / VxijIni(VMtl). 
Jm Jm Jm 

Expanding out the definition of F and taking into account that V is real, we thus 
obtain 

(52) ±ie XR\u\'^dg = - XR^^uf) dg + VxflIm(V«tZ). 

Jm Jm Jm 

The first claim now follows after letting R ^ oo, applying the Cauchy-Schwarz 

inequality and using the assumption on V. 

Similarly, integrating by parts with the function 7]^, where r]i{ is a cut-off function 
adapted to the shell of radius R, we obtain 

±ie r]ji\ufdg=- ??^Im(u/) dc/ + 2 / 77flVr?Klm(Vuu). 
Jm Jm Jm 

Therefore, 

[ ^l\u\^dg<c( {\f\^ + R-^\Wu\l) dg. 

Jm Jr<{x)<4r 

The estimate (29) now follows after multiplying the above inequality by R^" and 
summing over all dyadic R> Rq. This concludes the proof of Lemma 2.1. 



4.3. The abc method. Now we consider how to exploit the other conservation 
laws (48), (50), (51). The abc method of Priedrichs is a general way to exploit these 
identities. It proceeds by introducing arbitrary tensor fields a'', b, , and applying 
Stokes' theorem to the (momentum) vector field 

(53) = a^q + bw^ + c^Q*^™ 
to obtain (for a compactly supported cutoff function x) 

(54) / (divP)x dg = - [ P^dkX dg, 
Jm Jm 

and in particular (if P has sufficient decay) 

(55) / dbfPdg = Q. 

Jm 

From (48), (50), (51) we observe that the divergence divP = VkP'' can be computed 
explicitly as 

(56) 

divP = Vfca'^q + {2a'' + \lH)wk + {b- ^Vfcc'=)(e - j^l^q) + (VfcC„)Re(V^V"«) 

+ 6(1^1^ - A)q + CfcRc [{z^ - |zp)V^M^ - 6Re(uF) - CfeRe(V''uF). 

One then hopes to select a'^, 6, c'^, x so that many of the terms in (54) or (55) 
carries a useful sign, so that one obtains a non-trivial estimate on u. 
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In practice, the full generality of abc method is difficult to use, as there are too 
many possibilities for a*^, 5, c*^, x which are available, and too many ways in which 
one could hope to exploit positivity. Instead, one typically specializes to a sub-case 
of the Friedrichs method which has fewer parameters. For instance, if one sets 
a" = = and 6 = 1 one obtains the Lagrangean identity^ 

(57) / ex = / Aqx - / v'^VfeX + / MuF)x- 
Jm Jm Jm Jm 

A further integration by parts using (48) yields 

(58) / ex = / Mx+ \ I qAmX + / Re(wF)x 
J M Jm ^ Jm J m 

(or alternatively one can apply (55) with a*^ = — ^V'^X; ^ = X> c" = 0). We exploit 
this identity in Section 5. 



4.4. Pohozaev-Morawetz type identities. Another way to reduce the parame- 
ters is to introduce the ansatz 

c'=:=V'=VF; 6 = ^Am^F 

for some scalar real-valued function W. Then the formula (56) simplifies (because 
the (e — \z\^q) term disappears) to yield 

divP = Vfctt'^q + {2a'' + ^V''AMW)vk + Hess„fe(W)Re(V^V'=u) 

(59) + ^AWi\z\^ - A)q + VkWRe ((z^ _ \z\'')Vhiu^ 

- i(AMW)Re(uF) - {VkW)Re{V''uF) 

where HesSmfe(W^) := V^j, 14^ is the Hessian of W. If one then sets^ := -jV^'AmH^ 
to cancel the term, and introducing the modified inhomogeneous term 

G:^F^\-\z\^ 

one obtains the Pohozaev-Morawetz identity^^ 

(60) 

/ ResSmk{W){V"'u){V''u) - ];{Al,W)ci + l-{AMW)Re{uG) - {V''W){Re{W^G))x 
Jm 4 z 

= [ (-l(v'=AMW)q+ J(AMW)v'= + (V„W^)Q'"'=)Vfexrfff• 
JM 4 Z 



^The energy-momentum tensor Q'"*' is associated with the Lagrangean C = — 

^To continue our analogy with the wave equation, the above choice of a'' , b, is somewhat 
reminiscent of the construction of a modified momentum, associated with a conformal KiUing 
vectorfield, for solutions of the wave equation. 

^'^Particular cases of this identity have been used by Pohozaev to prove non-existence of so- 
lutions of certain nonlinear elliptic equations, [77], and by Morawetz in the study of long time 
behavior of solutions of a wave equation, [74]. 



QUANTITATIVE LIMITING ABSORPTION 27 

This identity is particularly^^ useful when W is geodesically convex (i.e. the Hessian 
Hess(M^) is positive definite) and the quantity A\,jW is non-positive, as this implies 
that the first two terms on the left-hand side of (60) become non-negative. We will 
also re-interpret this identity in terms of the positive commutator method, and give 
some consequences of this identity, in Section 6. 

To account for the the difference between A and \z\'^ we can modify the momentum 
P as follows 

pfe = -l(V*=AMH^)q+ ^(AmT4^)v'= + (V™W^)Q™'= -Flm^ v'= - Re^j^ 

Then 

divP = Hcss„fe(W^)Re(V^V'=M) - ^Al,Wq + ^AmW{\z\'' - A)q 
+ Imz(e - Aq) + Rezlm{z^)ci + 2lmz V^M^Re (j,z V'uuj 

- ^{AMW)Re{uF) - (VfeVF)Re(V'=uF) - Im(^uP) 
After some calculations we obtain 

divP = Hess„fc(W)Re(V^V'=w) - ^A|fWq ^AmVF(|^|^ - A)q 

+ Imz (e - \V''WWku\^ + \VkWV''u + izuf) 

- ^{AW)Re{uF) - {\/kW)Re{V''uF) - lm{zuF), 

which after integration over M becomes 

(61) 

/ f Hess„fc(H^)Re(V^V'=w) - ^A^^PFq + IamW{\z\^ - A)q + Imz (e - jV'^M^VfeWp) 
+Imz\VkWV''u + izuf - ^{AW)Re{uF) - (VfeVF)Re(V'=uF) - Im(^wP)^ x 

= -j^ (- J(V'=AMW^)q + ^{AmW)v'' + {VmW)Q""' + Imzv" - Rez}''^ VkX 

This identity is useful if the Hessian of W is positive definite, A|f < 0, AmW > 
and IWW] < 1. For instance in Euclidean space R" with the standard metric, the 
choice ofW = \x\ — (1 + la^l)^""^ satisfies all the above requirements and effectively 
leads to the proof of the limiting absorption principle. We will use a version of (61) 
in Section 6 and later in the derivation of the Sommerfeld radiation condition. 



4.5. Carleman type identities. Returning to (59), the choice a'' := jV'^AmW 
is not the only possible choice for a*^. If we set W := e^^ in the above ansatz, and 
then set 

a>' := l-{V''w)AMe^^ - hv''AMw)e^^, 



^^We should note that this identity is most appUcable when A = Izp, i.e. £ = 0, in which case 
it leads to interesting statements even for nonlinear F = F{u). 
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then after some computation using (59) we obtain"'^^ 

divP = {{AMwf + 4\Vw\IAmw + 4Hess„fe(w;)V"wV'=u; + 4|Vw;|^ - ^Al,w)e'''^<i 
+ {4{Vkw)AMW + 4Hess„fe(u;)V"w + 8\Vw\lVkw)e'^'^v'' 

- ^(Ae2"')Re(uG) - (Vfee2"')Re(V^G) 
which after completion of the square becomes 

divP = e^'"|2Vfc«;V*'w + (2|Vw|^ + Amw)u\^ + 2Hess„fc(«;)V™(e"'M)V*^(e'"u) 

+ 2Hess^fcH(V™H(V'He'"'|u|' - \{Aljw)e'^ \u\'' 
- i(Ae2«')Re(uG) - (Vfce""')Re(V^G). 

This identity may seem hke an algebraic miracle, but we can interpret it as another 
positive commutator estimate in the more general setting of Carleman inequalities, 
see Section 8. If « is compactly supported, then one can integrate this identity on 
M (as in (55)) to obtain the Carleman identity^^ 

[ |2VfewV'=u + (2|Vw|^ + AMw)w|^e2"' dg 

J M 

+2 / Hess„fe(u;)V'"(e"'u)V'=(e'"w) dg 

+2 [ Hess„fc(M;)(V'"tz;)(V'=w)|«| V- rf5 = / hAljw)\u\^e''^ dg 

+ / ^(AMe2'")Re(MG) + (Vfce2'")Re(V^G) dg. 
Jm ^ 

If we write U := (2VfewV'=u + (2| Vu;|g + AMw)w)e"' and note that 

^ A 2k; I /V7 2w\-r-7k Wtt 

-Amg u+(\/ke ) V M = e U 
we can rewrite the previous expression as 

llf^lli^(M) +2 / Hess„fe(tz;)V"(e-n)V'=(e-n) dg 
Jm 

+2 / BesSmkiw){d'"w){V''w)\u\^e^'" dg 
Jm 

i(A2,«;)|«|2e2- dg+{e^F,U)L^M) 



I 

Jm 



Note that no curvature terms appear here, because the Levi-Civita connection is torsion free, 
and at no stage do we need to commute the Laplacian Ajvf with a derivative. 

"'^'^It may seem remarkable at first that tiie energy A makes no appearance in this identity. 
But note that if one increments G by Xu then tiie resulting changes in the two terms involving F 
cancel each other out, thanks to (48) and integration by parts. In practice, the requirement that 
u be compactly supported means that we will have to apply a cutoff function to truncate u, and 
this will cause A to enter the identity in a non-trivial manner. 
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and hence by the Cauchy-Schwarz inequahty {e^ F,U) i,-^(^m) < 2l|e"'-P||L2(M) + 
l|f^l|L2(M)) we conclude 

(62) 

/ Hess„fe(w;)V'"(e'"«)V'=(e'"w) dg 

+ f Hess„,H(V™to)(V'=«;)|«|2e2-d5< / ]r{Al,w)\ufe^'"dg+he^F\\l,.^^ 
Jm Jm ^ » ' - 

This identity is useful when w is geodesically convex (so that both terms on the 

left-hand side is positive) and large (so that one can absorb some of the right-hand 
side terms into the left-hand side). We shall exploit it in Section 8. 



5. Elliptic regularity 



In this section we prove the various claims in Lemma 2.2. One of the key tools here 
will be the energy identities (57), (58). 



5.1. Proof of (30). Let the notation and assumptions be as in Lemma 2.2. We 
begin with the proof of the elliptic regularity estimate (30). By absorbing z^u into 

the / term we may assume that A = 0; similarly, by absorbing Vu into the / term 
we may assume that V = 0. From classical elliptic regularity we know that u is 
locally in H'^. Our task is thus to show that 

/ {xy"^{\Vu\l + \V^u\l)dg<Cm [ {xf^{\AMu\' + \u\')dg. 
Jm Jm 

Let us first verify this claim for the first derivatives Vu. Let i? ^ 1 be a large 
number, and let to the region {x) < R which equals one for (x) < R/2. From 
(58) we have 

/ \Vu\l{x)^"'v>R = l [ \u\^Am{{x)^'"^r)+ [ MuAmuKx^^^^r; 
Jm ^ Jm Jm 

applying Cauchy-Schwarz (and the very crude estimate Am{{x)'^™'^Pr) = Om((a;)^™)) 
and then letting i? — )• oo we obtain the estimate for first derivatives. Thus it will 
suffice to show that 



/ 

J^ 



{xf^\^''u\ldg<Cm [ {xf^\AMuf dg + CmWufH^.m^M)- 
Jm 



M Jm 

The standard Bochner integration by parts argument to give this estimate would 
require boundedness of the Riemann curvature tensor. Let us present a slight 
variant of that argument that only requires boundedness on first derivatives of 
the metric. To illustrate the method let us first work locally, using a compactly 
supported bump function We claim the estimate 

II|V^U|X||Z,2(M) < Cm{\\AMu\\HO,'n(^M) + ||w||iJi,m(M))- 
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To see this, wc integrate by parts, computing 

ll|V'u|xlli2(M) = (VfeV^W,x'V^V'=w)L2(M) 

= -(V^w,X^VfcV^V'=M)L2(M) - 2(vVx(Vfex)VfV'=M)i2(M) 
= -(VVx^V^AMu)i2(M) - (VVx^[Vfe, V^JVMl^Im) 
+ C(|||Vxls|Vu|g||i2(M)|||V^u|gxlU2(M)) 

= Om{\\^MU\\%o,m(^M)) +O(||u||Hi."»(M)||AMM||jf0,n.(M)) 

+ 0(|||V'u|gXl|i^(M)||w||Hi,m(M))-(VVx'[Vfe,V^]V'=u)i2(M)- 

Thus to verify the claim, it will suffice to show that 

KV^W,X^[Vfc, V£]Vfeu)L2(M)| < C'm(||w||?fi,™(M) + ll|V^u|x|U2(M)||w||Hi,r"(M))- 

If WC exploited the boundcdness of the Riemann curvature tensor at this point 
(recalling that g was assumed to be smooth), then we would be done (and we 
would not need the second term on the right-hand side). However, one can instead 
work in local co-ordinates, writing^^ 

for any vector field X, where T are the Christoffel symbols. Applying this iden- 
tity and integrating by parts as necessary to prevent any derivatives falling on the 
Christoffel symbols, we obtain the above claim, where wc only needed the boundcd- 
ness of the Christoffel symbols F (which are essentially one derivative of the metric 
g, in contrast to the Riemann curvature which is two derivatives). 

One can perform a similar argument with x replaced by (a;)™(pi{(l — x)) where 
1 — X localizes to the asymptotic region r » i?o; the point being that the conditions 
(6) ensures the appropriate boundcdness of the Christoffel symbols in this region. 
Indeed, there is some additional decay and vanishing properties of these symbols 
arising from the normal form co-ordinate structure, which we shall simply discard; 
similarly wc shall discard any gains of (x) when a derivative hits the (a;)^'" factor. 
We omit the standard details. This concludes the proof of (30). 

5.2. Proof of (33). We now localise (30) to prove (33). Let X be a smooth cutoff 
to the region {{x) > R/2} which equals 1 when {x) > R. Applying (30) to the 

function x^u wc have 

||X^w||ff2,™(M) < C(TO,x)(llxVlU<'.™(M) + ||X^|Vu|||i/o.™(M) + (l + A)||x^w||ffO.™(M)) 

which (by the product rule) implies 

||X^|V2«|||^o,^(M) < C(m,x)(||xVl|HO.™(M) + ||x'|Vw|||HO,™(M) 

By using the interpolating inequality 

||X^|Vu||||^o,m(M) < C'(m,x)||X^|V^u|||H°.'"(M)||Xw||HO,m(M) 

(which is easily proven by integration by parts and Cauchy-Schwarz) we conclude 
that 

||x'|V2w|||ffO.™(M) < C(x)(||/||ffO.™(M) + (1 + A)||xu|ko,™(M)). 

^^Schematically, what we are doing is observing that any integral expression of the form 
VuV^gVu can be rewritten using integration by parts as VuVgV^u plus lower order terms. 
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The claim (33) follows. 

5.3. Proof of (31), (32), (34). If we apply the identity (58) with x equal to a 

smooth cutoff function adapted to the annulus |i? < {x) < 3R which equals one 
when R< {x) < 2R, we obtain the estimate 

/ \Vu\ldg<cf {X + R-^)\u\^ + \u\^\V\ + \u\\f\dg, 

jR<{x)<2R Jr/2<{x)<3R 

and hence by Cauchy-Schwarz (and (9)) we establish (31) and (32). Similarly, 
applying (58) with x = {x)~^~^'^ we obtain (34). 

5.4. Proof of (35), (36), (37). Interpolating (34) ||u|1^o,-3/2-<t we have 

A3||u||ff0.-1-. < C{A) (||V?i||^0.-l/2-. + ||m||^0,-3/2-<, + 11/11^0,1/2+.) . 

Combining (34) with the charge estimate (28) we have that 

'^^^ll^*lli2(M) < C{A)\\f\\lfa.l/2-, (||VM||ffO,-l/2-. + ||u||h0,-3/2-<, + |l/|lff0.1/2+<,) . 

At this point we require the condition e < CX for some constant C > 0. Then 

e^ll"llL2(M) < C'(^) 11/11^0,1/2+,, (||Vw||j^0,-l/2-. + ||u||^0,-3/2-<, + 11/11^0,1/2+.) . 

The identity (58) (with X = 1) implies that 

||VU||Z,2(M) < X^\\u\\l2(^m)+C{A) (^X~i\\f\\H0.1/2+a + A3||u||jj0,-l/2-. + ||u||jyo,-i-.^ , 

Combining the above four inequalities and using that e < CX we conclude (35). 

By inserting smooth cutoffs localised to the region {x : (x) > R} to tlie above 
argument we also obtain (36) and (37); we omit the standard details. 

The proof of Lemma 2.2 is now complete. 

6. The positive commutator method 

We continue our study of the Helmholtz equation (27), and recall the well-known 
positive commutator method to analyze this equation. We let S be an arbitrary 
pseudo-differential operator, and consider the expression («[— Am, S]u, u)l2(m) , where 
i[S, B] := i{SB — BS) is the usual Lie commutator and {u, v) := Re uv dg is the 
real inner product. Then from the self-adjointness of —Am — A and (27) we have 
the commutator identity^^ 

{i[-AM,S]u,u)L2(^M) = {-il-^M - X,S]u,u)l2(^m) 

= {iSu, (-Am - X)u)l2(m) - {i^-^M - X)u, S*u)l2(m) 
= {iSu, F ± ieu)i,2(M) - {F ± ieu, {iS)*u)L2(M) 



^^In applications, u will be in H^(M), F will be in L^(M), and S will be first order, and so 
there is no difficulty justifying the manipulations below. 
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and in particular by Cauchy-Schwarz 
(63) 

{i[-AM,S]u,u)L2(^M) < {\\Su\\H0,-i/2-a(^M) + \\S*u\\fjo,-i/2-,,(^M)) \\F\\HO.i/^+<'{M) 

+ 2e||S'u||i2(M)||u||L2(M)- 

This identity is useful when S is chosen so that i[— Am, S] is a positive operator (at 
least to top order), which is why the use of this identity is known as the positive 
commutator method. 

There are a number of ways to generate an operator S whose commutator i[— Am, S] 
is positive. One is to make A itself equal to a commutator S := i[— Am, W], where 

is a real-valued scalar weight function, interpreted as a pointwise multiplication 
operator u i— )■ Wu. Then S is the self-adjoint first-order differential operator 

Su := -2i{V''W)VkU - i{AMW)u 

and i[— Am, S] can be computed to be the second-order operator 

t[-AM,S]u = -4Vfc(Hess'='"(VF)V„w) - {Al,W)u 

The positive commutator estimate (63) then becomes 

[ 4Hess''"(VK)Re(VlMV™u) - {Al^W)\u\^ dg 
Jm 

('^'^) < 2\\2{V''W)VkU + {AmW)u\\ho,-i/2-,^m)\\F\\hoa/2+.^m) 

+ 2e\\2i\/''W)VkU+iAMW)u\\L2(M)\\u\\mM), 

which should be compared with the x = 1 case of (60). Indeed (60) can be 
interpreted as the positive commutator identity applied to the operator S := 
Xi[— Am, W]; we omit the details of this calculation. 

6.1. A cheap proof of limiting absorption in the free case. Using this in- 
equality we can easily prove the limiting absorption principle^ (or at least the main 
estimates (1), (2)) for the operator Hq = —A in Euclidean space. We choose the 
function 

W(x) = |x|-(l + |x|)l-2-. 

Direct calculation shows that for dimensions n > 3 

\x\\AW\+\'7W\<C, -A''W>c{l+\x\)-^-^^, liesskm{W)>{l+\x\)-^-^^Skm- 
We then obtain from (64) 

/ ((a;)-l-2-|Vw|2 + {x)-^-^-\u\^) < (j|V«||H0,-i/2-. + \\\x\-\\\H0.-r/2-.) 11/11^0.1/2+. 

+ 2e {\\Vu\\l2 + \\\x\-^u\\l2)\\u\\l2, 



This can either be seen by expanding everything out using the Leibnitz rule, and then 
exploiting the symmetry properties of the Riemann curvature tensor, or alternatively noting that 
both sides of this identity are second-order self-adjoint operators with real coefficients and agree 
both at top order and at the constant term, and thus must be identical. See e.g. [51]. 

''^^This "cheap" proof of the limiting absorption principle applies to the region where e < C\ 
for some positive constant C. Analysis of the remaining region requires the use of an additional 
conservation law (61). This issue will be addressed in Section 11 in a more general context. 



QUANTITATIVE LIMITING ABSORPTION 



33 



Now combining Hardy type inequalities 

Ilkr^MlliS < C||Vm||l2, |||a:|"^M||j:i-0,-l/2-<. < C (|| Vu||^0,-l/2-^ + ||m||^0,-3/2-<,) 

with (35) and (34) we obtain 

||Vu||jJ0,-l/2-<, + A2 ||m||^0,-1/2-„ + ||u||^0,-3/2-<, < C 1 1 / 1 1 ^0,1/2+.. 

with second derivative estimates following from the elliptic regularity estimate (30). 

Remark 6.2. A proof of the Sommerfeld radiation condition and the estimate (3) 
will require the use of the identity (61) in place of (60). We will return to this point 
in a more general context when we establish Lemma 11.1. 

6.3. Proof of Lemma 2.3. Now we can prove the Pohozaev-Morawetz type esti- 
mate in Lemma 2.3. Let the notation and assumptions be as in that Lemma. 

For this argument, the function W is chosen to be 

W = x{{x)x){\x\-\x\'-^'^), 

where %('') is a smooth cut-off function vanishing for r < Tq and equal to 1 for 
r > 2ro for some Tq > Rq to be chosen later. On the support of W the metric g 
has the form 

g = dr^+ r^{Kb{u}) + r-^'"'eab{r,Lj))duj''duj^ 

The second fundamental form Qab of the hypersurfaces r =const for a metric in 
this form is equal to the expression 

Qab = rhab + ^dr [r^-^''° eab{r,u))] = rhab + 0(ri-2<^°) 
and its mean curvature 6 (i.e. the trace of Qab) has the asymptotic 

e = r-^{h{uj) + r-^'''e{r,uj)T''Qab = + 0{r-^-^^"). 

The Laplacian Am can be written in this coordinate system as 
(65) AM = d^ + Qdr + AaM^, 

where Aqm^ is the Laplace-Beltrami operator of the hypersurface r =const. In 
particular, Am applied to the function W{r) can be computed as 

AmH^ = Wrr + QWr = Wrr + ^ Wr + 0(r-'^-^''°)Wr. 

r 

Furthermore, 

AlW = {df + VLZld^fW + Am {0{r-^-^''°)Wr) . 
Finally the Hessian of W, 

B.eSSrr{W) = Wrr, HeSSar(W) = 0, HeSSa6(T^) = QabWr. 

From this we easily conclude that 
r\AMW\ + \\lW\g < C, 

-Al > Cxr-^-^'' + Am {0{r-^-^-'^)Wr) + C0(1), 
Hess„fc(W^)>Cr-i-2'^ + CO(l). 
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Here C is a smooth cut-off' function with, support in the region ro < (x) < 2ro 
introduced to account for the derivatives falling on x in the expression for W. Let 
Mr denote the set {x e M : {x) >r}. Then 



cf {{x}-'-^^\Wu\l + {x)-^-^''\uf)< [ Am {0{r-'-^'">)Wr) 

+ / (l^^ls + + (||Vu||ho,-i/2-.(m.„) + lkllff<',-3/2-<,(M )) ||/||ffO,i/: 

JMr-g\M2ro ^ ' 

Arguing as in the proof of the limiting absorption principle in Euclidean space 
above we can, with the help of (36) and (37), reduce the inequality further, 

c / ((a;)-i-2-|V«|^ + {x)-^-'"'\\u\^ + {x)-^-'"'\u\^) < [ Am {0{r-'-^-'')Wr) \uf 

JM2,-g JM 

+ [ (|V«|^ + H^) + ||/||^o./2+.(M.„,.)- 

JM,^/2\M2ro °' 

It remains to deal with the term Jj^ Am {0{r^^^'^'^")Wr) Wc avoid applying 
the Laplacian Am to the term 0{r~^~'^'^°)Wr, as this would require higher differ- 
entiability of the metric g than required in Theorem 1.7. Instead we integrate by 
parts to obtain 

/ Am (0(r- 1-2-0 )W^,) \u\^ = -2 [ (0(r-i-2-<')W,) (Re(AMWt7) + \Vu\l) 

JM JM 

= -2 [ (0(r-i-2-o)W,) (A|«|2 + Re{Fu) + |Vm|^) 

JM 



It then follows that 

'I 

JM2^a 



((x)-i-2-|Vu|^ + (a;)-i-2<^A|u|2 + {x}-^-^" \u\'') 

< I (|V«|^ + H^)+||/|lkv2+.(M.„,,)- 

Jm,^,2\M2.„ 

provided that ct < co and ro is sufficiently large. The desired conclusion now is a 
consequence of (32). This concludes the proof of Lemma 2.3. 

Remark 6.4. Prom the Hessian bound 

Hess„6(W) = QabWr = {rhab + 0{r^-^''°))Wr 

we observe that a more precise estimate is available on the angular part |Va;u|^ = 
r^dVulg — |9ru|2) of the the gradient of u: 



Is 

c 



JM2^„ JM,„,2\M2r„ °' 



We will revisit this analog of the angular Morawetz estimate for solutions of the 
time-dependent wave and Schrodinger equations in Lemma 11.1, where we will 
prove an even stronger result. 

Remark 6.5. As Proposition 1.14 shows, it is not possible to remove the error term 
from Lemma 2.3 in general. However, under the non-trapping assumption it is 
possible. For examples of such results see [100], [66], [102], [13], [82]. 
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7. High energy limiting absorption for non-trapping metrics 

We are now ready to prove Theorem 1.15. This case resembles the local-in-time 
theory of Craig-Kappeler-Strauss [30] and Doi [40], and indeed our main tool here 
will be the positive commutator method applied to a certain pseudo-differential 
operator, exploiting the non-trapping hypothesis to ensure that the symbol of the 
pseudo-differential operator increases along geodesic flow. As we shall now be 
working in the high frequency setting, we will not need to take as much care with 
lower order terms as in previous sections. It will be convenient to use the scat- 
tering pseudo- differential calculus, which is an extension of the standard pseudo- 
differential calculus which keeps track of the decay of the symbol at infinity. We 
briefly summarize the relevant features of this calculus hcirc;, referring the reader 
to [30] for more complete details. (This material will not be used outside of this 
section.) 

For any m,l £ R, we define a symbol s : T*M — >• C o/ order {m, I) to be any 
smooth function obeying the bounds 

|vsvfs(x,oi<c„,/3(0""""(^>-'-'"'; 

the function s(a;,^) = is a typical example of such a symbol. Note that 

we assume that each derivative in x gains a power of (x) , in contrast to the standard 
symbol calculus in which no such gain is assumed. We let S"^'\M) denote the space 
of such symbols. Given any such symbol s G ^'"''(M), we can define an associated 
pseudo-differential operator S = Op(s) by the usual Kohn-Nirenberg quantization 
formula 

Op(a)u(x) := {2Tr)-" J e'^''-y^^ls{x,£_)u{y) dy d^. 

We sometimes denote s by a-{S) and refer to it as the symbol of S. Heuristically 
speaking, we have S = a{s){x, jVx)- We refer to the class of pseudo-diff'erential 
operators of order (m, I) as Also, if /i : R — > C is any spectral symbol of order 

m/2, the corresponding spectral multiplier h{H) is a pseudo-differential operator 
of order (m,0). In particular, (1 -|- ff)"'/^ has order (m,0), and the Littlewood- 
Paley type operators Pio, Pmed, Phi have order (0,0). We caution however that the 
Schrodinger propagators e~**^ are not pseudo-differential operators. 

The composition of an operator S = Op(s) of order (m, /) with an operator of 
B = Op(6) order (m', /') is an operator SB of order (m + m', I + Z'), whose symbol 
a{SB) is equal to a(S)a{B) plus an error of order (m + to' — 1, Z + Z' + 1); note the 
additional gain of 1 in the decay index I, which is not present in the classical calculus. 
Similarly, the commutator i[S, B] will be an operator of order {m + m' — 1, — 1) 
with symbol <j{i[S, B]) equal to the Poisson bracket 

{a{S),aiB)} V,a(5) • V^aiB) - V^aiS) ■ V^aiB), 

plus an error of order {m + m' — 2,1 + I' + 2). We shall write the above facts 
schematically as 

a{AB) = a{A)a{B) + ^m-Fm'-i.i-Fr-M. 
a{i[A,B]) = {a{A),a{B)} + 5m+m'-2,;+i'+2 
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or equivalently as 

Op(s) Op(6) = Op(s6) + q,rn+m' -1,1+1' +1^ 
Z[0p(s), Op(6)] = Op({s, b}) + ^^,^+m'-2,l+l'+2_ 

In particular, since H has order (2, 0) and has principal symbol ^[^[^(j.) plus lower 
order terms of order (1, 1) and (0, 2), we see that if a € S""'', then we have 

i[H, Op(s)] =0p(Xs) + *™''+2, 

where Xa denotes the derivative of s along geodesic flow in the cotangent bundle 
T*M. 

Associated with the scattering calculus are the weighted Sobolev spaces if"''^(M) 
defined (for instance) by 



(M) 



(many other equivalent expressions for this norm exist, of course); when ^ = 
this corresponds to the usual Sobolev space H"^(M) and for m = 0, 1, 2 the space 
H"^'^{M) coincides with the weighted Sobolev spaces previously defined in (7). It 
is easy to verify that a scattering pseudo-differential operator of order (m, I) maps 
iJ-™'.''(M) to if™'-™.''+'(M) for any m',l'. 

In [30] (see also [40]) it was shown that the non-trapping hypothesis on M allows 
one to constriict a real- valued symbol s G S'^'" (depending on ip) which was non- 
decreasing along geodesic flow, Xs > 0, and in fact obeyed the more quantitative 
estimate 

Xs{x,0=^{xml + \b\'' 
for some symbol b of order (1, 1/2 -|- cr). The function (f){x) belongs to the class 
go,i+2a g^jj^ ^g^^ chosen to obey a lower bound 

0(x) > 

In Euclidean space, an example of such a symbol a is {x)-a{x)^-^'' ' ^ ^'^^ some 
sufHciently large constant C'^ and a small constant a. Quantizing this, we obtain 



1,2-|-2(T 



i[H, S] = VV(a;)Vj + B*B + ^ 

where S := Op(s) is a symbol of order (1,0), and B := Op(6) is a symbol of 
order (1,1/2 — a). Applying (63), discarding the positive term B*B and using that 
-^1^^+'^" maps iJi.-i/2-a to ^0,-3/2-a obtain 

/ V|Vu[^ dg < {\\Su\\H0.-i/2-a(M) + \\S*u\\H0.-i/2-a(M)) \\F\\ho.i/2+.(^m) 

Jm 

+ £||'S'w||l2(M)II"IIl2(m) + ||VM||/i-o,-i/2-<,(M)ll'"ll_f/o.-3/2-<'(M) 

Since S is of order (1, 0), it maps ii-i.-i/2-<^ ^ ^o,-i/2-<7 ^^^^ H\M) L^{M). 
Recalling that F = f + Vu with V satisfying the bound \V{x)\ < T{x)-'^-'^'"> and 
^{x) > {x)~^~'^'^, we have 

||Vu||^o,-i/2-<,(M) ^ (|I/IIho,i/2+„(m) + l|w||jfO,-i/2-<,(M) +£||w||jfi(M)||w||z,2(jv^)j 
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The term £|| VM||L2(jvf) ||w||l2(jv^) can be controlled with the help of (35), while the 
term £||w||i2(jvf) can be bounded from the charge estimate (28). As a consequence, 

||V'u||ho,-i/2-<7(m) < C'(^) (ll/llfJ'0.i/2+^(M) + II**IIho.-i/2-^(m)) • 

To prove the desired result it would be sufficient to show that 

||u||/j-o,-i/2-<,(M) < C (||Vu||/j-o,-i/2-<,(M) + ||'w||/i-o>-3/2-<'(M)) + II/IIho.i/2+<^(m) 

with a sufficiently small constant c. This follows immediately from (34) with c = 
C(r)A~^/^ provided that A is sufficiently large. We now have that 

||Vu|l^o,-i/2-„(M) + a^I|w|Iho.-i/2--(m) < C'(r)||/||/i-0,l/2+<.(M)- 

The remaining estimate for the second derivative follows from the elliptic regularity 
estimate (30). 

The proof of Theorem 1.15 under an alternative assumption that the potential 
V satisfies |V^(a;)| + {x)\\/V{x)\g < A{x)~'^^° follows simply requires running the 
positive commutator argument with H = —Am + V in place of —Am- We omit 
the details. 

8. Carleman inequalities 

In this section we prove the unique continuation estimate in Proposition 2.6. Let 
the notation and assumptions be as in that proposition. 

The standard way to prove such unique continuation estimates is via Carleman 
inequalities, which are inequalities of the form 

(66) [ {t^\uf+t\S/uf)e'^''^ dg<C{K,w) [ \ffe'^''^dg 

Jk Jk 

for some suitable smooth weight function w : K ^ R, and some large real param- 
eter t (typically t > C{K,w){l + VA)). However, as observed in Burq [12], such 
inequalities are not available when K has a non-trivial topology, due to the fact 
that w can be forced to contain stationary points. However, it is still possible to 
obtain a two-weight Carleman inequality of the form 
(67) 

/ {t''\uf+t\Vuf){e^'^^+e^'^^)dg<C{K,wi,W2) [ {\{-AM-X)uf){e^'^'+e^'''')dg, 
Jk Jk 

the point being that wo can choose wi and W2 to have stationary points at different 
locations (and furthermore that W2 > wi at the stationary points of wi, and wi > 
W2 at the stationary points of ^2). Such an inequality will allow us to obtain the 
above proposition. 

We now turn to the details. We begin by reviewing the standard approach to 

Carleman inequalities in the literature; for a more detailed survey sec [93]. Let u 
be a solution to the Helmholtz equation (27) which is supported on the compact 
set K. We rewrite the equation in the form 

(-Am -X)u = F := f± ieu - Vu 
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If we multiply this equation by a weight e"" for some smooth real- valued w, we 
obtain 

{e'^i-AM - A)e-"')e'"u = e"'F 
We now split b'^Amg^^ into real and imaginary parts 

e'"(-AM - A)e-'^' := Rc(e«'(~AM - A)e-"') + ilmie^{-AM - A)e-'"), 
where Rc(A) ^nd Im(A) . Using the general identity 

\\M\lhm) = l|R'e(^>|| 

i2(M) + ||Im(A)t;||i2(;i^-) + {i[Iie{A),lm{A)]v,v) l2(^]^'^ 

we thus obtain 

||Re(e-(-AM - A)e-)e-«|li.(M) + ||Im(e-(-AM - A)e-)e-zi||i.(^) 
(68) +(i[Re(e"'(-AM-A)e-'"),Im(e'"(-AM-A)e-'")]e"'u,e'"«>i2(M) 

= I|e"'/||L2(M); 

this identity is closely related to (62). In particular we have the inequality 



(z[Re(e-(-AM - A)e— ), Im(e-(-AM - X)^-)]e^u, e^u) ^^m) < Ik'/lli^tM)- 

Since w was arbitrary, we may replace w by tw for some arbitrary real parameter 
t (which we shall think of as being large and positive) to obtain 

(i[Re(e*-(-AM-A)e-*-),Im(e*-(-AM-A)e-*"')]e*"'«,e*-t.)i2(M) < We*"" fWl^^M)- 

The strategy is then to select the weight w so that the commutator i[Re(e*'"(— Am — 
A)e~*"'), Im(e*"'(— Am — A)e~*™)] is positive definite^^, at least to top order in t, to 
obtain a useful inequality such as (66) . 

It is certainly possible to adapt the above scheme to prove a two-weight inequality 
such as (67), and this was essentially what was done in [12]. Let us however pursue a 
slightly different (though closely related) approach, based on the inequality (62) as a 
substitute for (68), to emphasize the fact that certain applications of the Carlcman 
method can be viewed as special cases of the general abc method^^. Substituting 
w by tw as before, the identity (62) becomes 

(69) 

t [ Hess„/3(w)a"(e*'"w)a''(e*"'M) dg 
Jk 

+^3 / Ressa0{w){d'^w){d^w)\ufe''''" dg = t f ]{Aljw)\u\''e''''" dg + l\\e''"F\ 
Jk Jk ^ 8 



2 



It is easiest to apply this inequality when w is strictly geodesically convex (i.e. 
Hess(t(;) > 0) and non-stationary (i.e. Vw ^ 0) on the support K of u, as the 
left-hand side is then positive, and the first term on the right-hand side can be 



18a 

ctually, one only needs this commutator to be positive definite on the region of phase space 
where the operators Rc(e*"'(-AM - A)e-*'") and Im(e*"'(-AM - A)e-*'") vanish, thanks to the 
other two terms in the identity (68). 

"'^'^One advantage of doing so is that the abc method does not require more than one degree 
of regularity on the metric; in particular, the Riemann curvature tensor does not appear here, 
whereas this tensor will appear when computing the above commutator. 
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absorbed into the second term on the left-hand side if t is large enough. One would 
then obtain a Carleman inequality of the form 

i||V(e*-n)||i.(^) +t3||e*-«||i.(^) < CKi^)||e*-F||i.(^) 

for sufBciently large t, where the key point is that the constant C{w,K) is in- 
dependent of t. For instance, in the Euclidean setting M = R" one could take 
w = {x) = {1 + \x\'^y^^ to obtain an estimate of this form. 

For more general manifolds M , however, it is not possible to find a weight function 
w which is both geodesically convex and non-stationary, for two reasons. Firstly, 
if K contains a closed geodesic^", then it is clearly impossible for w to be strictly 
geodesically convex on this geodesic. Secondly, if K contains a non-trivial topology 
(e.g. K contains a handle), then from Morse theory we see that w must contain 
at least one stationary point. As mentioned earlier, the latter difficulty will be 
overcome by considering a pair of weight functions wi , W2 rather than a single 
weight w. To overcome the former difficulty, we can replace the weight functions 
'Wi,W2 by h{wi), h{w2) for "sufficiently convex" functions /i : R — )• R, however 
this only gives us convexity of the Hessian in the gradient directions Vwi, Vw2 
respectively (note though that this is the most important direction that we need 
convexity in, as one can already see from the second term on the left-hand side of 
(69)). To deal with the possible failure of convexity in the non-gradient directions, 
we will use energy identity (58). More precisely, we have 

Lemma 8.1 (Preliminary Carleman inequality). Let w : K ^ R be a smooth 
function and let B c K be a open subset of K. Suppose that w is non-stationary 
in the set K\B, so that the unit vector field h" := {V"'w)/\V°'w\g is well defined. 
Suppose that w is convex in the direction of ri" , and m,ore precisely 

Hess„„('u;) := HesSa/3 (w^)n"n^ > on K\B 

for some c > 0. Suppose also that for any unit vector field we have the estimate 

IlesSap{w)X"X^ > -^Hess„„(w;)|X|2 on K\B. 

Then for any solution u G H^{K) to the Helmholtz equation (27) which is supported 
in K, and any t > C{K, B, w){l + \/A) we have 

[ {t^\u\^+t\'^u\^)e^*'" dg<C{K,B,w)( [ iFl^e^'"" dg + [ {t^\u\^ + t\Vu\^)e^*'" dg 
Jk \Jk Jb 



'^^It is natural to conjecture that this implication can be reversed, i.e. that if K was geodesi- 
cally non-trapping and topologically trivial then there exists a weight function w which is geodesi- 
cally convex and non-stationary. This would allow one to replace the pseudo-differential calculus 
considerations in the previous section by a more elementary integration by parts argument. Unfor- 
tunately, there exist manifolds which are geodesically non-trapping, but for which no such weight 
function w exists; this can be seen by a minor modification of the example in [50, Section 6]. 
The point is that geodesic convexity is equivalent to the assertion that the function ^ av It? IS 
increasing with respect to geodesic flow on the cotangent bundle T* K = {{x,^a) '■ € T*M}, 
but this function ^aV"«; is necessarily linear in ^, which places additional constraints on solvabil- 
ity beyond the mere non-existence of trapped geodesies. If one uses pseudo-differential operators 
instead then one does not have this geometrically unnatural linearity constraint. 
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Proof. Suppose that we are in the region K\B. Then by hypothesis, we have 
Hess«/3H(5"u;)(5'^u;)|u|2e2*'" =Hess„„(M;)|Vw|^|tt|2e2*"' > 

and 

liesSa,fi{w)d^{e'^u)dP{e'^u) > -^Hess„„H|V(e*'"u)|2 
and hence 

iHess„^(w)a"(e*"'M)a'5(e*"'u) + t^Yi.esSoci3{w){d''w){d^w)\u\'^e'^^'" 

^5 

Integrating this on K\B and using (69), we see that 



\e j^^JiesSnn{w)\Vw\y''"\u\'' dg < Ct JjAlw)\u\'e^'^ dg + C\\e'^ F\\l,^^-^ 

+ U f EesSnn{wMw)\Vu\y^ dg 

+ C{K,B,w) f (t3|u|2 + t|Vu|2)e2*- dg, 
Jb 

where 93 is a cutoff function that equals 1 on K\B and vanishes near the stationary 
points of w. But if we apply the energy identity (58) with x = Hess„„(w)i^('u;)e2*"', 
and using the Cauchy-Schwarz inequality Re(uF) < t\u\^ + t~'^\F\^, we obtain 

\t ( B.esSnn{w)ip{w)\Vu\le'^*'" dg <\\t j Y^esSnn{w)ip(w)\u\'^ e^*"" dg 
^ Jk ^ Jk 

[ BesSnn{wMw)\Vw\l\u\^e^'^ dg 
3 Jk 



+ C{K,B,w)f [ I 
Jk 

+ \\<^^''''F\\hiMy 



\ufe''"' dg 



The contribution of the first two terms on K\B can be bounded by 
3 



/ Hess„„(«;)|V«;|2e2*-|«|2 dg 
Jk\b 

term since t"^ is assumed to be large compared with A. Thus we have 

[ Hess„„(«;)|Vu;|^e2*™|^.|2 dg < C\\e'- F\\l..^. 
o Jk\b 

+ C{K,B,w)t'^ [ lul^e^**" dg 
Jk 

+ C{K,B,w) [ {t^\u\^ + t\\/u\^)e^''" dg. 

B 



(absorbing the A|^w term into the C{K, B,w)t'^ Jj^ \u\'^e'^*'^ dg error). But since 
Hess„„(u')| Vw|g is non-zero on the compact set K\B, and t is assumed large com- 
pared with C{K,B,'w), we can absorb the C{K,B,w)t^ fj^ I'Upe^*'" dg error into 
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the right-hand side (plus the error on B) and eonclude that 

j e2*»|u|2d£,< C(ii:,B,«;)||e*'"i^|||2(M)+C(^,S,H j {t^\u\^+t\Vu\'^)e'^'''^ dg. 

Applying (58) again (with x = e^**", and using the assumption that t is large 
compared with \/A) we obtain the result. □ 

The above lemma has the drawback that there is an error term involving u on the 
right hand side. However, it is localized to a smaller set B than K. We can exploit 
this localization by using two weights instead of one, whose critical points are at 
different locations. We begin with a Morse theory lemma (first observed by Burq 

[12]): 

Lemma 8.2 (Construction of Morse function pair). There exist smooth functions 
ai : K ^ a and 02 : -ft' — )■ R which each have finitely many critical points on 

K, which are all in the interior of K . Furthermore, whenever x is a critical point 
of ai we have 02 (x) > ai{x), and whenever x is a critical point 0/02 we have 
ai{x) > a2{x). (In particular, the critical points of ai are at distinct locations from 
the critical points of a2)- 



Proof. By enlarging K if necessary we may assume that the boundary dK of K is 
a sphere dK = {{r,y) : r = R} for some i? ^ 1; in particular, K is now connected 
(since M is connected by hypothesis). We construct ai to be the solution to the 
Dirichlet problem 

Am^i = 1 on if ; ai = on dK 

which can be constructed for instance by a standard variational procedure (or the 
spectral theory of the Laplacian on a compact manifold with Dirichlet boundary 
conditions). By standard elliptic theory, this function is smooth on K, and all crit- 
ical points lie in the interior of K and are non-degenerate (since Hess(ai) is clearly 
non- vanishing), so in particular the number of critical points is finite. Furthermore 
none of the critical points can be local maxima, since AmOi is positive. In particu- 
lar if we enumerate the critical points of ai as Xi, . . . , Xm, and let £ > be a radius 
so small that the closed balls B{xi,e), . . . ,B{xm,s) are disjoint from each other 
and from the boundary dK, then we can find x'l,. . . , x'^ such that dM{xi, x'^) < e 
and ai{x'^) > ai{xi). 

We now \ct (j) : K ^ K he a. diffeomorphism which is equal to the identity outside 
of the balls B{xi,s), . . . , B{xm,£), and which swaps Xi and x'^ for each i. If we then 
set a2 := a\o (j) then it is clear that 02 has critical points precisely at x'^,... ,x'^, 
and the claims follow. □ 



We can now combine these two lemmas to obtain a two- weight Carleman inequality 
which avoids the B errors. 

Corollary 8.3 (Two weight Carleman inequality). There exist smooth functions 
wi : K and W2 ■ K ^ R with the property that for any solution u e H'^{K) to 
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the Helmholtz equation (27) which is supported in K , and any t > C{K, wi, W2)(l- 
\/A) we have 



[ {t^\u\^ + t\yu\l)ie^''"^ + e^'^-^') dg < CiK,wi,W2) [ \FWe^'^' + e^''^^') dg. 
Jk Jk 



Proof. Let ai, 02 be constructed by the previous lemma, and let a;i, . . . , Xm be the 
critical points of a\ and . . . , x'^, be the critical points of 02- If e = e{K, a\, 02) > 
is a sufficiently small radius, then we have 02 > ai on the set 

Bi := B{xi,s)U ...UB{xm,e) C K 

and ai > 02 on the set 

B2 := B{x[,e) U . . . U B{x'^,,e) C K 
for some c = c{K, ai, a^) > 0. 

We now let A = A{K, ai,a2,s) » 1 be a large parameter to be chosen later, and 

set wi := cxp(Aai) and u>2 := cxp(j4a2). Observe that for j = 1,2, wj has no 
critical points outside of Bj. Furthermore, on K\Bj we compute 

n,":=V"«;,V|Vu;,-|g=V"a/|Va|s 

and 

HesSc(/3(wj) = {A^S/aajV fiUj + AResSai3{a)) ex.p{Aaj) 
so in particular 

'ResSnjnji'Wj) = {A'^lVajfg + AHesSnj„^{aj)) exp{Aaj) 

and 

Bessxx{wj) = {A'^\X°'Vaajf + AHessxx(aj)) exp(Aoj). 
Thus if A is large enough, the hypotheses of Lemma 8.1 are obeyed, and we have 

/ (i3|u|2+t|Vu|2)e2*"'^- dg<C{K,Bi,B2,wi,W2)x 
for j — 1,2. Adding the two inequalities together we obtain 

j^{t^\uf + t\Vuf){e^'^' + e^'^') dg<CiK,BuB2,wi,W2)(^j \F\\e'''^- + e^'^') dg 

+ [ (t3|u|2+t|Vwp)e2*^i dg 

+ j (i^lwp + tlVuHe"*'"^ dg^. 

But observe that W2 > wi + c on Bi and wi > W2 + c on B2 for some c = 
c{K,wi,W2,Bi,B2) > 0, and thus for j = 1,2 we have e^*'^^ < e-2ct(e2t«'i +e2i«>2) 
on Bj . Thus if t is large enough, we can absorb the last two terms on the right-hand 
side into the left-hand side, and the claim follows. □ 



F|2e2*». dg+ I (t3|u|2+f|Vwp)e2*«'^ dg 

Bi 
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Remark 8.4. Morally speaking, the above two-weight Carleman inequality can be 
viewed heuristically as a single-weight Carleman inequality on the product mani- 
fold M X M, with u replaced by the tensor product u(Si u and w replaced by the 
tensor sum wi © W2 (and A replaced by 2A) . The point is that the critical points 
{xi, x'j) of Wi ® W2 lie off the diagonal, and the contribution of the weights at those 
points can be dominated by the contributions of the weights at the diagonal points 
(xijXi) and {Xj,Xj), for instance by the Cauchy-Schwarz inequality. Intriguingly, 
this perspective shares many similarities with the philosophy underying the interac- 
tion Morawetz inequalities, used for instance in [52] , where the positive commutator 
method of the previous section was also extended to a product setting. 

We can now quickly prove Proposition 2.6. 

Proof of Proposition 2. 6. Let x be a smooth cutoff function which equals 1 on K' 
and is supported on K . Then \u is supported on K and obeys the equation 

(-Am - X){xu) = -2(V„x)V"u - (Amx)m + x{f ± ^eu - Vu) + xVu. 

In particular, (—Am — A)(xu) is less than x(l/l + 0(^|m|)) on K' , and obeys the 
bound 

(-Am - \){xu) < C{x){A\u\ + |V«|g + |/|) 
on K\K' . Applying Corollary 8.3 to x^ we thus obtain 

/ (i^lup +i|Vu|2)(e2*«'i +e2*«'=) dg <C{K,wi,W2)C{x) I |/|2(e2*^i + e^**"^) dg 
JK Jk 

+ AC{x) f |w|^(e2*"'i +62*^"^) d5 

+ Cix) I (!«!' + |V«|^)(e2*-^ + e2*-=) dg 
Jk\k' 

for t > C{K,wi,W2)i'i- + \/A), and the claim follows by choosing f to be a large 

multiple of C{A, K, x, u>i, u>2)(1-I--\/A) (to absorb the second term on the right-hand 
side) and noting that wi and W2 are smooth and hence bounded above and below 
on the compact set ii'. □ 

Remark 8.5. Proposition 2.6 gives good control on solutions to the Helmholtz equa- 
tion on a compact set K. To obtain a limiting absorption principle, we will have 
to combine this proposition with more "global" estimates, such as the Morawetz 
estimates of the previous section, or the Bessel ODE analysis in Section 9. 

9. Conservation laws and differential inequalities of Bessel type 

The results of Lemma 2.3 and Proposition 2.6 provide us with the following very 
useful dichotomy. On one hand Lemma 2.3 shows that if the limiting absorption 
principle can be proved for the restriction, to a certain dyadic region ro/2 < (a;) < 
2ro, of a solution u of the Helmholtz equation, then it also holds on the set (x) > 2ro. 
Control of the region ro/2 < (x) < 2ro together with the unique continuation 
principle of Proposition 2.6 also imply that the limiting absorption principle can 
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be extended to the set (x) < ro/2. On the other hand Proposition 2.6 allows us an 
alternative scenario. To prove the limiting absorption principle on a compact set K, 
which can be thought of as the region {x) < ro/2, rather than proving unconditional 
control on the solution in the region ro/2 < (x) < 2ro it would be sufficient to show 
instead that the solution varies super-exponentially through that region, i.e., its 
rate of change is given by g-Cii+Vx) ^j^jj sufficiently large constant C. 

To show that either of these scenarios must be realized we need to perform analysis 
near the asymptotic end of the manifold M; this is the purpose of Lemma 2.8, 
which we shall prove in this section and in the next. 

Wo begin with an informal discussion. For simplicity consider the case when M 
is asymptotically Euclidean manifold. Then we can heuristically approximate the 
Helmholtz equation by its Euclidean version 

/ ^ n-1 d ^ 1 ^ ,_v + x]u--F. 

Assuming that V is radial, V = V{r), and applying the ansatz 

u{r,uj) = r-("-i)/2^(r)y;(a;), F{r,ij) = -r-^''+^^^^G{r)Yi{ij) 

where Yi is a spherical harmonic of order I on S"^^, normalized to have L'^{S^~^) 
norm equal to 1, the Helmholtz equation becomes the Bessel ordinary differential 
equation 

L{L- 1) 

Vrr o ^ ~ '^'^ + \v = G, 

r^ 

where L := 1 + As we shall see in Section 14 this equation will play an impor- 
tant role in a counterexample construction of Proposition 1.21. Despite providing 
good insight into behavior of solutions of the Helmholtz equation near infinity, the 
use of the Bessel equation approximation has several drawbacks. At first glance, 
it seems that this equation only emerges when the solution u has a specific struc- 
ture, namely that it decouples as the product of a radial function and a spherical 
harmonic. Of course, one could orthogonally decompose an arbitrary function u 
into such products and work on each harmonic separately (as is done in a number 
of places in the literature, e.g. [12]), but this becomes difficult if the metric and 
potential only decay slowly at infinity (although such analysis well suited for com- 
pact perturbations of the standard Euclidean metric). Also, such an approach often 
requires detailed analysis of the asymptotics of Bessel or Hankcl functions. Here, 
we present a more "energy-based" method to simulate differential equations (or 
differential inequalities) Bessel type for solutions to the Helmholtz equation (27), 
without requiring an explicit decomposition into spherical harmonics, and without 
requiring any knowledge of Bessel or Hankel functions (although such functions are 
in some sense lurking in the background in what follows). 

In this section we shall work purely in the asymptotic region r > i?o; thus the 
analysis here may be viewed as a "black box" analysis, requiring no knowledge of 
the manifold, solution, or potential in the interior region r < Rq. Eventually we 
will combine this black box analysis with the Carleman analysis in the near region 
r = 0(1) from previous sections to obtain the full limiting absorption principle. 
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52 n-1 d 4zdh[r] d 1 ^ . .r ■ . 

—u-\ —u + „r . -^u + -^Ah[r]U + Am = VuTieu- f. 

or r or an\r\ or r^ ' ' 



We write the metric in the forni 

g = dr'^ + r^hjk[r]{uj)du)^(Lj'' 

where 

f^3k[r]{i^) ■= hjk{uj) +r~'^'^''ejk{r,u) 
and work on the hypcrsurfaccs Sr ■= {{r,Lu) : uj € dM}, which are naturahy 
endowed with the metric h[r] and the corresponding measure dh[r] := ^yh[r]du; 
note that this differs from the induced measure dgls^ by a factor of r""^. In 
particular we have the co-area formula 

(70) / fdg=r ( f{r,uj)dh[r]{oj)r^-^dr. 

Jr>Ra JRo J 

We now rewrite the resolvent equation (26) in polar co-ordinate form as 

dh[r] dr 
Note that 

9jk{r,u>) = -drOjk{r,uj) 
is the second fundamental form of the surface Sr relative to the renormalized metric 

g = dr^ + hjkdij-' dui''. 

with mean curvature 

4-dh\r] 

It follows from (6) that 

(72) ^ = 0(r-i-2-o)^ |0|^j^j=O(^-i-2-xo)^ 

We expect solutions u of the Helmholtz equation to decay like r"^"'"-^)/^ as r — )■ 00. 
Thus we can renormalize u by defining 

(73) V := r("-i)/2u 

(cf. (124)) and observe that v obeys a Bessel-like equation 

(74) z;..+l(A,[,]- ~ ~ ^^ )^+A^; = -0Vr+iV+^e)vTiev+r^^-'y^f 
with the operator — A^[r] playing the role of the parameter 1(1 + n — 2). 

We now define the "spherical energies" 



Mass 


M[r] := 


= / \vf dh[r] 

JSr- 




Radial energy 


n[r] := 


-- [ \Vr\^ dh[r] 

JSr- 




Angular energy 


A[r] := 




-'f-'\vl^)dh[r. 


Mass flux 


^[r] := 


= / Ke{vvr) dh[r], 
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where VujV is the angular gradient. Note that the quantity M.[r\ was already 
introduced in Lemma 2.8. We also need the "forcing term" 

(75) G[r] := r^-^^/^ {\v\ + -J^^){\f\ + e\u\) dh[T], 

and record the following three "equations of motion" for the spherical energies, 
which arc closely related to the conservation laws (48), (50), (59). 

Lemma 9.1 (Equations of motion). We have positivity properties 

(76) M[r],H[r],A[r]>0 
the Cauchy-Schwarz estimate 

(77) \J^[r]\ < M[r]^/^n[r]^/^ 
and the equations of motion 

(78) = 2J^ + 0{r-^-^''°)M 

(79) J" = TZ + A- XM+ 0(r-^-^'"')M + 0(r-^-''°)X^/'^M + 0(G) 
dr 

^{'R + XM-A) = -A + 0{r-^-^^"){A + 'R + XM) 
dr r 

(80) +0(r-3-2'^°)A^ + 0(r-i+Ai/2)G. 
Here the implicit constants are allowed to depend on M and A. 



Proof. While in dimensions n > 3 positivity of the spherical energies and the 
Cauchy-Schwarz estimate are obvious, equations (78)- (80) follow from the iden- 
tities: 

-^M'^ = 2Re{vrv), 
dr 

d , , ,2 1 / A (^^ ~ - 3)x _ , , ,2 n _ 

— {v^v) = \Vr\ -:^[^h[r] 1 )VV-X\V\ - OVrV 

+ ({V + ^^-^e)v T iev + r("-^)/V)^^, 
d n ,2 M ,2 1,^ ,2 l)('^-3), ,2x 2/l_ ,„ (n - l)(n - 3) , 

-{\Vr? + X\V\' - ^\V^V\1 - ^ ^5 >-\v\^) = - [j^\V^v\l + ^ ^1 

+ ^diva;Re(Va,'y W) - 2e\vr\^ - ^^jfcV>V^w 

+ 2Re((y + ^^-j^O)v T lev + r("-i)/V)«r, 

identity (71) and the assumptions (72) and \V{x)\ < A{r-'^-'^'^° +X^/'^r-'^-'^'^°). □ 

Remark 9.2. It is helpful to keep in mind the model case (see (124)), where v solves 
the Bessel differential equation 

L(L-l) 

Vrr ^ — 2 + At; = 0, 
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in which case 

M[r] = |,;(r)|2; 7^[r] = \vr{r)\'; A[r] = ^^i^^\v{r)\^; 
T{r) = ^e(^Vr{r)); G[r] = 0. 

The reader may wish to verify the above equations of motion (with all error terms 
set to zero) in this special case. It may also be useful to keep in mind the dimensional 
analysis 

r ~ length^; M ~ length°; T ~ length"^; Tl, A, G ~ length"^, 

noting that the above equations then become dimensionally consistent up to errors 
involving ctq- 

Remark 9.3. Note that while we have four energies, we only have three equations 
of motion; we do not control the evolution of TZ and A separately, but only have an 
equation for a certain combination TZ — Aoi these two. However, we can obtain a 
lower bound on TZ from (77). This system of three equations and one inequality is 
thus still underdetermined, but we will still be able to extract enough control out 
of this system to establish all the estimates we need. 

Remark 9.4. Of course, the three equations of motion can also be interpreted in 
terms of the Friedrichs abc method, where the cutoff x is now the surface measure 
on a sphere Sr- We omit the details. 



Now we obtain some preliminary estimates on the above energies, in the setting of 
Lemma 2.8. 

Proposition 9.5 (Preliminary estimates). Let the notation and assumptions he as 
in Lemma 2.8. Then we have the integral estimate 

/•oo 

(81) / G[r] dr = 0{6). 

jRo 

and the boundary condition 

1 f^^ 

(82) lim - / M[s] + n[s] + A[s\ ds = 0. 



Proof. The boundary condition (82) follows from the normalization ||M||//o,-i/2-<T(jvf') = 
1 and (31). To prove (81), we first see from (75), (124) that 

G[r] < Cr"-i i\u\ + -i^!^)(|/| + e\u\) dh[r] 

and so by dyadic decomposition it would suffice to show that 

(83) E / (\u\+ }^^\ ){\f\+e\u\)dg = 0{S). 

From (41) we have the bounds 

/ \u\^dg<CRR'+^'^ 

Jr<{x)<2R 
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and 

\f?dg<5''cRR-^-^'^ 



L 



lR<{x)<2R 

where cr > Q are numbers such that X^/jCk = 0{1). From the charge estimate 
(Lemma 2.1) we also have 

£ / |«|^ dg < ||/||ho,i/2+<,(m)||'"||ho,-i/2-<,(m) < <^ 
J M 

and thus 



L 



'R<{x)<2R 

(after adjusting cr if necessary). Finally, from (31) we have 

2 



Jr<{x)<2R\ -K +A/ / Jr/2<{x)<4R^ (A + K ) / 

i so after adjusting cr a bit more we obtain 

+ iF^TAV^)' ^ (min(i^^+^^ V^) + 

(I/I + e|u|)2 dff < Ccr{S^R-^-^^ + mm{e5, ke^R^+^^)). 



L 



iR<{x)<2R. 

On the other hand, we have 



L 



'R<{x)<2R 

Applying Cauchy-Schwarz, we will obtain (83) if we can show that 

(^min(i?i+2-, 5/e) + (^'^^l'^^ ) iS^R~'~^^ + min(£^, ke^R'^^')) < CS^ 
for all R > Ro- This is clear for the terms involving mm{R^^^"' , 6 / e) . For the term 

^2^-l-2a 



rS'R- 



-1-2(7 



(A + i?-2)2 

we bound A + R~^ from below by i?"^ iq obtain a bound of S^R^~'^'' , which is 
acceptable since ct < 1/2 and 5 <C. Finally, for the term 

Si2 p-l-2<7 

A^^min(.<5,fc£2^i+2.^ 

we bound A + R~^ from below by A, and use the second term in the minimum, to 
obtain a bound of e^S^ j}? , which is acceptable since e < A. □ 



10. An ode lemma 



In view of Lemma 9.1 and Proposition 9.5, we see that Lemma 2.8 will follow 
immediately from the following ODE lemma. 

Lemma 10.1 (ODE Lemma). Let C\ ^ Rq be a large numher. and then let C2 3> 
Ci he an even larger numher. For all r > Rq, let Ai, TZ, A, T , G he real-valued 
functions obeying the differential inequalities in Lemma 9.1 as well as the properties 
(81), (82). Then if Ci is sufficiently large (but not depending on \), and C2 is 
sufficiently large depending on Ci (but not on X), then one of the following must 
he true: 
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• (Boundedness) There exists a radius C'l < ro < C(Ci,C2) such that (42) 
holds. 

• (Exponential growth) For all C\<r< lOCi, we have (43). 

Of course, the magnitude ofC\, C2, C{Ci,C2) will depend on the implicit constants 
in Lemma 9.1 and (81) (which in practice will depend on M and A). 

The proof of this Lemma is lengthy and will occupy the remainder of this section. 

10.2. Heuristics. We first describe in informal terms why one would expect there 

to be a dichotomy of the type asserted in Lemma 10.1. Wc first observe that (80) 
is an approximate monotonicity formula for the quantity TZ + AA^ — A. Since this 
quantity is zero at infinity by (82), we expect it to be negative (up to errors of size 
of 0{5), thanks to (81)) at other values of r. Thus we have a lower bound on A, 
hcuristically of the form 

(84) A>n + \M-0{6). 

This converts (79) into a monotonicity formula as well, roughly of the form 

(85) ^F>2n-0{5). 

ar 

Now the dichotomy in Lemma 10.1 rests on whether the forcing terms such as TZ or 
|>t in the monotonicity formulae are large enough to dominate the error terms such 
as 0{S). If this domination never occurs (or only occurs when r is relatively small), 
then one ends up in the "boundedness" scenario (42). On the other hand, if at least 
one of the forcing terms becomes large, one expects that this will eventually force the 
other forcing term to be large as well (as r decreases towards Ci), causing a positive 
feedback loop which will eventually lead to the "exponential growth" scenario (43). 
For instance, if TZ gets large (compared to A4 and 5), this should force M to be 
similarly large thanks to (85); from (78) one then expects A4 to grow exponentially 
(but slightly less fast than J^) as r decreases; from (77) one then expects TZ to stay 
large, thus creating a self-sustaining feedback loop. Similarly, if A gets large, then 
from (80) wc expect the quantity TZ + XA4 — ^ to get large and negative, which 
adds an additional positive term to the right-hand side of (85), which as mentioned 
earlier should cause J", M, and TZ to grow; using (84), this should eventually sustain 
the growth of A, thus creating another self-sustaining feedback loop. If these loops 
start far enough away from the origin (e.g. at r > C(Ci, C2)) then one might hope 
to expect the growth to become exponential with growth rate C2 by the time r 
reaches Ci, which is the "exponential growth" half of the dichotomy. 

Remark 10.3. An oversimplified model of this dichotomy can be seen by considering 
plane wave solutions to the Helmholtz equation Am = Am in the (flat) cylinder 
R/27rZ X R+ := {(6',r) : 6 G R/27rZ,r e R+}, which should be thought of as a 
caricature of polar co-ordinates; we assume some boundedness on u at r = +00 (e.g. 
u{r) — 0{S) for sufficiently large r) but not when r is small. One has "bounded" 
solutions of the form u{6,r) = Ce^"'^e'^ where a is an integer such that \a\^ < A 
and + 6^ = A and C = 0(S). Then there are "exponential growth" solutions of 
the form u{x,y) = Ce*°^e^''^ where |ap > A and a'^ — = X and C is arbitrarily 
large. Thus one expects the solution to stay bounded if the "angular energy" 
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stays smaller than A, and to grow exponentially otherwise. This dichotomy roughly 
corresponds in our setting to the case when A stays controlled hy TZ + XM (which 
will basically ensure the bounded scenario) or is larger than this quantity (which will 
ensure the exponential growth scenario). See [21] for some rigorous formulations of 
these heuristics, where the cylinder has now been replaced by a stadium. 



10.4. Step 1: a Pohozaev bound. We now begin the rigorous proof of Lemma 
10.1. We begin by giving a rigorous version of (84), based primarily on (80) and 
the positivity of A. 

Lemma 10.5 (Pohozaev bound). LetV denote the "Pohozaev flux" 

V[r] := XM[r]+Tl[r]-A[r]. 
Then for all r>C\, we have 

V[r] < 0{{r-^ + A^/2))(5 + 0{r-^-^'')M[r]. 

Proof. We rewrite (80) in terms of the Pohozaev flux as 

= + 0(r-i-2-o)(2yl + V) + 0(r-3-2-o)7Vl + 0{{r-^ + X^^^))G, 

To eliminate the 0{r~^~'^'^°)M. error^^ we shall consider the modifled Pohozaev 
flux 

r* ■.= r- Cor-^-^'">M 
for constant Co » 1, and observe using (78) that 

^r* = -A + 0{r-^-^'"'){2A + r* + Cor-'^-^''° M) 
dr r 

+ Cor-^-^"" {T + 0(r-^-2^° )M) 

+ (2 + 2a)Cor-^-'^'''>M + 0{r-^-'^''°)M + 0(r-i + X^/'^)G. 
We can use (77) to bound 

|jr| < j.-i-'^o^ r^+^o-jz = r-i-'^o(l + Co)M + 2r^+''"A + r^+^^V* . 

If Co is suitably large (and Ci < r is also suitably large) then the net M. term on 
the right-hand side is positive, as is the net A term. Thus we have 

^V* > 0(Cor-i-^»)P* + 0(r-i + X^^^)G. 
dr 

On the other hand, from (82) we know that 'P*[r] ^ as r — >■ oo. The claim then 
follows from (81) and Gronwall's inequality. □ 



^^In dimensions n > 4 we can use the A term to control tlie A4 error since {n— l)(n — 3)/4 > 
in that case; this leads to some minor simplifications in the proof of Lemma 10.1. 
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10.6. Step 2: Dimensionless formulation. In order to analyze our system lur- 
thcr it is convenient to make a number of changes of variable to a more scale- 
invariant or "dimensionless" formulation. At present we have three equations of 
motion (78), (79), (80) and one inequality (77) (as well as the positivity properties 
(76)) for four unknowns. We shall now use the inequality (77) to replace two of the 
unknowns 7?., ^ by a single unknown V, at the cost of replacing the equalities in 
(79), (80) by inequalities. Indeed, we can rewrite (79) using (77) as 

=2n-r + 0{r-^-^''")M + 0{X^/'^r-^-^''°)M + 0(G) 

- ^5 ~ ^ ^ 0(r-2-2'^°)A1 + 0{X^/^r-^-^''")M + 0(G). 
As for (80), we rewrite it (for Ci sufficiently large) as 

^P=Ia + 0(r-^-^^°)(2A + V) + 0(r-3-2-o)7W + 0(r-^ + X^/^)G 
dr r 

and then write 

A = n + XM-V> — + XM-V. 
M 

We thus have the new equations of motion 

-^7W =2J- + 0(r-i-2''°)A^ 
dr 

^J" >2^-V + 0(r-^-^''°)M + 0(X^/^r-^-^''°)M + 0(G) 
dr M 

^r>(- + O(r-i-2ao)\ (^^ + xM-V) + 0{r-^-^^°)V 
dr \r J M 

+ Oir-^-^"")]^ + 0{r-^ + X^'^)G. 

To analyze these equations, we now adjust the quantities J^, V slightly to handle 
the forcing terms involving G. Define 

j:*[r] := T[r] -C \G[s]\ ds = J^[r] + 0{6) 

J r 

and 

(s-1 + X^/^)\G[s]\ ds = V[r] + 0{r-^ + X^l^)5, 

where we have used (81). If the constant C appearing above is large enough, then we 
can dominate the G[r] forcing terms on the right-hand sides of the above equations 
to obtain 

^M = 2F + 0{r-^-'^'''')M 
dr 

^* > 2^ - P + 0(r-2-2<^o)A^ + 0{X^'^r-^-^'"')M 
dr M 



Ip* > + 0(r-i-2-°)^ + - P) + 0{r-^-^''°)V + 0{r-^-^''«)M. 
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Writing T = T* ^ 0{S) and so J^^ > {J^*Y + 0{5)F* , and similarly writing 
p = + 0(r-i + AV2)5, we thus obtain 



(86) 
d 

(87) 



= 2T* + 0(r~i-2'^°)7W + 0{5) 

ar 



d^. > ^ iJ"*) \0{5)T* _ ^ ^ ^ ;^l/2^)^-2-2.o_^ + + ;^l/2)^ 

dr M 



(89) + 0{r-^-^"°)M + 0{r-^ + X^/^r-^)6. 

To analyze this system of equations, it is convenient to work in the "dimensionless" 
co-ordinates^^ 

fj,{r) := r 



M[r] 

a{r) := -r-r-^ 
M[r] 

P{r) := —r 



,2^*M 



M[r] 

and to introduce the "dimensionless" derivative D := —r-^. We then have 

ar 



= -li- 2an + 0(r-2'^°)M + O(m^) 

and 

dr dr 

2 

> -a + (2(jr*)2 + 0(5)jr* - A^p* + 0(1 + rAV2-,^-2-2ao;^2 
= l3-a + 0{ix)a + 0(1 + rA^/^^r'^''" + 0(1 + Ai/V)At + 0(r-2'^°)c 



^^Admittedly, tlicsc co-ordinates have a singularity when M[r\ = 0, but this will not be 
relevant for us as we shall only perform the remainder of the analysis in the case when fj, is small 
(and hence M is large). 
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and 

Dl3 = -2p+^{M^V* -V*^M) 
dr dr 

> -2/3+-^ {^ + 0{r-^-^''°){{F*f + 0{5)F* + XM'^-MV*) 

+0{r-^-'^'"')MV* + 0{r-^-'^'"')M^ + 0{r-^ + X^/'^r-'^)5M 
-2V*F* + 0{r-'^-'^''°)MV* + 0{5)V*) 
= -2/3 + (2 + 0(r-2<^°))(a2 + 0{ii)a + Ar^ + /3) + 0{r-'^''°)l3 + 0(r-2<^°) 
+ 0(1 + \^/'^r)n - 2a/3 + 0(r-2<^°)/3 + 0(m)/3 

> 2a{a - /3) + Ar^ + 0(/^(l + A^V + |a| + |/3|)) + 0(r-2<^« (1 + |^| + a^)). 

Meanwhile, from (10.5) we have 

V* < 0(r-i + \^'^)5 + 0{r-^-^'')M 

and hence 

P > -0(l + Ai/V)/i - 0(r-2'^"). 
To summarize, the functions a(r), ^(r), /Lt(r) obey the differential inequalities 
(90) 

Da>p-a- Oiiii + r-^''°){l + A^/^ + |a|)) 
(91) 

Dp > 2a{a - /3) + Ar^ - 0(/i(l + A^/V + la] + |/3|)) - 0(r-2'^° (1 + + a^)) 
(92) 

DiJL = -ii- 2an + Oir-^^^ii) + ©(/x^) 
(93) 

/3 > -0(1 + A1/V)m - 0(r-2-o 



10.7. More heuristics. Recall that we are trying to establish a dichotomy be- 
tween boundcdncss and exponential growth. In our new co-ordinates, boundcdncss 
roughly corresponds to fi being bounded away from zero for small values of r, while 
exponential growth corresponds to a being large (as can be seen either from (92) 
or (78)). So, hcuristically speaking, we have to rule out the scenario in which /i is 
small and a is also small. If fjt is very small, however, then we expect to be able to 
ignore most of the error terms in (90), (91), and thus we reduce (heuristically) to 
the model equations 

(94) Da>p-a; D/3 > 2a(a - /3) + Ar^; /3 > 0. 

The intuition here is that as r moves backwards from oo to Ci, the first inequality 
will lift up a if a is below /3, while the second inequality will lift up /3 if /3 is below 
a (note that this forces both a and /3 to be non- negative, by the third equation); 
indeed, when Ar^ is large there some additional lift applied to /3 (and thus indirectly 
to a, by the first equation). These equations then suggest that if a and /3 are both 
large at some radius r, then they will stay large for all smaller radii also; this is 
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what will lead to the exponential growth scenario. On the other hand, if a and /3 
stay small for all time then we can hope to obtain the boundedness scenario. 

10.8. Step 3: A condition for igniting exponential growth. We now make 

the above heuristics rigorous. 

Lemma 10.9 (Exponential growth is self-sustaining). Suppose that r > Ci is such 
that 

li{r) < 1/Ci and a{r) > C|. 
Then for all Ci < s < r we have 

(95) /i(s) < 1/Ci and a{s) > 2C2. 

Proof. From (92) we have the somewhat crude estimate 

(96) Dn{s) < — ^/i(s) whenever a{s) > and /x(s) < 1/Ci. 

From (96) and the continuity method it will thus suffice to show that 

a{s) > 2C2 for all Ci<s<r. 
Suppose this claim is false, then there exists Ci < s* < r such that 

(97) a(s*) = 2C2 
and 

(98) a{s) > 2C2 for all s* < s < r. 

From (96) and Gronwall's inequality (and the continuity method) we conclude in 

particular that 

(99) ^i(s) < {s/rf'^ii{r) < {s/r)^^^/Ci < 1/Ci for all < s < r. 
It is convenient to introduce the quantity 

k{s) := a^{s)+(3{s). 

From (90), (91) we have 
Dk{s) = 2a{s)Da(s) + D(3{s) 

> -Oiifiis) + s-^^^Xais) + Ai/2sa(s) + a^{s))) 

+ Xs^ - 0(m(s)(1 + A^/^^ + a{s) + |/3(s)|)) - 0(s-2-«(l + |/3(s)| + a^s))) 

> As^ + 0{{fi{s) + .s-2-o)(i + ^(5) + ^2(^) ^ ^1/2^ ^ X^/^sa{s) + \(3{s)\)) 

> XS^ - 0((m(s) + S-2-o)(l + „2(^) ^ ^^2 ^ l^(^-)l)) 

> (1 - 0(m(s) + s-2-''))As' - 0((m(s) + s-'<^°)(l + «(s) + 2min(-/3(s),0))). 
But from (93) we have 

min(-/3(s), 0) < C/x(l + As^) + Cs-^"^" 

and thus 

Dk{s) > (1 - 0((m + s-^'">){l + /i)))As2 - 0((m(s) + s-^^°){l + k{s) + n{s))). 
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Using (99), we conclude 

(100) Dk{s) > ^Xs^ - 0((Cf i(.s/r)i/2 + s-^''"){l + k{s))) for all < s < r. 

On the other hand, from (93) we have the crude estimate 

p{r) > -0(1 + A^/V) 

and hence 

K{r) > C-^Cl - CX^I'^r. 
From this, (100), and Gronwall's inequality we see that 

(101) k{s) > C-^{Xr{r - s) + - CX^^^r) for all < s < r. 

On the other hand, from (90) and (98), and writing /3 = At — a^, we have the rather 
crude estimate 

(102) Da{s) > k{s) - 0(Cf ^(1 + A^/^s)) - 0{a'^{s)) for all s* < s < r. 

From (97) we have a{s^) < 3C2 < a(r). Thus we can find < < s such that 
a(r») = 3C2 and 2C2 < a{s) < 'SC2 for all s* < s < r». Then by (101), (102) we 
have 

Da{s) > C-^Cl + C-^Xr{r - s) - CX^I'^r 

Observe that the expression on the right is positive unless r > C2A~^/^ and s = 
r + 0(Ci A^^/^), in which case it is bounded below by — CA^/^r. From this and the 
fundamental theorem of calculus we see that 

-C2 = a{n) - a{s.) > / Da{s) — > -Ci, 

Jr, S 

a contradiction. The claim follows. □ 

Corollary 10.10. There exists C3 = C(Ci,C2) » C2 such that if there exists 
r > C3 for which ii{r) < 1/C\ and a{r) > Cf are both true, then we are in the 
exponential growth scenario (43) . 

Proof. Suppose first that we are in the low energy case A < 1. Then from the above 
Lemma we have a{r) > 2C2 for all Ci < r < lOCi. From (92) we then have 

-DAt(s) < -2C2m(s) for all Ci < s < lOCi 

which by definition of /i(s), yields the mass growth estimate (43). 

Now suppose that we are in the high energy case A > 1. In this case we observe 
from (102), (101) (with r now being replaced by C3) that 

Da{s) > C-^(AC3(C3 - s) + C|) - CC3X^/'^ - 0{a^{s)) for all Ci < s < C3 

and in particular 

Da{s) > C-'^XCi for all Ci < s < C3/2 such that a{s) < C|A^/^ 

We also have a{s) > 2C2 in this region. We thus conclude (if C3 is large) that 

a(s) > C|A^/^ for all Ci < s < lOCi 

and then by arguing as before we obtain (43). □ 
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10.11. Step 4: The case of no exponential growth. In light of the above 
corollary, we may assume without loss of generality that for any r > C3, at least 
one of 

(103) ii{r) > 1/Ci or a(r) < C| 
is true. 

We can now remove the second half of of the dichotomy (103) at large distances. 

Lemma 10.12. Suppose that (103) holds. Let C4 = C(Ci, C2, C3) > C3 6e a 

sufficiently large constant depending on Ci, C2 and C3. Then we have fi{r) > 1/C| 
for a//r>C4(l + A-V2). 

Proof. By (82) we have ^{r) ^ 00 as r — >■ 00. Thus if the claim is false, then we 
can find ri > ro > C4(l + A"^/^) such that 

(104) ^i{ro) = 

(105) A*(n) = C^' 

(106) < Ai(r) < CJ^ for all ro < r < n. 
In particular from (103) we have 

(107) a{r) < C| for all ro < r < ri 
From (105) we have 

Diiin) < 

which by (92), (105) forces 

(108) a(ri) > -1 

(for instance). Also, from (92), (106), (107) we have 

£>M(r) > -2Clfx{r) 

which by Gronwall's inequality and (104), (105) forces a certain largeness bound in 
the interval [ro,ri]: 

Jro 

which implies that ri > ro C4 ^ . In particular, if C4 is large enough then we 

have 

(109) ri-ro > C4(l + A-i/2)^ 

(for instance) . Now we control /3. If ro < r < ri, thenr > C4,X~^^'^ and/Lt(r) < CJ^, 
which by (91) implies the crude bound 

Dp > C-^c? - 0((1 + |a|)(l + I^D) + C-^Ar^; 

using the crude bound 0((1 + |a|)(l + |/3|)) = 0(1 + C-^\a}^ + C|^p) we conclude 

DP > C-'^Xr^ whenever |/3(r)| < C-VA^/^ 

On the other hand, from (93) we have 

> —CX^^^r whenever ro < r < ri. 
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Combining these two equations we see that 

/3(r) > C^^A^/V whenever tq < r < n - CA"^/^ 
On the other hand, from (90) and the estimates r > Ci\~^/^, ji < C^^ we have 
Da>(3- 0{a) - (^(Q"^''" A^/^r), 

and thus we have 

Da{r) > C-^rX^/^ - 0(a(r)) whenever ro < r < n - CA"^/^ 

and 

Da{r) > -0{C\^l'^r) - 0{a{r)) whenever n - CA"^/^ < r < n. 
Moreover, since r > Ci\~^/'^ and a{r) < C|, we can replace the first estimate by 

Da{r) > C-VA^/^ whenever ro < r < n - CA"^/^ 

Prom these estimates and the initial condition (108), applying Gronwall's inequality, 
we see that 

o:{r) > —C whenever n — CA~^/^ <r<r\ 
and then by a further application of Gronwall we see that 

a(r) > C-i(ri - r)A^/^ - 0(1) whenever ro < r < n - CA'^/^ 

Since ri — ro > C4A~^/^ this contradicts (109) and (107), and the claim follows. □ 

In the high-energy case A > 1 this lemma immediately gives the boundedness half 

(42) of the dichotomy, by choosing r := 4C4 (for instance) and using the definition 
of ji. In the low-energy case one observes from (103) and (92) that 

D^i{r) > -2Cifi{r) 

whenever C3 < r < 6*4(1 + A~^/^) and /x(r) < 1/C'i. From this. Lemma 10.12, and 
Gronwall's inequality we see that 

(110) M(r) > (:7(Ci,C2,C3,C4)-U^(^i'^^) 

for all C3 < r < C4(1+AV2), and by definition of fj, we are again in the boundedness 
half (42) of the dichotomy. The claim follows. This completes the proof of Lemma 
10.1, and thus Lemma 2.8. ■ 



11. Additional equations of motion 

As we have already seen in Section 3.1, the combined results of Lemma 2.1, Lemma 
2.2, Lemma 2.3, Proposition 2.6, and Lemma 2.8 are already sufficient to prove 
Theorem 1.7. But to prove the Sommerfeld radiation conditions of Proposition 
1.12, it turns out that we must study further spherical energies in addition to the 
energies M., TZ, A, T introduced in Section 9. We take a quick detour to define these 
energies, study additional equations of motions and derive their consequences, in 
particular establishing Lemma 2.4. The results of this section will also be useful in 
the proof of Proposition 1.27. 
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Using the notations of Sc;c,tion 9, we define the following additional energies: 
Outgoing null energy 7V[r] := / \Vr T izvl"^ dh[r], 

Complex flux Z[r\ := I VrV dh[r], 

Jsr- 

where z = a±ih]s the complex number such that = \ ±ie and Im(^;) = ±6 > 0. 

Lemma 11.1 (Additional equation of motion). Let u G H^[M) he a solution to 
the resolvent equation 

(H -{X± ie))u = f 

for some A, £ > with the Hamiltonian H = —Am + V, satisfying assumptions of 
Proposition 1.12. Then (with the notation of the previous two sections) 

^{M - A) ^-A + 2b {M + A) + 0(r-2-2-o 
ar r 

(111) + o(r-i-2-o (7^ + ^ + kpM)) + OiFlrY'^N^I^) 

where F[r] := r"-i J^^ \f\'^dh[r]. 

Remark If. 2. A similar equation of motion appeared in [57] for the proof of a 
qualitative limiting absorption principle for a Schodinger operator with magnetic 
potential in R". This equation can be interpreted as a special case of the identity 
(61). 



Proof. For simplicity assume that {H — X — ie)u = f. Note that = Re.E. We 
have as in Lemma 9.1 

-^M = 2ReZ + 0(r-i-2'^°M), 
dr 



d_ 

dr 



Z = n + A-z'^M + 0{r-'^-'^''°)M+r'^ [ {f + Vv)vdh[r] 

JSr 



Note that 

\vr — izv\^ = + Izpl^'P + 2Im(z?;iv) 

= \vr\^ + z^\v\^ + {\z\^ - z^)\v\^ + 2lm{zvv^). 

Thus, 

J\f = n+\z\'^M + 21m{zZ). 

Since 

^(tZ-A)^ -A-2Re{z^Z) + 0{r-^-^^°)in + A) 
dr \ / r 

(112) +2 / Re((/ + yv)iv)d/i[r] 
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we derive the following equation for M — A: 
^(M-A^= ^^ + 2Re((|z|2-z2)Ar) +2Im(z)(7^ + ^) 

+ 2|z|2 luy{z)M + 0(r-l-2-«)(7^ + A+ \z\^M) + 0{F[rY/^N^'^) 
+ 0(r-2-2-o)im(^)A^ + 0(r-2-2-o + \^/\-^-^-^)M^/^N^/^ 

= + 2Re (d^l^ - z'^)Z) - 4Ini(^) lm{zZ) 
+ 2\m{z) {N + A) + 0(r-i-2'^°)(7e + A+ \z\'^M) + 0{r-'^-'^''°)\m{z)M 



Since 



we have 



2Re {{\z\^ - z^)Z) = 4Im(^) Re{-izZ) = 4Im{z) lm{zZ), 



= ^A + 2b {Af + A) + 0{r-^-^'">){'Jl + A + \z\^M) + 0{r-''-^''°)bM 

+ 0(F[r]l/2^V2) + O(r-2-2-0 + ;^l/2^-l-2ao);^l/2^1/2_ 

□ 



11.3. Proof of Lemma 2.4. We are now ready to establish Lemma 2.4. Let the 
notation and assumptions be as in that lemma. 

Multiplying equation (111) for (AA— A) by r^'^ we obtain 

l-r^^iAf _ ^) = _1_ ((2 - 2a)A + Af)+ 2br''' {Af + A) 

+ 0(r-i-2'^<'+2'^)(7e + A + \zfM) + r'^''0{F[r]^/^Af^/^) 

+ 0{r-^-^^°+^'')bM + O(r2^(r-2-2<^0 + ;^l/2^-l-2ao)-);\4l/2^1/2_ 

The presence of the 2b r^'^ (A^ + ^) term will allow us to ignore the boundary term 
at r = oo arising after integration in r. Therefore, for some universal positive 
constant c. 



/>00 /'OO 

c {r-^+^'' +br^^) {J^ + A)dr< r-^-^^°+^'' (TZ + \z\''M) dr 

J 2ro J ro 

p^ro POO POO 

(113) +r-o^+^'' / Adr+ r-^-^'^o+^'^A^ dr + / r^+'^T[r]dr. 

Jro Jro Jra 

Note that in view of Lemma 2.3, the bound (113) already gives the estimate (40) 
in the region s < CX. 

We now claim that u verifies the Poincare type inequality 
(114) 

+ b\\u\\lo,-.+. 

(M2„) - C'(''o) / Af[r] + ||w|'^ 

\J2ra 



L=(M2,„\M4,„) ) • 



Note that (114) together with (113) and Lemma 2.3 imply all the statements of 
Lemma 2.4 for s < CX. 
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To establish (114), observe that 

M[r] = [ \vr- izvfdh[r] = [ e-'^^\dr{e-"''v)\'^dh[r] 

Our assumptions on the metric g imply that the area form dh[r], which corresponds 
to the metric {habis^) + r~'^'^°eab{r, dcj^, is equivalent to the area form dh of 

the metric hab{'jj)duJ°' duj^. Integrating by parts the expression 

/°°dr r-2+2-((l-a)r-i +6) / \v\''dh 

J2ro J Sr- 

= r dr r-i+2'-e-2*'' / Re (drie'"'' v)(F^)) dh 

+ i(2ro)-2+2"((l - a)(2ro)-^ +b) j \v\^dh 
<2(l-a)-i/ drr-^+^^e-^'"' [ \dr{e-''''v)\'^dh 

J2ro JSr- 

+ h2ro)-'+'mi-a){2ro)-'+b) j \v 



S2ro 



and averaging over ro we immediately obtain (114). 



It remains to consider the case £ > CA. Under this assumption we have \z\ < 2e 
(for instance). In view of (114) and a trivial inequality 

the bound (113) can be reduced to the estimate 

+ C(ro) (||u||l2(m,„\M4,j + I|Vm||l2(m,„\M4,„) + II/II//o,i/2+<,(m,„)) • 

According to (29) 

e||u||fl-0,-l/2-,xo + <.(M2,^) < C|l/ll//t).-l/2-<'0+<'(M,o) + C'l|Vw|li/0>-3/2-<'0+<'(Mr(,)) 

which implies that 

||Vu||jJ-0,-l/2+<,(M2^lj) + ll'"ll/i-0>-3/2+<'(M2ro) - ^NIII'"ll/i'0.-l/2-<^0+<'(M2r(,) 



provided that ro is sufficiently large. The bounds (38), (39) follow immediately, 
while (40) can be recovered from (113). The proof of Lemma 2.4 is now complete. 
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12. The low energy regime 



We now prove Proposition 1.27. We assume that the potential V obeys the bounds 
\V{x)\<A{x)-^-^'^'^, [ |y_(a;)|t</3 

JM 

with a small constant /? = I3{M) and show that for all sufficiently small £, |A| < 
Ao(M) and a < min(l,cro), 

(115) ||i?(A±i£)/||H2,-i/2-.(M) < C(M,A)A-i/2||/||^,o,i/2+.(M), 

(116) ||i?(A±ze)/||j^2,-3/2+.(M) < C(M,A)||/||^o,i/2+.(M). 
Wc start by recalling the Friedrichs inequality: 

(117) l|w||z,2(M) < C'«|suppu|i|||VM||z,2(M), 

which holds for any smooth function of compact support on M. The Friedrichs 
inequality is a direct consequence (by Holder's inequality) of the Sobolev inequality 
(with the same constant Cg): 

(118) <C^«I|V«||l2(m), 

which holds for all smooth functions of compact support on M with the constant 
Cs = Cs{M) is related to the the isoperimetric constant X{M) for M: 

I(M) = inf = X{M)l (^^)\ 

^ ' N |7V|"-i ' ^ ^ ^ \2{n-l)J 

The infimum above ranges over all open submanifolds N with compact closure and 
smooth boundary dN, and A{dN) denotes the area of dN. A proof of this Sobolev 
inequality can be found for instance in [25]. 

Let u := R{X ± ie)f, thus 

(-Am + V{x) -At i£)u = /. 

Let V = Reu then 

(119) (-Am + Vix))v = F = Re(/ + (A ± ie)u). 
We normalize 

(120) ||w||jjO,-3/2+<,(M) = 1, ||/||^0,l/2+<,(M) = 

with a small constant S. We will then show that a sufficiently small Aq and all 

|A|,e < Ao 

(121) ||u||jjo,-3/2+.(M) < ^+C{M,A)5, ||u||^o,-i/2-.(M) < C{M,A)X-^/^6. 
By Lemma 2.4 it will suffice to show that for a sufficiently large^^ compact set K 

||Vw|U2(^) + \\u\\l2^k) < {2CK)-'+CiM,A)S. 

with the constant Ck from (38). Let /i > be a small constant to be determined 
later and set ^ 



^^In what follows the set K will be fixed and the small constants Ao , <5 will be allowed to depend 
on K. 
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Define a sequence of increasing nested sets 

B(2^i?o) n {x : v{x) > f^} 

Let xe be smooth cut-off functions on M adapted to B{2'^Rq), i.e., Xf(x) = 1 for 
X e B{2^Ro) and xeix) = for a; e B%2^+'^Ro). 

Multiplying (119) by^^ x?(^ ~ and integrating by parts we obtain 

/ xl(|V(w - + V{x)u{v - iii)+)dg = -2 / xt^Xi ■ ^{v - l^e)+{v - fie)+ dg 

Jm Jm 

+ / XeF{v - t^i)+ dg 
Jm 

By Cauchy-Schwarz and a simple rearrangement, 

/ xU\^{v - + V+{x)v{v - ne)+)dg < 4 / \Vxef{v - /x^)^ dg 

Jm Jm 

+ 2 / x^\F\iv - dg 
Jm 

+ 2i^e [ xeV- ix){v - i^e)+ dg + 2 [ xeV- {x){v - l^ef+dg 
Jm Jm 

Using the smallness assumption on the negative part of the potential V we have 

/ xlV-{x){v-fxe)ldg<p(f {xi{v ~ f^e)l)^A~ 
Jm ^Jm ' 

< pc[jjx}Y^{v - + \^xA\v - W)+) dg), 

where the last line follows from the Sobolev inequality (118). Also, choosing a 
sufficiently small (universal) constant a, 

xlV-{x){v - Hi)+dg < ^nP\Ki+i\'^ ( I {xi{v - 
Jm ^Jm ' 

<a-^Cii\0'\K,^yf^ -,a(^ Jj^j\V{v - fie)+\^ + \Vxi\\v - dg). 

Similarly, by Cauchy-Schwarz, 

/ Xe\F\{v - l^e)+dg < a-'\Ki+i\i [ x}\F? dg + a\Ki+^\- i j xl{v - ,it)\dg 
Jm Jm Jm 

Thus, assuming that /3, a are sufficiently small, 
(122) 

c/ x||V(w-M«)+Pc^5< / \^Xe?'{v- IJ'i)\dg + a\Ki+i\-^ xl{v - l^if+dg 
Jm Jm Jm 



+ a-^\Ki+^\^ [ Xi\F\'dg + a-'n',f3^\K,+i\ 
Jm 



im 

for some universal constant c > 0, which depends only on the Sobolev constant Cg. 



^^Here of course we use the notation a;+ := max(a;,0). 
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By Fricdrichs' inequality (117), wc have 

Jm \Jm Jm / 

Therefore, substituting the bound for /jv^X^|V(t; — iJLi)+\'^dg from (122) and using 
smallness of a, we obtain 

c/ Xe{v - l^ef+dg < \Ki+i\i \Vxef{v - l^e)l dg 
Jm Jm 



Jm 



We can further simplify this by using the bound < CRq2 

ciRo) [ Xe{v-f^e)ldg<\Ki+i\^ [ \Vxe\''iv - fxe)l dg 
Jm Jm 

+ 2«/ x'e\F\'dg + n',P'2-'. 
Jm 

It is important however to keep the factor |if;+i|n in front of the first term! 

By the Chebyshcv inequality we have 

\Ki+l\ < J- 3- r2 / Xi+2iv - W+2) + - 

Observe that 

|Vx«| < 2-^R^'^X(+2, (w+i-/U^+2) = ^/i7(l+7)"^"^ {v-He+2)+ > {v-l^e)+- 
We then obtain 

c(i?o) / Xt{v-l^e)ldg<C{i^,7){l+l)^2-^'( [ x'+2(^ " ^+2)^+^5)'^" 
Jm ^Jm ' 

+ 2"' [ Xe\Ffdg + n^p^2^^. 
Jm 

The constant C{ii,j) is essentially (/Z7)~^/". In particular, C{fx,j) becomes large 
as /X, 7 ^ 0. 

We now define 

^, = 2-(3-^-)^ / x'Av-^e)ldg 
Jm 

so that, with the normalization (120) 

Ae < c(it!o)||u||^o,-3/2+<,(M) = c{Ro) < 1- 

Then 

A, < C{i^,j,Ro) (^(l + 7)"2(-2+*(3-2-))^^i+| +2(1+2-)^ J Xe\F\^dg + fj,^P^2<-^-^+^''^' 
Observe that for n > 3 

-2+-(3-2a)<-^, 
n 3n 
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which ahows us to conclude that with an appropriate choice of a sufficiently small 
7 = 7(cr) (explicitly, we need 61n(l + 7) < C7ln2) there exists a positive constant 
CO > such that 

Ai < C{i^, 7, Ro) ( 2-2"%'+| + 2(1+2-)^ / XelFfdg + ^-^pH^-Wi^ 

\ JM 

We now recall our normalizations (120). In particular, 

/ x'AF\'dg<[ xK\f\' + {^'+e')\um<CS' + 2^'-'^^'C{Ro){X'+e^) 
Jm Jm 

and Ae<l. Thus 
(123) 

Ai < C{ll, 7, Ro) (2-'"^^^+2 + 2(1+2-)^ + 2(3-2-)^(a2 + £2)^ + ^2^22(n-3+2a)£^ _ 

Iterating (123) we obtain that for any k <^ L 

Ak < C{ii, 7, i?o, fc)^(2-'^^' + ,52 ^ a2 + £2 + 

The constants 7 and Rq have been fixed, independently of A, /3, 5. Constant /x will be 
chosen to depend only on the compact set K. Therefore, using the smallncss of A, 5 
and /3 we can find a large integer Lq = Lo(M; 7j -^o) such that for all k < ko ^ Lq 

(e2 + A2+/32)<C'(M,7,i?o,A;)^»2-i™"-^°, C{iJ.,'y,Ro,k)^''2-'^^o < 2-ioO"-£'o. 

Thus, for all k < ko 

Ak < 2-1™"-^" +C(Lo)<5^. 
We can assume that the set K is contained in the ball B(2*^°i?). Therefore, 

f (V- ^lko)l ^ 2(3-2-)fco(2-100nLo + C{Lo)S^). 
JK 

The above inequality implies that 

f vl< An^\K\ + 2(3-2<-)feo(2-ioonLo ^ C{Lo)8'^) < {SCk)'^ + CS^, 

JK 

provided that fi is chosen to be sufficiently small relative to the size of K. Repeating 
the above arguments for the function —v = — Reu and similarly for Imu we conclude 

that 

/ |up < (2C;^)-2 + CT2_ 

Jk 

as desired. The remaining estimate for Vu is straightforward and follows immedi- 
ately by integrating the equation (—Am + V{x) — A =F = / against u over the 
set K. We omit the details. 



13. The quasimode counterexample 

We now present the (standard) counterexample in Proposition 1.14 which shows 
that the losses of exp(C\/A) which arise for instance in the unique continuation 
estimates in Section 8 are in fact sharp. For more sophisticated versions of this 
type of construction in the more general context of an elliptic trapped geodesic, see 
[23]. 
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We begin by reviewing some basic facts about spherical harmonics. Let {S^,gn) C 
R"+^ be the unit sphere with the standard metric, which we parameterize in Euler 
polar co-ordinates as 

:= {((sin 6*)^, cos 6*) : < < 7r;a; e S*""^} 

Observe that the Laplace-Bcltrami operators on and on S"~^ are related by the 
formula 

(n - 1) cos 6* 9 1 



Asn = 



56*2 sin 6* d0 sin^ 



where A^n-i is of course applied to the u) variable. Thus if Z > is an integer, and 
Yi{oj) is a spherical harmonic of order I on S^~^, i.e. 

Asn-iYi{uj) = + n- 2)Yi{uj), 

then the sectorial harmonic 

Ui{9,Lo) —sin^e) cos{e)Yi{oj) 

is a spherical harmonic on S*", in fact wc have 

(_A5„ - X)Ui = 0, where X:={1 + + n). 

We shall think of being I as being much larger than n, so X ^ P . Because of the 
sin'(^) factor, Ui will be highly concentrated near the equator = 7^/2 (which is a 
stable trapped geodesic of the sphere S"). Indeed, if we normalize Yi to have 
norm equal to 1 on S"~^, then a simple computation shows that 



I 



/S":|6l-7r/2|<7r/4 

while 

/ \Ui\^ + \VUi\^ = 0{e-^')=0{e-^^) 

JS":|e— 7r/2|>7r/4 

for some constant c = c„ > depending only on the dimension. 



We now transfer this phenomenon to the setting in Proposition 1.14. Let M C R""*"^ 

be any smooth n-dimcnsional manifold isomctrically embedded in R"^"'^, which is 
equal to the plane R" x {0} C R""*"^ outside of a compact set, and contains the 
portion of the sphere 

S^^^^ := {((sin6»)a;,cos6') : 6 > 7r/8;a; e 5""^} C 5" C R"+^ 

In other words, M is formed by gluing a large portion of the sphere (which contains 

the equator 9 7r/2) to the Euclidean space R". Note that M thus inherits the 
equator of 5" as a stable trapped geodesic. Now let I be a large integer parameter, 
set A = A; := {I + + n), and consider the "quasimode" u : M ^ C defined by 
u := Uix, where x is a smooth cutoff supported on the set S'^^^g which equals one 
on S^^i^. Then from the previous computations we see that 

II II \-l/2 



and 



-Am - A)u||^o,i/2+.(M) = 0(6-^=^), 



and thus for £ = e; sufficiently small 

II (-Am - (A ± i£))u||H0,i/2+.(M) = 0(e-^^). 
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Proposition 1.14 follows (with a small value of Co) by setting /; :— (—Am — (A ± 
ie))u. To obtain a larger value of Co, we simply scale this example, replacing the 
unit sphere 5" by a sphere of larger radius; we omit the standard details. ■ 



14. The Bessel matching counterexample 

We now present the proof of Proposition 1.21. We begin by considering the 
Helmholtz equation {H — X)u = / in polar co-ordinates (r,w) on R", thus 

^ n - 1 o> ^ 1 ^ , _v + X)u - -F 

where wc ignore the singularity at r = for now. If we now assume V to be radial, 
V = V{r), and apply the ansatz 

(124) u{r,uj) = r-("-i'/2t;(r)y;(w), F{r,u) = -r-^''+^^/^G{r)Yi{uj) 

where Z > is a large even integer parameter and Yj is a spherical harmonic of 

order I on S"~^, normalized to have L^(S'"^^) norm equal to 1, then the Helmholtz 
equation becomes the Bessel ordinary differential equation 

(125) Vrr - -Vv + Xv = G, 
where L := 1 + Also, we observe that 

/>oo 

||w||ho,-i/2-<,(m) ~ / r-'^'^\v{r)\^ dr 
Jo 

and 

POO 

l|i^llHo.i/2+.(M) - / r'^mr)]^ dr. 
Jo 

Let us temporarily ignore the contributions of the A, V, and G factors, and consider 
the ODE 

(126) Vrr - ^^^^2 = 0- 

This equation has two linearly independent solutions, and r~^+^. The former 
function decays quickly as r 0, and the latter decays quickly as r — > oo. To 
exploit this, let us define v by fiat to be a smooth function on (0, cjo) which equals 
r^ when r < 1/2, equals r~^+^ when r > 1, and is smooth and positive in between 
(of course, the smoothness bounds on v will depend on L). Then the function 
Vrr — ^^\.2^^ v is smooth and supported on the annulus {1/2 < r < 1}. Note also 
that u is then equal to r'F;(w) near the origin, which is a smooth function (indeed, 
it is a harmonic polynomial of degree I). 

Now let us define the spherically symmetric potential = Vi by 

L(L-l) 
y Vrr- r'i ^ 

V 
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By construction, V is smooth and supported on the annulus {1/2 < r < 1} (though 
the bounds on this potential will get worse as I increases), and we have 

L(L-l) 

Vrr 2 ^ — Vv = Q. 

Reversing the above steps, this means that the corresponding function u is an 
eigenfunction of the Schrodinger operator H := — A + with eigenvalue zero, thus 
A.U = Vu. Also, u decays at infinity like \u{r,w) \ = 0(r-("-i)/2--^+i) = ©(r"'""), 
and in particular will be bounded in the space iI°'^/^+'^(M) for I large enough. 

Next, we let m > 1 be an arbitrary integer and set := uXm, where Xm is a 
smooth cutoff to the region {r < 2m} which equals one when {r < m}. We also 
set \ = \m m~'/^°, and set Vm -—V + Xxm- Then we see that 

(-A -Vm- X)Um 

vanishes outside of the annulus m <r < 2m, and has magnitude ©((m"'"*"*^^^^) on 
this annulus, thus 

||(-A + Vra- A„)m„||ho,i/2+.(m) < 0,(m-'+°(i)). 
In particular for sufficiently small 

||(-A + Vm- {Xm ± i£m))«m||H0,i/2+.(M) < 0;(m-'+°W). 

Also, by construction of Um and Vm we see that (if / is sufficiently large) 

||'«m||/i-o,-i/2-<T(jv^) > C; > 0. 

and 

sup||(a;)i+"°V^y™|Uoo(R„) < oo for all a > 0, 

m 

and the claim follows (by taking I large enough depending on Cq). To complete 
the construction we need to show that the hamiltonian Hm = —A + Vm does not 
contain a small negative eigenvalue or an eigenfunction or resonance at zero. (Of 
course, the potentials Vm are very close to a fixed potential V which does have an 
eigenfunction at zero, which is indeed the cause of the bad behavior of low frequency 
limiting absorption for the perturbed potentials.) By the Cwikel-Rozcnblum bound, 
the Hamiltonian H = —A + V has only finitely many negative eigenvalues labeled 
Ai < . . . < Afe < c < 0. Their number and location is independent of m. The 
operator norm of the perturbation Hm — H is bounded by m^'/^'', which implies, 
by perturbation theory, that for sufficiently large values of m the hamiltonian Hm 
has at least k negative eigenvalues A'^ < . . . < AJ, < c + m~'/^° < 0. Moreover, 
denoting the linear span the associated eigenfunctions of Hm by Pm, we have 

/ H(P(P>0 

for any (j> € i?^(R") in the orthogonal complement of Pm- Recall that the potential 
V is supported in the annulus 1/2 < r < 1. We can also assume without loss of 
generality that |V^(a;)| < 2. As a consequence, for any <p e iI^(R") in the orthogonal 
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m 1 



complement of P„ 

>j[ + V\(l>f + XXm\<f>\') + f >^Xm\<l>\' 

(127) > (^IV^I^ + ^xJ0l^) > C-A /^^ 

This estimate already implies that Hm has precisely k negative eigenvalues X'l < 
. . . < X'f. < c + m^'/^" < and that is not an eigenvalue. On the other hand if 
(j) € na>o-ff^'~2^"(R") is a zero resonance it is formally (the corresponding eigen- 
functions are exponentially decaying) orthogonal to the subspace Pm- Moreover, 
Vcj) e na>oi?^'5-"(R"), Acj) e na>oi?°'5-"(R"), which implies that (127) holds 
and thus (p vanishes identically. ■ 

Remark 14.1. It is possible to sharpen this example a little bit by making the 
potentials compactly supported in the annulus 1/2 < r < 1. This is achieved 
by replacing the approximating ODE (126) by the minor variant Vrr — + 
(Am ± i£m)v = 0, and replacing and r~^~^ by Hankel functions (which have 
similar decay behavior at zero and at r ~ m respectively). We omit the details. 



15. Applications 



15.1. Spectral applications. We now prove Proposition 1.30 and Proposition 
1.31. 



Proof. (Proof of Proposition 1.30) Suppose first that we have an outgoing resonance 

u for some A = for some k > 0, thus u € H'^'~^^'^~'^{M), Hu = \u, and 
{drU - inu) € H'^'-^/'^+''' {M\Ko) for some cr' > and all cr > 0. By shrinking 
a, a' if necessary we can take Q < a < a' < 1/2. By elliptic regularity, u e 

Write K := VX, let < (5 <C min(l, k), and let z := k + iS. Consider the function 

where (r) := 

Since Hu = Xu, one has 

{H — z^)u = —2iSKU + 6^u, 

while from (65) one has 

{x) 

in the exterior of Kq. A straightforward application of the product rule then gives 
the equation 

(if-^V=^e-*<^>/5 
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where fs takes the form 

fs = 2{drU - iKu) + O ^■^l"!^ + ^ ((x) 
in the exterior region M\Kq, and takes the form 

fs = 0{u) + 0{\Vu\) 

in the compact region Kq; here all implied constants are allowed to depend on M 
and K. Applying the limiting absorption principle (Theorem 1.7) we conclude that 

||M5||ffO.-i/2-„(M) < C(A)J(||e"''<^>(a^U-ZKM)||^o.i/2+.(M\Xo) 
+ ||e-^<^>M||Hi,-i/2+<.(M))- 

Now observe that 

||e-*<^>u||Hi,-i/2+.(M) < C5-^^\\u\\h^,-u2-.(m) 

and 

We thus have 

||w*||/i-0,-l/2-<,(Jln) < C{X){5^~'^'' \\u\\jJl,-l/2-a (^1^)+^" "'^ II ^r^"* || JJO, - 1/2 + ^' (^f ) ) 

Taking limits as 5 and using monotone convergence we conclude that u is 
identically zero, which is absurd. 

The same argument rules out incoming resonances and eigenfunctions at any posi- 
tive energy A > 0. □ 

Remark 15.2. If the Hamiltonian H obeyed the limiting absorption principle (17) 
(note in particular that the constant here docs not blow up as A — > 0), a modification 
of above argument would also rule out an eigenvalue or resonance at zero; we omit 
the details^^. This provides a converse to Proposition 1.23. 



Proof of Proposition 1.31. The absence of singular continuous spectrum in an in- 
terval (a, h) is guaranteed by the condition that 

,.b 

sup / \lm.{R{X + ie)(t),(t))L'^(M)\^ d\ < CO 

0<€<lJa 

for any function (p E C^{M), sec [80, Theorem XIII.20]. (Indeed, the spectral 
measure associated to (p will have an density with respect to Lebesgue measure.) 
The result now easily follows from the resolvent estimates established in Proposition 
1.7. □ 



In fact, by taking z to be a number such as 2S + iS, it suffices to have the limiting absorption 
principle (17) in a sector such as {z = K + ie : < £ < k}. 
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15.3. Local smoothing estimates and integrated local energy decay. We 

now prove Propositions 1.34 and 1.38. As the arguments are standard (dating back 
to [63]; see also [16], [27], [28]), we shall be somewhat brief in our discussion. In 
particular, we will work formally (assuming that all exchanges of integrals, etc. are 
justified); one can make these formal computations rigorous by standard approxi- 
mation arguments which we omit here. 

We begin with the proof of (22) from Proposition 1.34. Henceforth we fix M, V, A 
and allow all constants to depend on these quantities. 

Prom Duhamel's formula we have 

Jo 

If we can establish the retarded estimate 

(128) II / ill/2p^e^(*-*')^F(t')||L?ffO,-i/.-.(RxAf) < C||PL2^o.v.+.(RxM) 

Jt'<t 

then by time reversal we also have the corresponding advanced estimate in which 
the constraint t' < t in (128) is replaced by t' > t. Summing, we obtain 

II / if^^^Pj/e'^*"* ^■^P(t')||j;,2^0,-l/2-<,(RxM) < C'I|-P'IIl2H0,1/2+<,(RxM); 

taking the inner product of the expression inside the left-hand norm with F and 
rearranging, we obtain the inhomogcneous estimate 

II / H'/^PHe-^''"F{t') dt'h^^M) <C\\F\\^.HO.m+.iR^My 
Jr 

which by duality gives 

||ViI^/^PH-e'*^'Uo||L|/J'0>-i/2-a(RxM) < C\\H^^'^PhUo\\l^(M) 

for any test function uq. From all these estimates we sec that to prove (22) it 
suffices to establish the retarded estimate (128) and the standard elliptic estimate 

WH'^^'^PhuoWlhm) < C||"o||ffi/2(M)- 

To prove the latter estimate, we observe from the TT* method that it is equivalent 
to H^/'^Ph mapping iI^/^(M) to H-^f^iM), which by interpolation and duality 
follows from H^^^Ph mapping H^{M) to L'^{M), which by another application of 
the TT* method (and the boundedness of Ph on L'^{M)) follows from the Dirichlet 
form {Hu,v) being bounded in H^{M). 

Now we turn to (128). By a limiting argument, it suffices to prove the damped 
retarded estimate 

(129) 

liminf II ^^^^i?^/'PHe*(*-*')(^+^^)P(t')llL?H0,-V2-.(RxM) < C||F||i2^o,i/2+.(R,,M)- 

Following Kato [63], we perform a Fourier transform in the time variable 

F{t) = [ e'^*F(A) dX. 
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The expression inside the left-hand norm of (129) can then be expressed as^^ 

e'^*H^/^PHR{X - ie)F{X) d\. 



Applying Plancherel's theorem (and Fatou's lemma) to both sides, we see that it 
suffices to establish the estimates 

limsup \\H^/^PhR{X - ie)/||//o,-i/2^<,(M) < C||/||^o,i/2+^(m) 

e->0+ 

for all A e R and all test functions /. If A is negative and bounded away from 
the origin, then H^^^PhR{X — is) is bounded on L'^{M) (thanks to the spectral 
theorem), so we may assume that A is either positive, or close to the origin. 

Using the elliptic estimate 

||F^/^w||ho,-i/2-<,(m) < C||w||hi,-i/2-.(m) 

it suffices to show that 

\\PhR{\ - i£)/||H-l,-l/2-<,(M) < C'||/||jJ0,l/2+a(M) 

whenever A e R and e is sufficiently small. 

For A sufficiently close to the origin, this follows from Proposition 1.23. For A 
positive and bounded away from the origin, the claim follows instead from Theorem 
1.15. This completes the proof of Proposition 1.34. 

Remark 15.4. It is clear from the above argument that the global-in-time local 
smoothing estimates are in fact equivalent to the limiting absorption principle (with 
the appropriate decay in A in the constants). 

Proof of Proposition 1.38. As before, we allow all constants to depend on C, M, A. 

We rewrite the wave equation in the form 

$j - ^$ = F, 

where ["Vf:= (^V and A is the matrix operator 



A- 



1 
-H 



We begin with the proof of (24). As in the proof of Proposition 1.34, it suffices to 
consider the solution of the retarded inhomogeneous problem 

$(f)= f e(*-^)^F(s) ds. 

Js<t 

Applying Plancherel as before, we see that bounds on the retarded solution are 
equivalent to bounds on the resolvent {iA — /x — iO~)~^. More precisely, if ^ = 
(V'i;V'2) is a solution of the equation 

{iA - ^i + ie)'i' = {0,g) 



^^One can also shift the contour, taking advantage of the projection away from any negative 
eigenvalues to work with resolvents in a neighbourhood of the positive real axis, if desired. 
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then we need to show that uniformly in e > and G R that 

||VV'l||_f/0,-l/2-^(R,„) + ||'i/'2||ffO,-l/2-<,(R_„) < C'||5||^0,l/2+<,(R,„). 

Since 

(130) M - M + ^e)-' = R{{i^ - lef) ( ^'.'^^ ^ ^ ) 

we have 

and the desired estimate (24) follows immediately from Theorem 1.7. 

The proof of (25) follows from a similar argument using Lemma 2.4; we omit the 
details. □ 



15.5. The RAGE theorem. We now prove Proposition 1.32. We establish the 
claim just for the Schrodinger equation; the claim for the wave equation is analogous 
and is left to the reader. 

By stability of the Schrodinger equation it suffices to establish (20) on a dense 
subset of iy^(M). Thus we may assume that u{0) e H^{M). Moreover, Proposition 
1.30 together with the assumption that 'u(O) = / is orthogonal to all eigenfunction 
of H implies continuity of the measure dfifiX) = {dE\{H)f,f) defined from the 
spectral measure dEx{H) associated to H. As a consequence, by a density argu- 
ment, we can assume that u(0) is spectrally supported on an interval A e (a, h) 
with 0<a<6<ooof energies A. For such functions we have both the limiting 
absorption principle and the local smoothing estimate 



f 

Jo 



||w(i)||^2,-i/2-.(M) dt <oo 

•) u 

for some cr > 0. 

Let £ > 0. Applying the monotone convergence theorem, we can find Tq > such 
that 

/■OO 

/ II"(*)IIh2,-i/2-.(m) 
Now for any T > To + 1 and g e i?°'V2+<"(M) we have 

/ «(r,x)<?(a;) = / / ^Xu{t,x){t-T^\)g{x)) dxdt 
JR" Jt-1 Jm "I 



{iHu{t, x){t-T + l)g{x) + u{t, x)g{x)) dxdt 

T-l JR" 



< C||ff||H0,i/2+<,(M) ij^ Ju{t)\\%2.-i 



/2-<7(M) dt 



< C'£||5o||ho,i/2+^(m)- 

This implies that u{t) -J> in H°'~^/^~'^ (M) for some cr > 1/2, and the result 
follows. 
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15.6. The limiting amplitude principle. Wc now prove Proposition 1.33. We 
will follow an approach of Eidus [41] and deduce the limiting amplitude principle 
from the limiting absorption principle, combined with a Holder continuity property 
of the resolvent. More precisely, we will use the following fact: 

Proposition 15.7 (Holder continuity). Let the assumptions be as in Theorem 1.1. 
Let K he a compact subset of M , and let F be a compact subset of the right half- 
plane {z : Rez > 0} that avoids 0. Then there exists C > and a > such that 
one has the Holder continuity bound 

\\Riz)f - R{z')f\UiK) < C\z - z'rWfh^K) 

for all z,z' G F with Imz, Imz' > 0, and all f e L'^{K). 



Let us assume this proposition for the moment and prove Proposition 1.33. For 
simplicity of notation we will assume that H has no eigenfunctions; the general case 
is analogous and proceeds by inserting a projection Ph to the absolutely continuous 
portion of the spectrum of H (which commutes with all spectral multipliers of H , 
or of related operators such as the matrix operator A introduced below) throughout 
the argument. By Proposition 1.32 and linearity, it suffices to handle the homo- 
geneous case uq = ui = 0. As in the proof of Proposition 1.38, we introduce the 



vector $ = ( ^ ) and rewrite the wave equation in the form 



1 anrl A no fVio i r 





where F := [ . ] and A is the matrix operator 



1 
-H 

\llJ,V 



A:= 

It will sufiice to show that 

in Hi{K) X L'^{K) as t +oo. 
By the Duhamel formula, we have 

e-v*$(t) = f e^^-^^^^-'^'UsF 
Jo 

ft 

lim / e'^'-'^'^-^^'^'-''>dsF 







^ ^t{A-in) 



{A-iiJ.- 0+)-^F -{A-ifi- 0+)-^F, 



writing {A — i^i — 0+) ^ as the weak limit of the {A — ijj, — e) ^ as e — ^ 0+. Prom 
(130) we have 

-{A 0-)-'F = Rii^^ - ^0+f) (^^^^) 

and so from Proposition 1.12 (and elliptic regularity) one has 



(^-iM-0+)-F^(.;^) 
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in. (K) X (K) . Discarding the bounded phase e , it thus remains to estabUsh 
that 



e 



tA 



converges strongly in H^{K) x L'^{K) to zero as t — +00. 

Let C be a semicircular contour of the form 

{iy : |y - mI < ^} U + re*^ : 7r/2 < 6* < 37r/2} 

for some radius < r < (e.g. one can select r := /u/2). From the Cauchy integral 
formula and the spectral theorem (using the fact that the spectrum of A is the 
imaginary axis) we have 

e^-^iA-in-O+y^F = lim lim J- [ — (A-z-e')-^fdz 

+ e'^l^^-^^,l>riA-^fl-0+)-^F 

strongly in H^{K) x L'^{K). From the spectral theorem, l\A-iiJ,\>r{A — ifJ, — 0~^)~^F 
lies in H^(M) x L'^{M), and so by the RAGE theorem, the second term on the 
right-hand side goes to zero strongly in H^{K) x L'^{K) ast ^ ±00. It thus suffices 
to show that 

lim lim lim / -. (A-z-e')-^fdz = 

Strongly in H'^{K) x L'^{K). 

Prom the Cauchy integral formula we have 

dz^Q 



I (J z — 111 + e 
for £ small enough, and so it suffices to show that 



(131) lim lim lim / Ge'(z)dz = 

t-s-+cx) £^-0+ e'->-0+ Jc z — in — e 

strongly in H'^{K) x L'^{K), where 

G,,{z) := iA-z-s')-^F-iA-in-s')-^F. 

Prom (130) and the Holder continuity property in Proposition 15.7, G^' is Holder 
continuous in H^{K) x L'^{K), uniformly in e'; in particular, as G'^ vanishes at i/i, 
we have the bound 

(132) \\Ge'{z)\\H^K)>,LHK) = 0{\z - 

for some ^ > 0. This is already enough, when combined with Minkowski's integral 

inequality and the dominated convergence theorem, to control the semicircular 
portion of the contour C; it remains to demonstrate that 

lim lim lim / Gr'(iu + iv) dy = 

e->-0+ s'^0+ J_r ly - e 

strongly in H^{K) x L'^{K). Prom (132) and Minkowski's integral inequality we 
have 

-G,,{in + iy) dy =0{s^) 



i: 



ly-s 
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for any < s < r, so it suffices to show that 

lim lim lim / Ggiiiu + iy) dy = Q 

t^+oos^O+ e'^0+ Ji ty - e 

for any compact interval / in [— r. r] avoiding the origin. However, from (132) we 
see that for any unit vector w in H^(K) x L'^(K), the scalar function 

y H> (— — Ge'{ifi + iy),w\ 

is uniformly Holder continuous in e, e' on /, and in particular (by the Arzela-Ascoli 
theorem) is precompact in the uniform norm. Prom the Riemann-Lebesgue lemma 
we thus have 

lim lim lim // t- Gs'Ufj, + iy) dy,w\ =0 

t^.+cx> e^o+ e'^0+ \ J i ly - S / hHK)xLHK) 



uniformly in w, and the claim follows by duality. 



It remains to prove Proposition 15.7. We first observe from the limiting absorption 
principle that 

\\R{z)f\\H^^K) < C\\f\\L.^K) 

uniformly for all z in a compact set avoiding the origin. Thus, by interpolation, 
it suffices to obtain an LP'i^K) — > L'^{K) Holder continuity bound, thus we need to 
show that 

\\R{z)f - R{z')f\\L^^K) < c\z - z'nnWiK) 

for all z,z' in F with positive imaginary part, and some sufficiently small o" > 
(which may be different from the one in Proposition 15.7). Clearly, we may assume 
that \z — z'\ is smaller than any given absolute constant. 

Let cr > be a small exponent to be chosen later. By the triangle inequality, it 
suffices to prove this claim under the additional assumption that Im(z') > 0.1|z — z'| 
(say). Write ry := Im(z'), then 77 > is small and \z — z'\ = 0{ri), with \z\, \z'\ 
comparable to 1, and with Re(2:), Re(z') positive and bounded away from the origin, 
and our task is to now show that 

\\R{z)f - R{z')fU.^K) <Cr^''\\f\W(.Ky 
Using the resolvent identity 

R{z) - R{z') = {z- z')R{z)R{z') 
it thus suffices to show that 

\\R{z)R{z')f\\L2(K) < Cr?-i+ni/|U=(x), 
so by duality it suffices to show that 

\{R{z')f,R{z)h)L.^M)\<Cr^^'' 
whenever /, /i e L'^{K) have unit norm. 

Fix /, h. Recall that z' has imaginary part ry. This effectively localises R{z')f to 
the region {x : {x) = 0{l/rj)}. Indeed, if we write z' := — w^, where Imw > 0, then 
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|a;| is comparable to 1 and Imw is comparable to rj (with constants depending on 
the compact region F). Using the elementary identity 



1 



e cos tx dt 



we see that 

1 r°° , 

R{z')f = - / e"'"* cos{t^/^)f dt. 
^ Jo 

The wave propagators cos{t\/—H) are contractions on Lp'[M), and cos{t\/—H)f is 
supported in the region where (x) < t + 0{1). Prom the exponential decay in the 



(133) / {xy°°\R{z')f{x)\'' dg < Cri'°° 

J|z|>r,-i+<' 



(say), for any fixed a > (allowing C to depend on a). Thus, it will suffice to show 
that 



(134) I / i;iR{z')f)R{z)h dg\ < Cr^' 

JM 



-1+a 



for a smooth cutoff t/j to the region (x) < 2ri^^'^ that equals 1 when (x) < rj^~'^ , 
as the error term caused by 1 — ^ can be estimated by (133) and the limiting 
absorption principle (with plenty of powers of r] to spare). 

Write u := R{z')f and v := R(zh). From the limiting absorption principle we have 

(135) ||m||//o,-i/2-^(m), \\v\\ho.-i/^—{m) < C. 

If we apply these bounds and Cauchy-Schwarz to estimate (133) directly, we obtain 
a bound of 0{r]~^~'^'^), which barely fails to be adequate for our purposes. To 
obtain the additional powers of 77 needed to close the argument, we take advantage 
of the fact that z' and z lie on different sides of the real axis, and so R{z')f and 
R{z)h obey opposing radiation conditions. Indeed, if we write u := R{z')f and 
V := R{zh), then from Proposition 1.12 we have 

(136) \\{dr - iz')u\\HO,-i/2+3.(^M\Ko) < ^ 

and 

(137) \\{dr + iz)v\\HO,-i/2+3.^M\Ko) < ^ 

(say), if a is small enough. To use these facts, denote the integral in (134) by /. 
Prom (135) and Cauchy-Schwarz, we can write 



/ {tp-'il)o)uv dg + 0{l), 



where t/jq is a cutoff to a fixed compact region that equals 1 on Kq. Applying (136) 
(and using (135) to estimate the error), we can then write 



= {tp - '4'o){dru)v dg + 0{r]~ 

I'Z Jm 
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We integrate by parts to then obtain 

Jm 

+ 0{[ -l-i;\u\\v\ dg). 

Using (135) and Cauchy-Schwarz, the final error term is 0(r?~^+'^) for a small 
enough. For the main term, we use (137), estimating the error again using (135) 
and Cauchy-Schwarz, to obtain 

^ = -A/ (V- - V-o)"!^^ ^5 + 0(r?-i+") 
iz' Jm 

which simplifies (using (135) and Cauchy-Schwarz one last time to remove the '^o 
cutoff) to 

I = --,! + 0{n-'+n. 

z 

As z,z' are both in the upper half-plane, have magnitude comparable to 1, and 
are within r\ of each other, we see that — p- is a bounded distance away from 1 for 
T] small enough. Wc thus conclude that / = 0(ry~^+'^) as required. This proves 
Proposition 15.7 and thus Proposition 1.33. 

Remark 15.8. It should be clear from the above argument that one also has a similar 
limiting absorption principle for the Schrodinger evolution (replacing H^{K) x 
L'^{K) by L'^{K)). Wc leave the details to the interested reader. It should also be 
clear from the argument that one can strengthen the convergence in H^{K) x L^{K) 
to convergence in H^~^{M) x H^~^{M) for some sufficiently large s, and dually 
one can also relax the hypothesis / e L'^{K) to / e i?°'*(M). 

Remark 15.9. The hypothesis in Proposition 15.7 that F lie in the right half-plane 
and avoid zero was needed in order to invoke the limiting absorption principle. If 
one assumes that H has no eigenfunctions (or if / is assumed to be orthogonal to 
such eigcnfimctions) , then one can extend this proposition to the left half-plane as 
well (and indeed the claim follows easily from the spectral theorem in that caase), 
and if H has no eigenfunction or resonance at zero, then one no longer needs to 
avoid the origin (thanks to Proposition 1.23). It is likely that one can upgrade 
the Holder continuity bound to a stronger bound, such as difi'erentiability, under 
further regularity and decay hypotheses on the metric g and potential V, but we 
will not pursue this matter here. 

Remark 15.10. Wc sketch here an alternate approach to the limiting amplitude 
principle that does uses wave equation energy estimates instead of Holder continuity 
properties of the resolvent, but requires the uo parameter to be large (in particular, 
potential needs to be strongly short-range), and also requires H to have no bound 
states. The main task, as noted above, is to establish decay of e*-^(^— i/x— 0+)~^F in 

H\K)xL^{K) ast^+oo. If the data (^^^ = (yt-i/i-0+)"i^^ had finite total 

energy, then this would follow from the RAGE theorem; however, from the spectral 
theorem we see that we do not expect this data to have finite energy. However, one 
can use Lemma 2.4 (or more precisely, a variant of this lemma) to show (if ctq is 

large enough) that ( J has finite incoming energy, in the sense that r~^'VuV and 



78 



IGOR RODNIANSKI AND TERENCE TAO 



Vr + ifJ,v lie in Lp'{M). It turns out that these types of estimates, together with wave 
equation energy estimates formed by contracting the stress-energy tensor against 
a weU-chosen vector field (essentially outgoing vector field dt+dr), shows that the 
energy of e*-^(^ — i/i — 0+)~^F on a forward light cone is bounded, which by energy 
estimates implies that e*-^(>l — ijj, — 0+)~^F stays bounded in {K) x L'^{K). By 
truncating away a compactly supported component of {A — i/J. — Q~^)~^ F (whose 
contribution decays by the RAGE theorem), one can upgrade this boundedness to 
decay. 

16. Decay estimates for the time-dependent Schrodinger equation 

In this section we prove Proposition 1.40. We may assume that t > 1. since the 
case < t < 1 follows from Sobolcv embedding and H^(M) energy estimates. Our 
initial analysis will be valid in all dimensions three or greater, but we will eventually 
specialize to the three-dimensional case for sake of concreteness. 

Let V be a solution of the time-dependent Schrodinger equation 

idtiJ + AmV' = 0, 
V'|t=o = V'o- 

Consider the following second order self-adjoint operator 

P := t^AM - itxirdr + ^ + ^0)- ^x^' 

with a smooth cut-off function x(r) supported in the region {x) > vq and equal to 

1 for r > 2ro for some sufficiently large tq. As we shall see later, this operator can 
essentially be viewed as a conjugate of ^^Am, via the heuristic^"^ 

(138) P « e'^'-'^^H^AMC-'^'^'/^K 

We will establish the decay estimate by commuting the Schrodinger equation with 
P. 

Recall that in the region (x) > rg the metric g on M takes the form 

g = dr^ + r^h[r]abdu}'' dJ' , h[r\ab = hab{tc) + r^'^"'' eab{r,uj) 
and 9 = \h[r]°'^drh[r]ab- By hypothesis, we are assuming that 
\{rdrf{yi)h\<Cka, k<3, \a\<2. 
Let (j) := Ptjj, then ^ solves the forced Schrodinger equation 

idtct) + AM(t> = F 

where 

F:= {idtP+[AM,P])^. 



In the EucUdcan case H = Aj^n , setting x equal to 1, P is the Laplacian conjugated by the 
pseudoconformal transformation u{t,x) i-> ^2/2 f)' which explains why we expect P 

to approximately commute with the Schrodinger operator idt + ^M- 
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To obtain bounds for F, we compute 

Tl T 

idtP = 2itAM + xirdr +2+2^^ 

and 

[Am,P] = -it[AM, xirdr + I + ^^)] - \[AM,X^r^] = -itiAMX){rdr + | + 

- 2x'itdr - 2x'itrdf - x'itdrire) - 2x'it{rdr + | + - xit[AM,rdr + 

- i(AMx')r2 - 2rxx " Xx'r^dr - ^x^[AM,r^] 

with the understanding that all expressions involving x or its derivatives vanish on 
Ko. 

In polar coordinates (outside of Kq), the Laplace-Beltrami operator has the follow- 
ing representation 



Am = 5^ + (^^^-j^ + 0^ dr + ^Ah[r]- 



As a consequence, we may compute commutators: 



[AM,rdr] := 2Am - (r^^^ + e)dr + -driA^r]), 



r 

„2l 



[Am, := 4:rdr + 2 + 2(n - 1) + 2^r, 

[Am, re] := 2{0 + r9,^)9, + - Va,^V<. + (^^—^ +e]0 + 29^^ + r(AM^). 

r V ^ / 

Putting all these estimates together, we obtain the pointwise bound 

\F\<C{{l + t)C\dr^P\+tC\d'riA 

+ r-2-2-ot (|v2 v,| + IV.^M) + (1 + t)r-2-2-o 1^1) , 

where is a smooth cut-off supported in the region ro < (x) < 2ro. 

Fix T > and let tit be the characteristic function of the interval [1,^]. Then 

(idt + Am)vt(I) = VtF + mi)6{t - 1) - mS{t - T)) 
By the global in time local smoothing estimate for cr > 1/2, 

ll»?T?!'|L2H0,-i/2-<.(M) < C (||??T-F||i,2jjo.i/2+a(M) + ||?!'(l)||l,2(M) + U{T)\\l^M)) , 

which means that 

ll?!'IUfj ,j,,H0,-i/2-.(M) < C (|1F|1l2^ ^jH0,i/2+<,(M) + ||'?^(l)||l,2(M) + \\(l>iT)\\L^M)) ■ 

We also have the standard estimate, which implies that for any t G ['^,T], 

\\4>it)\\L^M) < C (^||F||i2^ _^^^o.i/2+.(M) + ||0(1)||l2(M)) ■ 

Adding the two estimates above we obtain the standard bound 
(139) 

ll<?^(i)IU^(M) + Uhf^^^^HO'-y^-'iM) < C (ll-P'llLfj .^,H<'.V2+<.(M) + 11^(1) I|l2(M)) • 
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To relate (j) = -P^ back to ^, we now work on developing the heuristic (138). We 
compute 



A 



M e 



A 



We expand the final term Am ''^ « ^ : 

Am (e-'>^^) = -^a, (e-^'^5^(2rx + xV^)) - ^( 



i ,71— 1 



At 



-»X4t _ 



4t 



(2rx + x!r^) + 4xV + 2x + x'V^ 



n — 1 

+(^+^)(2^X + xV2) 
We thus obtain the pointwise estimate 

(j> - e^^i t'^AM (e-'>^5^ V) I < C{tC\drij\ + (1 + t)C|V'l)- 
Inserting this bound into (139) and letting 

u := e~^-'^^ijj 

we obtain 

||f^AMM||L2(M) + l|i^AMu||L2^ ,^j/fO,-i/2-^(M) < C* (||-F||L2^_^j/f-o,i/2+<.(M) 

+ t \\C\dr^\ + (1 + i)CI^I IUf,_^,H0,-i/2-<.(M) + ||0(1)||l2(M)) 

We now observe that 

IV'I = \u\, IVcu^i-l = |Vu|, \drll)\ < \drU\ + j\u\, 



IV^VI = \ylu\, < \d'rU\ + -\drU\ + ^\U\'. 



Therefore, 



\F\+tC\dr^\+tC\^\<Ct {C\dru\+r 



-2 — 2(To \yj2 



+ C\dru\ + |v^i,| + |y|) 



We thus have 

J.2 



II* Amw||z,2(m) + \\t Amw||z,2^ .^jH0,-i/2-<,(M) 

< C'^||ta^u||j^2^^^jjo,-i/2-<,(M) + ||i9rW||j:,2^ ^^jyO,-l/2-<,(M) 

+ ||ir~ V„u||j;,2^^^^^0,-l/2-2<,o + <'(M) + ll*''~^^w'"llz,fj .j,jH''.l/2-2<'0 + <'(M) 

+ ll*'"llz,fj .j,,HO,-3/2-2<,o+<'(M) + II^!'(1)I|l2(M))- 

We can couple this bound with the local smoothing estimate for the original solution 
ip, which implies (as u has the same magnitude as ip) that 

II^^IU2(M) + ll'"llz,2^ .j,jfi-o.-i/2-<.(M) < C'l|w(l)||l,2(M). 
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By combining these two estimates, we arrive at the estimate 
(140) 

||t^AM'u|U2(M) + Ili^AMMllLSj ^jH" -i/2-<.(M) 

+ ll"(l)llL2(Af) + !|(/>(1)||l2(M)) 

The first two terms on the right hand side require special care. We introduce a 
radius parameter R and spUt 

+ ||tr~^V^u|ji2^ ^^^0,l/2-2„o+-(M-) 

< CiR)\\tr-^^lu\\L2^ ^^H0-^;2-.^M) 

te[i,T] 

For ao > (T > 1/2 the second term can be absorbed by the ||t^AAf?illL2(M) in the 
left-hand side of (140), while the first term, via Bochner identities, can be partially 
absorbed by the second term^^ on the left hand side of (140) with the remaining 
residual term of the same form as ||tu||i,2^ ^^ho,-3/2-2,,o+'{m)- 

Therefore, we have that 
(141) 

+ I|w(i)I1l2(m) + II'/'(i)IIl2(m)) 

In dimension n = 3 we now proceed as follows. Using a local smoothing estimate 
for u we get the bound 

\\t'^AMU{t)\\ L'^{M) < Ct{\\u{l)\\i,-2(M) + ||?!'(1)|U2(M)), 

which gives the preliminary bound 

(142) \\AMu{t)\\mM) < y(||M(l)||L2(M) + mi)W(M)) 

for alH > 1. By the Sobolev estimate, 

l|w(i)lk~(M) < C||AM«(t)||i(M)II^^WIli2(M) < |-(II«(1)IU^(M) + ||<^(1)|U2(M)) 



Note that Bochner identity and our assumptions on the metric g imply that 
I|V^w|Il2(m) + ll(a^>"^|V«|9||i2(M) < C\\Amu\\^2. 

The dependence of constant C on the manifold M is imphcit as this inequality proved with the 
help of compactness arguments. This bound allows us to prove Sobolev multiplicative inequalities 
of the form 

l|w||z,P(M) <C||AM«|l22(M)ll«lli2(M)> 2a = (|--), 0<a<l, (a,p) 7^ (1, oo). 
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Therefore, 

P^ILfi ^,_ff"-3/2-2-o + -(M) < r^||tu||L~^,L~(M) < Ti{\\u{l)\\L2^M) + HWWl^M))- 

We can thus bootstrap the prehminary bound (142) to the improved bound 

||AM«(i)|U2(M) < 5-(||«(1)||l2(M) + \m)\\LHM)) 

for any t>l. 

We now iterate this process. If on the n**^ step the decay rate of ||Amw||l2(m) is 
then 

3 1 

Oin+l = -^OLji + 2 

and thus a„ 2. It then follows that 

IIV'WIUoo(M) < Cet-^+^(||V(l)|U^(M) + ||</'(1)||l=(M)) 

for any e > and t>\. 

In higher dimensions the norm ||to(f)||jjo,-3/2-2<7o+<T(M) is -t'^ integrable and can be 
dealt with directly in one step by interpolating between || AMw||i2(jv^) and ||«||z,2(m)- 
We omit the details. 

17. Decay estimates for the wave equation 

In this section we prove^^ Proposition 1.43. For t<2 the claim is easily established 
from H"^ X energy estimates and Sobolev embedding, so we shall limit ourselves 
to the case when i > 2. 

Let u be a solution of the wave equation 

□mu := utt - Am« = 0, 

u|t=0 = Wo, '«t|t=0 = Ml- 

In Section 4 we defined the energy-momentum tensor Qmfc associated with the 
Helmholtz equation (iJ — z^^u = 0. Here we use its spacetime counterpart^*^ 

Qa/3 := daudpu- ^gafiS'"'9f,ud„U. 

As is well known, the energy-momentum tensor Qq,^ is divergence free: 

D/^Qap = 0. 

Let K = K'^da be an arbitrary (smooth) vcctorfield. We form the quantity 

Pa ■■= QapK^ 

^^Thc argument below is an optimal but lias been retained for illustrative purposes. The 
following decay rates can be improved via an alternative approach as in [35]. 

■"^Greek indices a, /3 = 0, . . . , n with index corresponding to the t coordinate. Operations 
of raising and lowering of indices are done with respect to the space-time metric g = —dt^ + g. 
Finally, D will denote the Levi-Civita connection of metric g. 
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with the property that 

where TTap '■= ^jC-KQafi is the deformation tensor of K. We note two more identities: 

D°'{dauu) = d°'udaU. 
Then choosing additional smooth function a" and 6, we have that 

(acM^ + bdauu + Qa^K^) = D'^UaU^ + (2ac« + dab) + bd^ud^u + Qa/STr"'' 
or more exphcitly 

{aaU^ + bdauu + QaffK'^) = D'^a^u^ + (2a„ +dab)d°'uu 

+ (^TT^^ + (6 - ^tr7r)g«'5^ daudpu, 

where trTr = g'"^7r^i/. 

It is clear that one of the choices of a", 6 to simplify the right hand side is 

1 

a„ = - 

whence 

Z?" (^-^dabu'' + bdaUU + Qia0K^^ = □m&u^+ (^7r«^ + (6 - ^tr7r)g"'5^ daud^u 

Integrating this expression in the space-time slab [t, 0] x M we obtain 
(143) 

/ (-l-dtbu'^ + bdtUU + QopK^) dg- [ ( -^-dtbu^ + bdtuu + QopK^ 

Jmx{0} \ ^ / JMx{t} \ ^ 

= J (oMbu'' + + (6 - ^tr7r)g«'5^ daudffU^ dgdt 

We now make the remaining choices of the vector field K and function b. For K 
we choose a modification of the Morawetz vectorfield (t^ + r^)9( + 2trdr, known to 
play an important role in the study of the decay properties of solutions of the wave 
equation in Minkowski space and obstacle problems (see [74]). Define 

K:= {t''+xr^)dt + 2trxdr 

with a cut-off function x of the previous section siipportcd in (x) < ?'o and equal 
to 1 for (x) > 2ro- The deformation tensor of K can be computed as follows: 

TToo = -2t, TTrr = 2tx + 2trx' , 

T^ab = 2tx{h[r]ab + r9ab), T^or = -x'r^, 

tvK = 2t + 2tx{n - 1) + 2txr0 + 2trx' 

We then choose 

b={n-l)t + 2txre 
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SO that 

(144) (n^P + (6- itr7r)g"^) VauV^u = 2tr-\0abVZuVtu+O{O{l+t)\Vu\\ 

where C is a smooth cut-off supported in ro < (a;) < 2ro- Furthermore, 

(145) nMb = 2tAg{xr0) = t 0{r-'^-'^'"') 

We now analyze the expression (^—^dtbu^ + bdtuu + QopK^). 
Lemma 17.1. Let n > 3. Then for sufficiently large ro and all t^ro 
(146) 

/ (-l-dtbu^ + bdtUU + Qo0K^] dg>C{ro)t^ [ {\dtuf + \Vu\'' + u'') dg 
+ C(ro) / {{t + rf\{dt+dr)u\^ + {t-rf\{dt-dr)u\^) dg 
+ C(ro) / {{t^ + r^)\r-^Vu.u\^ + (1 + t\-^y) dg 



C(ro) / 0{r-^^'^){\dtu\^ + \Vu\^) dg. 



Remark 17.2. Similar arguments can be found in [64] for the wave equation in 
Minkowski space and [33] for the wave equation on Schwarzschild background. 



Proof. First observe that with our choices of b, K°' 
Q := —dtbu'' + bdtuu + Qo/jK^ = ^{t^ + xr^Xldtuf + |Vw|^) + 2trxdtudrU 

+ t{n-l + 2xr9)dtuu - i(n - 1 + 2xr9)u'^ . 
On the complement of the set M2ro can thus bound Q as follows. 

(147) Q> ^(|aH' + |Vu|^)-^7(ro)w^ onMt^ 

provided that t^ ro. Let x be a cut-off function supported in M2ro ^^^^ equal to 
1 on Miro ■ To prove the desired result it will be sufficient to show that 

/ X \\{t' + r''){\dtu\^ + \Vu\l) + 2trdtudrU + t(n - 1 + 2re)dtuu- ^(n - 1 + 2r9)i 
Jm \^ 2 

>C(ro) / {{t + rf\{dt+dr)uf + {t-rf\{dt-dr)uf) dg 
+ C{ro) ( {{t' + r'')\r-^y^u\^ + (1 + 1^-^^) dg 



C{ro) [ 0{r-^-"){\dtuf + \S7uf)dg 



Note that the t^v?' term on M^^^ can be obtained via a Poincare inequality from 
the estimate above and (147). 
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Define 

Su = {tdt + rdr)u, 5 = {tdr + rdtu). 

Then 

l{e + r-'){\dtu\^ + \Vu\l) + 2trd,udrU=\{{Suf + {S^^^^ 
Furthermore, 

1 t 
tdtuu = Suu — -rdr{u)'^, tdfuu = -Suu — —dr{u^). 

As a consequence, 

t{n-l) I x9tuu={n-l)l xSuu+^^^—^j {x{n + rO) + rj(!) , 
Jm Jm ^ JM 

r f f — I r y2 

t{n-l) xdtuu^in^l)- xSuu+—-e {x{n - 2 + rO) + rx!) —. 
Jm ^ Jm ^ Jm f 

It is important to note that x' > and (n — 2) » r9 on the support of x provided 
that ro is sufficiently large. It is not difficult to show that for n > 3 one can find 
constants A, B such that 

A + B = n-l, A^+B'^ <{n- ^)(n - 3) + ^(n - 1)A. 

(in dimension n = 3 the choice oiA = 3/2, B = 1/2 is sufficient). As a consequence, 



X (* + ^ WtuY + |Vu|^) + 2trdtudrU + t{n-l + 2r6)dtuu- ■^{n-l + 2re)u' 

+ ^ ^ X (((n - 2)(n - 3) + (n - 1)A - - + ^i'i?) ^5 

xO(r-2-'')(u2+t|a,u||u|) dg. 



IM 

Running the argument again we obtain that 

1 



X ( -(r + r^){\dtuY + iVwIg) + 2trdtudrU + t{n-l + 2re)dtuu- -(n - 1 + 2re)u' 



2 

> c 



;/ x({Suf + {Suf + {t' + r^)\T-^Vu^u\l + {l + %)u'\dg^( xO{r-^'"'){dtu) 
Jm \ J Jm 

The result now follows immediately from the identity 

(i + rf{{dr + dr)uf + {t- rf{{dt - dr)uf = {{Suf + {Suf) 

□ 
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We introduce notations for the energy associated with the vectorfield K and a 
conserved energy associated with the vectorfield d/dt: 

(148) EK{t) = Ek,u{T) ■■= I (|9tu|' + |V«p + u") dg 

(149) +/ {{t + rf\{dt + dr)u\^ + {t-rf\{dt-dr)u\^) dg 

(150) + / {{t" +r'^)\r-^V^u\^ + {l + t'r-^)u') dg 
E{t)^E^{T) / {\dtu\'' + \^u\'')dg. 



M 



Combining Lemma 17.1 and (143), we derive the following 

Proposition 17.3. Let u he a solution of DgU = and let T be sufficiently large. 
Then there exists a sufficiently large R <^T and a smooth cut-off function C, sup- 
ported in the region \ Mr = {x : r^ < {x) < R} such that 



Ek{T) < Giro) (^Ek{0) + ^ + ^''^"1' '^^dt'^ . 



Proof Lemma 17.1 together with (143), (144) and (145) imply that all T > ro > 1 
EK{T)<C{ro){EK{0)+E{T)) 

+ C{ro) r [ (0(r-2-o) t {\r-'V^u\l + r-^) + 0{Q{1 + t)\Wu\'') dgdt. 

Jo JMrg 

Since the expression for the energy Exit) contains both t'^\r~^Vuu\'^ and t^r~'^u'^ 
we can find a sufficiently large R such that the above inequality can be simplified 
to 

Ek{T) < C{ro){EK{0) + E{Q)) + C(ro) / {1 + t) [ C|Vwp dgdt 

Jo Jm 

with ( supported in Mrg \ Mr. Note that to get to the last inequality we also used 
a Poincare inequality to convert the term into |Vw|g, and conservation of the 
energy E{t). To conclude the proof of the proposition it remains to observe that 
the energy Ek{^) easily dominates -B(O). □ 



Henceforth we specialize to the three-dimensional case n = 3. 

The integrated local energy decay estimate of Proposition 1.38 gives the bound 

/ {1+t) [ Cl'^uf dgdt < {1 + T) [ [ Cl'^uf dgdt<C{l + T)E{0). 
Jo Jm Jo Jm 

This immediately implies that 

Ek{T)<C{1 + T)Ek{0). 
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In particular, as can be seen from (150), inequality above leads to the 1/i decay of 
the local energy. Precisely, for any compact set -B C M 

(151) X^^^*"^' + - ^^^W- 

Applying in addition the same argument to the function v = dtU we then have 



IB - t 

Prom the Gagliardo-Nirenberg inequality we thus obtain uniform decay for u: 
for any e > 0. 

This however is not an optimal result (in particular, as suggested by our experience 
with the wave equation in Minkowski space in dimension n = 3). To improve on 
it one observes that Proposition 17.3 can be iterated, [33]. Heuristically (for more 
details see [33]), the argument is as follows. First partition the time interval [0,T] 
dyadically into subintervals of the form [t, Ft] with some F > 1 sufficiently close to 
1. By finite speed of propagation, on the set [t, Tt] x M^, solution u is completely 
determined by its values at time t on the set -^^+(7(r-i)t fo'" some universal constant 
C dependent on the manifold M. In particular, F can be chosen in such a way that 
R + C{T-l)t < l/2t. Now, the estimate Exit) < C{l+t)EKi'0) in fact gives more 
than just the local energy decay. One can easily show that (151) holds for any set 
B such that B C M^^^^ with a constant C{B) independent of B (in particular of 
t). In view of the above discussion the integrated local energy decay estimate for u 
on the interval [t, Tt] should take the form 

r [ aidtu)' + \Vu\l)dgdt<C [ {{dtuf + \Vu\l)dg<^EK{0). 

Adding these estimates over all such time subintervals we obtain that 
^^(1 +i) j aidtuf + \Vu\l) dgdt < logTEKiO). 

Using Proposition 17.3 now implies that 

EK{T)<C{ro)logTEK{0) 

and 

l^mum' + \Vu{t)\l) dg < 9^^^Ek{0) 

for all t <T. Repeating the whole argument one more time eliminates the remaining 
logt term and shows that 

Ek{T) < C{ro)EKiO) 

Finally, combining the bound for u with the similar bound for v = dtU and using 
Sobolev inequalities we arrive at the desired estimate 

\Ht)\\L^ < Cet-'+' {EkA^) + EK,a,um ■ 
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